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INTRODUCTION 

Bridges  of  the  United  States  are  under  constant  attack  from  the  environment.  It 
may  be  sudden  yet  violent  effects  of  floods  or  earthquakes  or  the  long-term  effects 
associated  with  sunlight,  rain,  deicing  salts  and  freeze/thaw  cycles. 

The  quality  of  the  environment  has  always  been  of  concern  to  the  general 
populace.  Concerns  over  the  potential  pollution  of  the  environment,  by  lead  and  other 
heavy  metals,  during  the  removal  of  older  paint  systems  from  bridges  has  resulted  in 
regulations  severely  limiting  the  options  available  for  paint  removal.  Prior  to  1985  open 
air  blasting  with  grit  was  the  standard  method  of  paint  removal.  While  inexpensive,  it 
resulted  in  paint  debris  scattered  over  a  wide  area.  Today,  open  air  blasting  is  illegal  and 
has  been  replaced  with  removal  in  containment.  Containment  means  that  the  bridge  or  a 
portion  of  the  bridge  is  enclosed,  thereby  containing  the  paint  debris. 

Additional  regulations  have  been  enacted  limiting  Volatile  Organic  Compounds 
(VOCs)  for  architectural  and  maintenance  coatings.  These  regulations  have  resulted  in 
the  application,  into  bridges,  of  modified  paint  systems  with  limited  available  information 
associated  with  their  durability.  The  use  of  these  paints  to  replace  or  repair  existing 
systems  is  a  risk  that  the  states  must  take  in  order  to  meet  regulatory  compliance. 

The  economic  health  of  a  nation  is  dependent  on  its  ability  to  engage  in 
commerce.  This  ability  is  directly  related  to  the  capability  of  its  infrastructure  to 
efficiently  and  safely  respond  to  the  demands  placed  upon  it,  not  only  by  its  users,  but 
also  by  the  environment.  A  recent  survey  indicates  that  of  the  nearly  600,000  bridges 
tabulated,  just  over  1 90,000  bridges  were  considered  substandard.  While  the  reasons  for 
this  classification  are  varied,  a  growing  number  are  the  result  of  the  presence  of  lead- 
containing  paints  previously  applied  for  corrosion  protection.  Both  recently  adopted  and 
proposed  future  regulations  have  resulted  from  a  growing  awareness  of  the  need  to 
protect  the  environment  from  uncontrolled  pollution,  and  to  safeguard  the  health  of 
workers  engaged  in  renovation  as  well  as  that  of  the  general  populace.  Under  this 
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research  project,  INDOT  needs  to  evaluate  the  existing  and  new  paint  systems  from  cost 
point  of  view.  Therefore,  a  life  cycle  cost  analysis  was  done  for  the  existing  and  some 
new  paint  systems  to  select  the  best  paint  type  that  is  convenient  to  Indiana. 

STUDY  OBJECTIVES 

The  purpose  of  this  project  was  to  review  aspects  of  the  rehabilitation  process  in 
light  of  the  above  mentioned  limitations,  assess  to  the  degree  possible  the  state  of  the  art, 
arrive  at  conclusions  and  make  recommendations  where  applicable. 

The  objective  of  this  study  was  to  perform  an  economic  study  on  INDOT  steel 
bridge  paint  maintenance  problems  and  life  cycle  cost  analysis.  A  major  goal  of  this 
study  was  to  develop  economic  models  that  can  be  used  to  provide  a  rational  framework 
for  the  evaluation  of  alternatives  in  the  paint  maintenance  of  steel  bridges.  To  accomplish 
the  objective,  an  extensive  study  of  steel  bridge  maintenance  practices  was  conducted. 
The  purpose  of  this  effort  was  to  acquire  cost  data  and  detailed  information  on  practices 
and  performance  experience,  and  gain  a  better  imderstanding  of  the  bridge  maintenance 
problems  as  viewed  from  the  owner's  perspective.  This  study  has  included  a  literature 
search  and  a  series  of  meetings  and  discussions  with  various  groups  and  individuals 
within  the  bridge  maintenance  community,  including  various  state  highway  department 
persormel,  representatives  of  the  paint  industry,  the  different  departments  of 
transportation  (DOTs)  and  bridge  painting  contractors.  The  data  and  experience  from  the 
bridge  paint  maintenance  study  were  used  to  formulate  the  models  and  to  provide  input 
data  for  the  completed  models. 

The  sub-objectives  of  this  study  are: 

1  -  Study  the  deterioration  models  for  the  existing  paint  systems  using  the  deterministic 
and  probabilistic  methods.  The  deterministic  method  such  as  the  regression  analysis 
and  the  stochastic  method  such  as  the  Markov  chains  process. 

2-  Compare  the  different  existing  paint  systems  according  to  life  cycle  cost  analysis 
using  economic  traditional  methods  and  the  Markov  decision  process  method. 
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3-  Compare  the  best  existing  system  with  some  new  paint  systems  using  the  traditional 
economic  methods.  This  step  will  identify  the  best  solution  that  INDOT  can  use  in  the 
future. 

4-  Analyze  the  new  system  according  to  other  states'  experience  for  that  type  of 
paint.  This  analysis  includes  a  deterioration  model  and  life  cycle  cost  analysis  for 
different  rehabilitation  policies. 

5-  Establish  an  INDOT  maintenance  plan  according  to  the  analysis  that  is  made  in  the 
previous  step. 

STUDY  STEPS 

The  steps  that  are  applied  to  this  study  are: 

1-  Data  is  collected  from  INDOT  on  the  condition  rating  and  paint  age. 

2-  Data  is  classified  into  different  existing  paint  types  such  as  Lead  based  paint  and 
ZincA^inyl. 

3-  Regression  analysis  is  done  in  terms  of  paint  types  and  categories  such  as  Interstate 
roads  and  State  roads. 

4-  Markov  chains  process  is  used  to  represent  the  deterioration  models  for  data  as  a 
stochastic  method. 

5-  Markov  process  as  a  stochastic  method  is  carried  out  based  on  regression  analysis, 
where  regression  concludes  the  best  fitted  models  for  data. 

6-  Paint  type  cost  data  is  collected  from  INDOT  for  a  life  cycle  cost  analysis  of  the 
existing  systems. 

7-  Economic  traditional  methods,  such  as  present  value  (PV)  and  equivalent  uniform 
annual  cost  (EUAC),  are  used  as  a  deterministic  method  to  compare  different  paint  ■ 
rehabilitation  scenarios  and  different  paint  types. 

8-  The  Markov  decision  process  (MDP)  is  used  as  a  stochastic  method  to  compare 
different  paint  rehabilitation  scenarios  for  each  paint  type. 

9-  A  conclusion  was  drawn  to  select  the  best  existing  paint  system  according  to  the 
previous  life  cycle  cost  analysis. 
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1 0-  Data  is  collected  from  different  neighbor  states  about  some  other  new  paint  systems 
that  can  be  useful  in  Indiana.  Cost  data  is  also  collected  for  these  new  types. 

11-  Comparison  between  the  existing  and  the  new  systems  is  made  to  select  the  best  paint 
system  quantitatively  or  according  to  the  economic  analysis. 

12-  Different  rehabilitation  scenarios  for  the  best  paint  system  were  analyzed  to  draw  the 
maintenance  plan  for  INDOT. 

13-  Based  on  the  previous  economic  analysis,  a  maintenance  plan  for  the  new  paint 
system  is  recommended  for  INDOT  to  be  used  in  the  future. 

14-  Conclusions  and  recommendations  are  made. 

REPORT  OVERVIEW 

This  report  includes  four  chapters.  The  overview  of  these  chapters  is  shown  in 
Figure  1 .  This  figure  indicates  that  Chapter  One  discusses  the  derivation  of  deterioration 
models  using  the  deterministic  method,  such  as  regression  analysis,  and  the  stochastic 
method,  such  as  the  Markov  chains  process.  It  also  includes  the  comparison  between  the 
regression  and  Markov  process  results.  Chapter  Two  discusses  the  life  cycle  cost  analysis 
for  the  existing  and  new  paint  systems  using  the  economic  analysis  methods  and  Markov 
decision  process  (MDP)  method.  The  life  cycle  cost  analysis  for  the  best  paint  system  is 
discussed  in  Chapter  Three.  This  chapter  indicates  also  the  maintenance  plan  for  the  new 
system.  Chapter  Four  shows  the  cited  references  that  are  used  in  this  study.  Appendix  A 
shows  the  regression  analysis  results  for  the  existing  and  new  paint  systems.  The  Markov 
transition  probability  matrix  values  are  shown  in  Appendix  B.  Markov  decision  process 
(MDP)  calculations  and  results  for  one  paint  system  is  indicated  in  Appendix  C  as  an 
example.  Appendix  D  includes  the  economic  analysis  results  for  the  existing  and  new 
systems. 
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SUMMARY  AND  CONCLUSIONS 

The  application  of  life  cycle  cost  analysis  to  the  maintenance  coating  of  steel 
bridges  represents  a  major  departure  from  the  current  practice  of  basing  coating  decisions 
on  lowest  initial  cost.  The  economic  analysis  models  and  the  Markov  Decision  Process 
(MDP)  can  be  used  in  a  variety  of  possibilities  and  leads  to  some  interesting 
observations.  A  major  problem  becomes  apparent  upon  a  careful  examination  of  the  data 
that  feeds  the  model.  The  cost  and  performance  data  is  highly  variable  and  considerable 
doubt  exists  about  the  validity  of  comparisons  between  different  sources.  The  reasons  for 
these  problems  are  many.  The  costs  are  seldom  available  in  a  detailed  breakout  format. 
Typically,  the  cost  data  are  simply  a  lump  sum  that  may  be  transformed  to  a  unit  cost  per 
area  if  the  surface  area  is  provided.  fNDOT  still  uses  the  practice  of  expressing  coating 
costs  as  a  cost  per  ton  of  steel.  Given  the  wide  range  of  geometry  used  in  bridge 
construction,  conversion  of  cost  per  ton  to  cost  per  square  foot  is  difficult.  Another  major 
factor  in  the  variability  of  the  cost  data  is  the  presence  of  hidden  costs.  For  example,  we 
may  consider  two  similar  bridges  with  the  same  siarface  preparation  and  coating  system 
that  have  different  accessibility  limitations.  In  one  case,  access  for  painting  is 
unrestricted,  while  in  the  other  case,  painting  can  only  be  performed  during  periods  of 
minimal  traffic  disruption.  For  the  first  bridge,  scaffolding  and  containment  can  be 
erected  and  left  in  place  until  the  job  is  complete,  while  the  second  bridge  requires  the 
scaffolding  and  containment  to  be  assembled  and  disassembled  before  and  after  each 
daily  painting  period.  If  all  other  costs  remain  identical,  these  two  jobs  could  easily  differ 
by  two  or  three  times  in  cost  per  square  foot.  Another  source  of  hidden  costs  results  fi"om 
typical  contractor  practices  that  stem  from  cash-flow  problems.  Contractors  will  typically 
shift  labor  costs  to  purchased  materials  so  that  they  can  front-load  their  invoicing.  This 
practice  introduces  large  variations  in  hardware  and  material  costs.  Performance  data  is 
also  badly  clouded.  The  judgment  of  when  a  coating  system  has  failed  can  be  very 
subjective,  particularly  if  no  uniform  quantitative  standards  exist  forjudging  failure. 
Additional  consideration  needs  to  be  given  to  contracting  practices.  If  life-cycle  cost 
methods  are  to  provide  meaningful  comparisons  of  alternatives,  reliable  cost  data  must  be 
available.  One  approach  that  might  aid  the  situation  would  be  uniform 
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coating  contract  standards  that  require  detailed  cost  breakouts  of  the  elements  of  the 
project.  These  cost  elements  need  to  hold  up  to  audit  standards  so  that  hidden  costs  cease 
to  exist.  This  process  will  not  be  popular  with  painting  contractors,  but  if  we  are  going  to 
employ  life-cycle  costs  as  a  criteria  forjudging  coating  alternatives,  we  must  be  sure  that 
the  data  (cost  and  performance)  are  accurate  and  that  we  are  indeed  comparing 
alternatives  on  a  rational  basis. 

Based  on  the  previous  discussion,  the  existing  and  new  paint  systems  are  analyzed 
and  come  up  to  the  conclusion  that  the  3 -coat  system,  Inorganic/Organic  Zinc  Epoxy 
Urethane,  is  the  best  and  comparative  paint  system.  INDOT  just  started  implementing  the 
3-coat  system.  Therefore,  there  were  inadequate  data  for  life  cycle  cost  analysis. 
However,  based  on  MDOT  data  set,  life  cycle  cost  analysis  seemingly  indicated  that  the 
optimal  policy  for  3-coat  system  rehabilitation  was  doing  spot  painting  every  15  years  or 
when  paint  condition  rating  reaches  7,  regardless  of  the  age  of  the  paint  on  the  steel 
bridge.  Accordingly,  a  detailed  maintenance  plan  was  developed  based  on  the  economical 
life  cycle  cost  analysis  results.  Mainly,  INDOT  steel  bridges  could  be  painted  whenever 
the  paint  age  is  1 5  years  and  the  paint  condition  rating  reaches  7  because  this  plan  will 
result  in  the  lowest  cycle  cost  for  INDOT.  However,  this  result  was  derived  from  MDOT 
data  set.  Condition  rating  7  is  defined  in  Table  E.l,  Appendix  E,  in  MDOT  description. 


RECOMMENDATIONS  AND  FUTURE  WORK 

1  -    Multimedia  could  be  used  as  a  training  tool  for  the  inspectors  to  apply  the  inspection 
criteria  to  the  steel  bridges  paint. 

2-   It  is  recommended  to  use  Markov  Decision  Process  for  life  cycle  cost  analysis  of 
Steel  Bridge  paint  in  future.  This  method  is  based  on  the  stochastic  method,  Marko\- 
chains,  where  it  needs  more  data  to  analyze  the  life  cycle  cost.  Consequently,  a  design 
for  a  program  that  uses  these  methods  together  to  provide  the  best  rehabilitation 
scenario  is  recommended  to  INDOT. 
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3-   Maximizing  the  benefit  of  budget  allocation  for  INDOT  is  very  important  issue. 
Consequently,  INDOT  steel  bridge  paint  rehabilitation  budget  could  be  allocated 
every  year  to  different  bridges  according  to  their  importance.  Dynamic  Programming 
or  Mixed  Integer  Programming  techniques  are  recommended  to  be  used  based  on  the 
Markov  chains  results  to  accommodate  or  to  maximize  the  benefit  of  INDOT  budget 
allocation  for  steel  bridge  paint  rehabilitation  projects. 
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CHAPTER    I 

PAINT  PERFORMANCE  ANALYSIS:  DETERIORATION  MODELS 

A  performance  function  is  the  relationship  between  bridge  paint  condition  rating 
and  age  that  reflects  the  level  of  service  of  that  paint.  The  performance  functions  for  steel 
bridge  paint  in  INDOT  were  developed  using  regression  method.  Bridge  paint 
performance  prediction  models  were  also  developed  using  Markov  chains.  The 
probabilistic  model  that  was  developed  by  Markov  chains  was  done  to  reflect  the 
stochastic  nature  of  bridge  paint  conditions.  This  model  can  be  used  to  predict  the 
condition  rating  of  a  bridge  paint  at  a  given  age. 

This  study  used  the  techniques  of  regression,  Markov  chains,  non-linear 
programming  and  a  combination  of  those  techniques  to  analyze  bridge  paint  performance. 
The  results  exhibited  the  power  of  those  techniques,  especially  of  Markov  chains 
approach  in  predicting  or  estimating  future  bridge  paint  conditions. 

1.1.  DEVELOPMENT  OF  PERFORMANCE  FUNCTIONS  USING 
DETERMINISTIC  METHOD 

I.l.l.  Deterministic  Method  (Regression  Analysis) 

There  are  many  factors  that  affect  the  performance  of  steel  bridge  paint.  Some  of 
these  factors  are  climate,  age  and  traffic  conditions  (interstate  or  state  road).  The  data 
available  is  divided  by  road  type:  interstate  steel  bridges  and  state  steel  bridges. 
Regression  analysis  by  using  SAS  1996  is  done  to  check  if  climate  conditions  have  an 
effect  on  the  performance  function  or  not.  Indiana  State  is  divided  into  v^\o  different 
climate  regions:  north  (INDOT  Districts  1  and  3)  and  south  (INDOT  Districts  2,  4,  5,  and 
6).  Since  we  have  no  data  on  weather  conditions,  the  climate  effect  is  analyzed  by  putting 
it  in  as  a  class  variable  with  one  indicator  variable  in  the  SAS  analysis.  This  indicator 
variable  is  a  qualitative  variable  that  has  value  1  if  it  is  in  the  north  and  0  if  it  is  in  the 
south.  Different  types  of  functions  are  used  to  express  the  performance  functions  to  get 
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the  best  fit.  Data  is  checked  for  Hnear  and  polynomial  (quadratic  and  cubic)  functions  and 
climate  as  a  class  variable  is  also  checked  for  different  paint  types  and  road  conditions. 
The  F-test  and  t-test  are  performed  to  determine  the  best  model  and  whether  climate  is 
significant.  The  assumptions  for  the  regression  models  are  checked  to  select  the 
appropriate  model  to  fit  the  data.  SAS  results  are  indicated  in  Appendix  (A)  that  contains 
linear,  quadratic  and  cubic  models  and  their  F-test  and  t-test  results.  A  lack-of-fit  test  is 
performed  for  all  the  models  because  there  are  replications  in  the  data.  This  test  is  done  to 
check  the  adequacy  of  these  models  to  fit  the  data. 

Deterioration  models  information  is  indicated  in  Table  1. 1 .  It  contains  all  the 
models  that  are  built  using  SAS,  along  with  some  vital  information  about  them  to  make 
comparison  among  them  and  select  the  best  model  that  can  represent  the  data  or  fit  the 
data.  <'^' 

For  paint  type  (1),  interstate  roads  have  their  performance  fiinction  in  linear, 
quadratic  and  cubic  formulas.  The  r-squared  value  for  each  formula  is  calculated  and 
indicated  in  the  table,  but  the  difference  among  them  is  very  low.  The  lack-of-fit  test  is 
done  for  each  case  and  it  is  not  significant  for  any  formula.  This  means  that  all  the 
models  are  adequate  to  fit  the  data.  The  assumptions  for  the  regression  models  are 
checked  for  the  error  constant  variance,  normality  and  independence.  The  results  indicate 
that  all  the  models  are  good  according  to  these  assumptions  as  indicated  in  Appendix  A. 
The  climate  effect  is  measured  when  this  factor  is  included  in  the  regression  analysis  as  a 
class  variable  with  one  indicator  variable.  The  significance  of  that  class  variable  or 
indicator  variable  is  checked  for  all  the  models  of  that  paint  type  and  indicates  that  it  is 
not  significant  enough  to  be  included  any  model.  The  p-value  of  this  test  of  significance 
of  the  climate  factor  is  indicated  in  Table  I.l  where  all  of  the  models  indicate  negative 
answers  for  the  significance  of  that  factor.  In  addition  to  this  information,  the  significance 
of  all  models'  parameters  is  checked  to  select  the  best  model  that  can  fit  the  data. 
According  to  these  tests  and  arguments,  the  best  model  is  selected  for  paint  type  (1) 
Interstate  roads  is  a  quadratic  model:  ''^^ 

Paint  Rating  =  9,06-  0. 0821  *  Age  -  0. 001 78  *  Age^2. 
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For  Paint  type  (2),  interstate  roads  have  their  performance  function  in  linear  and 
cubic  formulas.  The  r-squared  value  for  each  formula  is  calculated  and  indicated  in  the 
table,  but  the  difference  among  them  is  very  low.  The  lack-of-fit  test  is  done  for  each  case 
and  it  is  significant  for  10%  and  5%  levels  of  significance,  but  it  is  not  for  1%.  This 
indicates  that  all  the  models  are  adequate  for  fitting  the  data  at  only  1  %  level  of 
significance.  The  assumptions  of  the  regression  models  are  checked  for  the  error  constant 
variance,  normality  and  independence.  The  results  indicate  that  all  the  models  are  good 
according  to  these  assumptions.  The  climate  effect  is  measured  when  this  factor  is 
included  in  the  regression  analysis  as  a  class  variable  with  one  indicator  variable.  The 
significance  of  that  class  variable  or  indicator  variable  is  checked  for  all  the  models  of 
that  paint  type  and  indicates  that  it  is  not  significant  enough  to  be  included  any  model. 
The  p-value  of  this  test  of  significance  of  the  climate  factor  is  indicated  in  Table  I.l 
where  all  of  the  models  indicate  negative  answers  for  the  significance  of  that  factor.  In 
addition  to  this  information,  the  significance  of  all  models'  parameters  is  checked  to 
select  the  best  model  that  can  fit  the  data.  According  to  these  tests  and  arguments,  the 
best  model  is  selected  for  paint  type  (2)  Interstate  roads  is  a  cubic  model:  ''^^ 

Paint  Rating  =  9.06-  0.201  *  Age  +  0. 0103  *  Age''2  - 
0.000348  *Age''3. 

For  Paint  type  (1).  state  roads  have  their  performance  function  in  linear,  quadratic 
and  cubic  formulas.  The  r-squared  value  for  each  formula  is  calculated  and  indicated  in 
the  table,  but  the  difference  among  them  is  low.  The  lack-of-fit  test  is  done  for  each  case 
and  it  is  not  significant  for  quadratic  and  cubic  formulas  and  it  is  only  significance  at  1% 
level  for  linear  formula.  This  indicates  that  all  the  models  are  adequate  to  fit  the  data.  The 
assumptions  of  the  regression  models  are  checked  for  the  error  constant  variance, 
normality  and  independence.  The  results  indicate  that  all  the  models  are  good  according 
to  these  assumptions.  The  climate  effect  is  measured  where  this  factor  is  included  in  the 
regression  analysis  as  a  class  variable  with  one  indicator  variable.  The  significance  of  that 
class  variable  or  indicator  variable  is  checked  for  all  the  models  of  that  paint  type  and 
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indicates  that  it  is  not  significant  to  be  included  in  any  model.  The  p-value  of  this  test  of 
significance  of  the  climate  factor  is  indicated  in  Table  I.l  where  all  of  the  models  indicate 
negative  answers  for  the  significance  of  that  factor.  In  addition  to  this  information,  the 
significance  of  all  models'  parameters  is  checked  to  select  the  best  model  that  can  fit  the 
data.  According  to  these  tests  and  arguments,  the  best  model  is  selected  for  paint  type  (\) 
State  roads  is  a  quadratic  model:  ^'^^ 

Paint  Rating  =  9.06-0.007*  Age  -  0. 0051 7  *  Age  ^2. 

For  Paint  type  (2),  state  roads  have  their  performance  function  in  linear  and  cubic 
formulas.  The  r-squared  value  for  each  formula  is  calculated  and  indicated  in  the  table, 
but  the  difference  among  them  is  very  low.  The  lack-of-fit  test  is  done  for  each  case  and 
it  is  not  significant  for  all  formulas.  This  means  that  all  the  models  are  adequate  to  fit  the 
data.  The  assumptions  of  the  regression  models  are  checked  for  the  error  constant 
variance,  normality  and  independence.  The  results  indicate  that  all  the  models  are  good 
according  to  these  assumptions.  The  climate  effect  is  measured  where  this  factor  is 
included  in  the  regression  analysis  as  a  class  variable  with  one  indicator  variable.  The 
significance  of  that  class  variable  or  indicator  variable  is  checked  for  all  the  models  of 
that  paint  type  and  indicates  that  it  is  not  significant  to  be  included  in  any  model.  The  p- 
value  of  this  test  of  significance  of  the  climate  factor  is  indicated  in  Table  I.l  where  all  of 
the  models  indicate  negative  answers  for  the  significance  of  that  factor.  In  addition  to  this 
information,  the  significance  of  all  models'  parameters  is  checked  to  select  the  best 
model  that  can  fit  the  data.  According  to  these  tests  and  arguments,  the  best  model  is 
selected  for  paint  type  (2)  State  roads  is  a  cubic  model:  ^'^* 

Paint  Rating  =  9.03 -0.0753  "^  Age  -  0.00489  -^  Age''2  + 
0.000054  *Age^3. 

For  Paint  type  (3),  interstate  and  state  roads  are  combined  in  one  model  and  ha\e 
their  performance  function  in  linear  formula  only  due  to  lack  of  data.  The  r-squared  value 
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is  calculated  and  indicated  in  the  table.  The  lack-of-fit  test  is  done  and  it  is  not  significant 
for  this  model.  This  means  that  the  model  is  adequate  to  fit  the  data.  The  assumptions  of 
the  regression  model  are  checked  for  the  error  constant  variance,  normality  and 
independence.  The  results  indicate  that  the  model  is  good  according  to  these  assumptions. 
The  climate  effect  is  measured  where  this  factor  is  included  in  the  regression  analysis  as  a 
class  variable  with  one  indicator  variable.  The  significance  of  that  class  variable  or 
indicator  variable  is  checked  for  the  model  of  that  paint  type  and  indicates  that  it  is  not 
significant  to  be  included  in  any  model.  The  p-value  of  this  test  of  significance  of  the 
climate  factor  is  indicated  in  Table  I.l  where  the  model  indicates  negative  answers  for  the 
significance  of  that  factor.  In  addition  to  this  information,  the  significance  of  model's 
parameters  is  checked  to  select  the  best  model  that  can  fit  the  data.  According  to  these 
tests  and  arguments,  the  best  model  is  selected  for  paint  type  (3)  is  a  linear  model:  ''^* 


Paint  Rating  =  8.88-0.123*  Age. 


1.2.  STOCHASTIC  METHOD  (MARKOV  CHAINS  PROCESS) 

Stochastic  processes  are  processes  that  evolve  over  time  in  a  probabilistic  manner. 
A  stochastic  process  is  defined  to  be  an  indexed  collection  of  random  variables  (St), 
where  the  index  t  runs  through  a  given  set  of  non-negative  integers.  One  special  type  of 
stochastic  process  is  called  Markov  chain.  A  stochastic  process  is  a  Markov  chain  if  it  has 
the  Markovian  property;  the  conditional  probability  of  any  future  event,  given  any  past 
event  and  the  present  state  St  =  i,  is  independent  of  the  past  event  and  depends  only  upon 
the  present  state.  This  property  can  be  written  as:  '"'^^  ^'^ 

P  (St+i  =  it+i  /  St  =  it ,  St-i  =  it-i Si  =  ii.  So  =  io)  =  P  (St+i  =  it+i  /  St  =  it). 

Many  processes  fit  this  description  including  this  study  of  steel  bridge  paint 
conditions.  In  the  model  development,  to  reduce  the  complexity  of  the  analysis,  the  future 
condition  of  bridge  paint  is  assumed  to  depend  only  on  the  present  state,  and  independent 
of  the  past  conditions.  It  is  further  assumed  that  for  all  states  i ,  j  and  all  t,  P  (St+i  =  it+i  / 
St  =  it)  is  independent  oft.  The  probability  ,  Pij .  that  bridge  paint  is  in  state  i  at  time  t 
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and  it  will  be  in  a  state  j  at  time  t+1  does  not  change  (  remains  stationary)  over  time.  This 
stationary  assumption  is  expressed  by  the  following  equation:     P  (St+i  =  j  /  St  =  i)  =  Pij 

The  term  transition  is  used  when  the  system  moves  from  state  i  during  one  period 
to  state  j  during  the  next  period.  Accordingly,  the  probabilities,  Fij's  ,  are  referred  to  the 
transition  probabilities.  The  transition  probabilities  are  commonly  displayed  as  n  x  n 
matrix  called  the  transition  probability  matrix  P.  In  this  study  there  are  five  states 
associated  with  the  five  possible  conditions  of  bridge  paint  ratings.  State  1  corresponds  to 
the  best  condition  and  state  5  corresponds  to  the  worst  condition.  Then,  the  transition 
probability  matrix  can  be  written  as: 


Pll  Pl2  Pl3  Pl4  Pl5 

P21  P22  P23  P24  P25 

P31  P32  P33  P34  P35 

P41  P42  P43  P44  P45 

P51  P52  P53  P54  P55 


The  above  transition  probability  matrix  is  for  one-step  (one-period)  transition. 
The  n-step  transition  probability  matrix,  P^"',  of  the  process  that  is  in  state  i  and  will  in 
state  j  after  n  periods  is  computed  by  Chapman-Kolmogorov  equation:  P^"^=  P"  .  The  n- 
step  transition  probability  matrix  is  obtained  by  taking  the  n-th  power  of  the  one  step 


transition  matrix, 


(22)  (1) 


The  Markov  chain  as  applied  to  bridge  paint  performance  prediction  is  based  on 
the  concept  of  defining  states  in  terms  of  bridge  paint  condition  ratings  and  obtaining  the 
probabilities  of  paint  condition  changing  from  one  state  to  another.  These  probabilities 
are  represented  in  a  matrix  form  that  is  called  the  transition  probability  matrix.  Knov-ing 
the  present  state  of  bridge  paint,  or  the  initial  state,  the  future  conditions  can  be  predicted 
through  the  multiplication  of  initial  state  vector  and  the  transition  probability  matrix. 

(13),(12) 
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1.2.1.  Development  of  Markov  Prediction  Model 

Steel  bridge  paint  begins  its  life  in  a  near-perfect  condition  and  is  then  subjected 
to  a  sequence  of  duty  cycles  that  cause  the  paint  condition  to  deteriorate.  In  this  study  the 
paint  condition  rating  is  defined  by  an  inspector's  judgement.  The  paint  rating  ranges 
from  0  to  9,  where  9  is  the  perfect  condition  and  0  is  the  worst  condition.  A  duty  cycle  for 
a  paint  type  is  defined  as  one  year's  duration  of  weather  and  traffic.  A  state  vector 
indicates  the  probability  of  bridge  paint  being  in  each  of  the  10  conditions,  0-9,  in  any 
given  year.  Figure  I.l  shows  a  stochastic  representation  of  state,  state  vector  and  duty 
cycle  for  pavement  conditions.  ^'**''^'  The  description  of  each  of  these  conditions  is 
indicated  in  Table  E.l,  Appendix  E. 

It  is  assumed  that  all  bridge  paint  conditions  are  in  state  1  (condition  9)  at  an  age 
of  0  year.  Thus  the  state  vector  in  Duty  Cycle  0  (age=0)  is  given  by  (1,0,0,0,0)  because  it 
is  known  that  all  the  bridges'  paint  must  be  in  state  1  at  an  age  of  0  year  with  a 
probability  of  1 .0.  ^'''  '^'  To  model  the  way  in  which  bridge  paint  deteriorates  with  time,  it 
is  helpful  to  establish  a  Markov  probability  transition  matrix.  In  this  research,  the 
assumption  is  made  that  the  bridge  paint  condition  will  not  drop  by  more  than  one  state  in 
a  single  year.  Thus,  the  bridge  paint  will  either  stay  in  its  current  state  or  degrade  to  the 
next  lower  state  in  1  year.  It  should  be  desirable  to  examine  historical  data  for  each 
bridge,  and  determine  how  often  a  paint  condition  dropped  by  more  than  two  units 
between  biennial  bridge  inspections.  However,  only  recent  data  on  bridge  paint  condition 
have  been  preserved.  For  an  alternative  analytical  method.  Table  1.2  is  used  to  support  the 
assumption  that  the  bridge  paint  condition  may  not  drop  by  more  than  one  state  in  a 
single  year.  Based  on  fNDOT  data,  the  deterioration  rate  index  (year(s)  per  unit  change  in 
condition  rating)  is  calculated.  This  index  is  found  by  dividing  the  paint  age  by  the 
difference  between  the  highest  condition  rating,  which  is  9,  and  the  current  paint 
condition  rating.  After  calculating  this  index  for  each  bridge,  an  average  of  the  indices  for 
bridges  that  have  the  same  current  condition  rating  is  calculated.  Table  1.2  shows  these 
averages,  and  their  corresponding  standard  deviations.  For  the  189  bridges  in  the  ESIDOT 
data  set  that  zinc -based  paint  at  rating  6,  for  example,  the  mean  index  is  4.254  years  per 
unit  change  in  paint  rating.  Our  assumption  is  that  the  index  is  never  (or  seldom)  below 
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Table  1.2:  Markov  chains  Assumption  Validation 


Condition  Rating 

Paint  type                                                        | 

Lead  Based 

Zinc  Based                j 

Ave.  Rate* 

STDEV 

Under 1yr 

Ave.  Rate* 

STDEV 

Under  1yr 

5 

5.43 

1.065 

None 

4.208 

0.681 

None 

6 

7.644 

0.676 

None 

4.254 

1.719 

None 

7 

10.66 

1.962 

None 

4.031 

2.294 

3.90% 

8 

NA 

NA 

NA 

4.443 

3.149 

None 

NA  =  Not  Available. 

'  These  rates  are  calculated  by  dividing  the  bridge  age  by  the 
difference  between  the  highest  condition  rating,  which  is  9,  and 
the  current  condition  rating  for  that  bridge,  as  follows: 
Deterioration  Rate  Index  =  Age/(9  -  current  condition  rating) 
After  calculating  the  rates,  an  average  is  calculated  for  these 
rates  for  bridges  having  the  same  condition  rating.  All  these  values 
are  based  on  INDOT  data. 
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one.  Looking  at  the  individual  indices  that  went  into  the  mean  value,  we  find  that  only 
3.9%  of  the  index  values  of  condition  rating  7  in  zinc-based  paints  were  less  than  1.0. 
This  justifies  our  assumption,  which  is  invoked  to  reduce  the  computational  effort 
associated  with  the  Markov  method.  Consequently,  the  probability  transition  matrix  has 
the  following  form: 


P  = 


p(l)  q(l)      0       0       0 

0  p(2)  q(2)    0       0 

0  0  p(3)    q(3)   0 

0  0  0        p(4)  q(4) 

0  0  0        0         1 


where  p(j)  is  the  probability  of  a  bridge  paint  staying  in  state  j  during  one  duty 
cycle,  and  q())  =  1  -  p(j)  is  the  probability  of  the  bridge  paint  transiting  down  to 
next  state  (j+l)  during  one  duty  cycle.  The  entry  of  1  in  the  last  row  of  the 
transition  matrix  corresponding  to  state  5  indicates  an  "absorbing"  state.  The 
bridge  paint  can  not  transit  from  this  state  unless  repair  action  is  performed.  The 
state  vector  for  any  duty  cycle  n  is  obtained  by  multiplying  the  initial  state  vector 
S(0)  by  the  transition  matrix  P  raised  to  the  power  of  n.  Then,  ^'^'  ''^^ 

S(l)=  S(0)*P. 

S(2)=  S(l)*P  =  S(0)*Pl 


S(n)=  S(n-l)*P  =  S(0)*P 


*  pn 


With  this  procedure,  if  the  transition  probability  matrix  can  be  obtained,  the  future 
state  of  the  bridge  paint  can  be  predicted  at  any  duty  cycle,  n.  It  should  be  noticed  that 
the  lowest  rating  number  before  a  bridge  paint  is  repainted  again  is  5  (as  indicated  by 
FHWA  and  DSfDOT).  Consequently,  the  corresponding  transition  probability  p(5)  is  equal 
tol. 
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This  Markovian  model  provides  a  reliable  mechanism  for  developing  prediction 
models.  The  markov  process  imposes  a  rational  structure  on  the  deterioration  model 
because  it  explains  the  deterioration  as  uncertain  issue  and  it  also  ensures  that  the 
projections  beyond  the  limits  of  data  will  continue  to  have  worsening  condition  pattern 
with  age.  This  model  has  been  successfully  used  in  other  types  of  infrastructure 
deterioration  modeling  such  as  pavement  and  bridges.  ^'''  *''  *^'  *^^'  ^'' 

1.2.2.  Transition  Probability  Matrix  Determination 

To  estimate  the  transition  probability  matrix,  a  non-linear  programming  approach 
is  used.  The  objective  is  to  determine  the  five  parameters,  p(l)  through  p(5),  that  would 
minimize  the  absolute  distance  between  the  actual  data  points  and  the  expected 
(predicted)  bridge  paint  condition  for  the  corresponding  age  generated  by  the  Markov 
chain  using  these  five  parameters.  Due  to  the  large  number  of  data  and  the  complexity  of 
applying  this  method  for  each  data  available,  the  regression  models  output  values  are 
used  instead  of  this  data  to  construct  the  transition  probability  matrix.  This  is  because  the 
regression  models  values  are  the  output  of  least  square  method  which  is  the  best  way  to 
predict  the  data.  Based  on  the  fact  that  regression  model  is  the  best  fit  of  data  using  the 
least  square  method,  we  are  confident  that  these  models  are  the  best  we  can  use  in 
predicting  the  paint  conditions.  Consequently,  if  markov  probability  matrices  are 
constructed  according  to  these  models,  it  will  give  the  best  fit  for  the  data  using  the 
stochastic  approach.  Then,  the  objective  function  of  the  non-linear  programming  is  to 
minimize  the  absolute  difference  between  the  average  condition  rating  ate  time  t 
estimated  by  the  regression  function  and  the  estimated  value  of  condition  rating  using 
Markov  chain  at  time  t.  ^'*'  *^' 

The  objective  fiinction  has  the  following  form:  ^^'  ''^  "^^  ^'^'  ^'^  ^^^^ 

A' 

Minimize    Z  \Y(t)  - E(t,P)\ 
t=i 

Subject  to:  0<p(i)<1.0,  i  =  1, 2, ,1. 
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Where: 

N  =  4,  the  number  of  years  in  one  age  range. 
1  =  4,  the  number  of  unknown  probabilities. 

P  =  [p(l),  p(2), ,p(I)],  a  vector  of  length  I. 

Y(t)  =  the  average  condition  ratings  at  time  t,  estimated  by  regression  function. 
E(t,P)  =  the  estimated  value  of  condition  rating  by  Markov  chain  at  time  t. 

The  solution  to  this  function  is  obtained  by  using  GAMS  program  for  solving 
non-linear  programming.  This  program  is  available  on  the  Engineering  Computer 
Network  system  (ECN)  at  Purdue  University.  The  values  of  the  corresponding  regression 
function  are  taken  as  the  average  condition  rating  to  solve  the  non-linear  programming. 

The  bridge  paint  life  is  divided  into  age  zones  or  ranges  of  4  years.  These  4-year 
ranges  are  selected  to  facilitate  the  calculation  of  probabilities  and  optimization 
procedure.  It  is  assumed  that  each  zone  or  range  has  a  constant  rate  of  deterioration  and, 
hence,  that  a  constant  duty  cycle  has  been  assumed  within  each  range.  The  rate  of 
deterioration  is  assumed  to  vary  from  range  or  zone  to  another;  therefore,  different  duty 
cycles  have  been  assigned  to  different  ranges  or  zones.  '^'-^*'^' 

Because  the  duty  cycle  within  a  range  is  assumed  to  be  constant,  a  homogeneous 
Markov  chain  has  been  used  for  each  range  and  a  separate  transition  matrix  has  been 
developed  for  each  range.  The  duty  cycle  varies  from  range  to  another.  Therefore,  a  non- 
homogeneous  Markov  chain  has  been  used  for  transition  from  range  to  another. ''"'  ^'■'^ 

The  maximum  rating  of  bridge  paint  condition  is  9  and  it  represents  a  near-perfect 
condition  of  bridge  paint.  It  is  almost  always  true  that  a  new  bridge  paint  has  condition 
rating  of  9.  In  other  words,  abridge  paint  at  age  0  year  has  condition  rating  of  9  with  1 .0 
probability.  Therefore,  the  initial  state  vector  S(0)  for  a  new  bridge  paint  is  always 
[1,0,0,0,0],  where  the  numbers  are  the  probabilities  of  having  condition  rating  of  9,8,7,6 
and  5  at  age  0  year  respectively.  That  is  the  initial  vector  of  the  first  range  for  developing 
the  bridge  performance  curve  or  transition  probability  matrix.  Second  range  takes  the  last 
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State  vector  of  first  range  as  its  starting  state  vector.  Similarly,  range  or  zone  n  takes  the 
last  state  vector  of  range  n-I  as  its  starting  state  vector.  The  rest  of  the  work  to  obtain  the 
overall  bridge  paint  performance  curve  is  nothing  rather  than  conducting  the  following 
matrix  multiplication:  ^'-^^'^' 

S(l)=  S(0)*P. 

S(2)=  S(1)*P  =  S(0)*P-. 

S(3)  =  S(2)  *  P  =  S(0)  *  P^ 


in-l 


S(n-1)=  S(n-2)  *  P  =  S(0)  *  P" 
S(n)=  S(n-l)*P  =  S(0)*P". 

Where:  S(n)  represents  the  condition  state  vector  at  age  n. 
Let  R  be  the  column  vector  of  condition  ratings. 


9 
8 
R=       7 

6 

5 


Then,  the  estimated  condition  rating  at  age  n  by  Markov  chain  is, 

E(n,P)  =  S(n)  "^  R 
For  illustration,  let  the  performance  function,  using  regression,  for  the  steel  bridge  paint 
gives  the  values  of  Y(l),  Y(2),  Y(3)  and  Y(4)  for  the  predicted  condition  ratings  in  the 
first  four  years  of  paint  life.  The  corresponding  values  of  prediction  of  condition  ratings 
by  Markov  chain  method  can  be  expressed  by  the  following  equations:  ''^*^'-^' 

E(1,P)=S(1)  *R=S(0)  *P*R 

E(2,P)  =  S(l)  *R=  S(0)  *p^  *R 

E(3,P)  =  S(l)  *R=S(0)  *P^  ''R 

E(4,P)  =  S(l)  *R  =  S(0)*P'  *R 
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Since  the  E(n,P)s  are  functions  of  p(l),  p(2),  p(3)  and  p(4),  the  non-Hnear 
programming  objective  function  should  be  solved  to  get  the  values  of  these  probabilities. 
Based  on  these  probabilities'  values,  the  Markov  transition  probability  matrices  could  be 
constructed  using  these  values.  *'^^<'^' 

1.2.3  Steps  of  NLP  Application 

1-  Paint  age  is  divided  into  4-year  ranges  with  four  unknown  probabilities  in  each 
range. 

2-  Maple  mathematics  program  is  used  to  prepare  the  file  of  matrix  multiplication  as 
illustrated  in  the  previous  section.  The  output  of  this  file  is  the  equations  that  contain 
p(l),  p(2),  p(3)  and  p(4)  as  unknowns.  Therefore,  there  are  four  equations  with  four 
unknowns  that  can  be  solved  to  get  the  values  of  p(l),  p(2),  p(3)  and  p(4). 

3-  The  output  of  the  previous  step  is  substituted  to  the  non-linear  programming 
equation  as  E(t,P),  which  is  subtracted  fi-om  the  part  Y(t).  This  part,  Y(t),  is 
calculated  using  the  regression  formulas  for  the  various  models'  types. 

4-  After  constructing  these  equations,  they  are  put  into  a  GAMS  program  input  file  to 
solve  the  NLP  problem  based  on  the  constraints  that  are  illustrated  in  the  previous 
section.  This  input  file  requests  from  GAMS  to  make  minimization  for  the  absolute 
value  of  Y(t)-E(t,P)  for  each  four  years  of  paint  age.  The  output  of  this  step  are  the 
values  p(l),  p(2),  p(3)  and  p(4)  that  minimize  the  NLP  objective  function. 

5-  By  knowing  the  values  of  these  probabilities  in  each  4-year  range,  the  last  year  state 
vector  for  each  range  can  be  calculated  by  using  the  Maple  mathematics  program. 
This  state  vector  is  calculated  to  be  used  as  the  initial  state  vector  for  the  ne.xt  range 
of  four  years. 

6-  Steps  from  1  to  5  are  repeated  for  each  4-  year  range  of  paint  age. 

7-  Steps  from  1  to  6  are  repeated  also  for  various  paint  categories  such  as:  Interstate 
roads  paint  type  (1)  and  (2),  State  roads  paint  type  (1)  and  (2)  and  paint  type  (3). 

8-  Figures  from  1.2  to  1.4  shows  the  input  file  used  for  Maple  to  get  probability- 
equations,  the  input  file  for  GAMS,  and  the  input  file  for  Maple  to  get  the  last  state 
vector.  ^•''^'''('''^''^<'^' 
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Figure  1.2:  Maple  Input  File  to  get  Probabilities  Equations. 

maple 
with(linalg): 

P:=matrix(5,5,[pl,l-pl,0,0,0,0,p2,l-p2,0,0,0,0,p3,l-p3,0,0,0,0,p4,l-p4,0,0,OAl]); 

S(0)  :=matrix(l,5,[.8422907597  ,  .03693100731  ,  .03855028105  ,  .04024055068  , 

.0419874012]); 

R:=matrix(5  J, [9,8,7,6.5]); 

P2  :=multiply(P,P); 

P3  :=muItiply(P2,P); 

P4  :=multiply(P3,P); 

E(1,P)  :=multiply(S(0),P,R); 

E(2,P)  :=multiply(S(0),P2,R); 

E(3,P)  :=multiply(S(0),P3,R); 

E(4,P)  :=multiply(S(0),P4,R); 

Output: 

E(1,P):= 

.842290759  *pl  +  7.639284574  +  .0369310073  *p2  +  .038550281 1  *p3 
+  .0402405507*  p4; 

E(2,  P)  := 

7.580616837  *pl**2  +  6.738326078  *pl  *(1  -  pi)  +  6.738326078  *(1  -  pi)*  p2 
+  .2954480585*  p2**2  +  5.896035318  *(1  -pl)*(l  -  p2) 
+  .2585170512  *p2  (1  -  p2)  +  .2585170512  (1  -  p2)  p3  +  .2698519674*  p3**2 
+  .2215860439  *(1  -p2)*(l  -  p3)  +  .2313016863  *p3  *(1  -  p3) 
+  .2313016863  *(1  -  p3)  *p4  +  .2414433041*  p4**2  +  .4111397594 
+  .1927514053  *(1  -p3)*(l  -  p4)  +  .2012027534  *p4  *(1  -  p4) 
- .2012027534  *p4; 

E(3,  P)  := 

7.580616837  *pl**3  +  6.738326078  *(pl**2)  *  (1  -  pi) 
+  6.738326078  *(pl*  (1  -pl)  +  (l  -  pi)  *p2)  *p2  +  .2954480585  *p2**3 
+  5.896035318  *(pl  *(1  -pl)  +  (l  -pl)*p2)*(l  -  p2) 
+  5.896035318  *(1  -*  pi)  *(1  -  p2)*  p3  +  .2585170512*  (p2**2)  *(1  -  p2) 
+  .2585170512  *(p2  (1  -  p2)  +  (1  -  p2)*  p3)*  p3  +  .2698519674  *p3**3 
+  5.053744558  *(1  -  pi)  *(1  -  p2)*  (1  -  p3) 
+  .2215860439  *(p2  *(1  -  p2)  +  (1  -  p2)*  p3)  *(1  -  p3) 
+  .2215860439  *(1  -  p2)  *(1  -  p3)  *p4  +  .2313016863  *(p3**2)  *(1  -  p3) 
+  .2313016863  *(p3  *(1  -  p3)  +  (1  -  p3)  *p4)*  p4  +  .2414433041*  p4**3 
+  . 411 1397594 +  . 1846550366*  (1  -p2)*(l  -  p3)*  (1  -  p4) 
+  .1927514053  *(p3  *(1  -  p3)  +  (1  -  p3)  *p4)*  (1  -  p4) 
''  +  .1927514053  *(1  -  p3)  *(1  -  p4)  +  .2012027534*  (p4**2)  *  (1  -  p4) 
+  .2012027534  *p4  *(1  -  p4)  -  .2012027534  *p4; 
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E(4,  P)  := 

.1927514053  *(1  -  p3)  *(1  -  p4)  -  .2012027534  *p4  +  .2012027534  *p4*  CI  -  p4) 

+  .1846550366*  (1  -  p2)*  (1  -  p3)*  (1  -  p4) 

+  .1927514053  *(p3*  (1  -  p3)  +  (1  -  p3)  *p4)  *(1  "  p4j 

+  .2012027534  *(p4**2)  *  (1  -  p4)  +  6.738326078*  (pl**3j  *  (1  -  pi  j 

+  6.738326078  *((pl**2)  *  (1  -  pi)  +  (pi*  (1  -  pi)  +  (1  -  pi)  *p2)  *p2)*  p2 

+  7.580616837  *(pl**4)  +  .2012027534* (p4**3)  *  (1  -  p4)  +  5.896035318* 

((pi  *(1  -pl)  +  (l  -pl)*p2)*(l-p2)  +  (l  -pl)*(l  -p2)*p3)*p3 

+  5.896035318  *((pl**2)*(l  -pl)  +  (pl*(l  -pl)  +  (l  -  pi)  *p2j  *p2)*  (1  -  p2j 

+  .2414433041  *p4**4  +5.053744558* 

((pi  *(1  -  pi)  +  (1  -  pi)*  p2)  *(1  -  p2)  +  (1  -  pi)  *(1  -  p2)  *p3)*  (1  -  p3j 

+  5.053744558  *(1  -  pi)  *(1  -  p2)  *(1  -  p3)*  p4  +  .2954480585  *p2**4 

+  4.21 1453799  *(1  -  pi)  *(1  -  p2)  *(1  -  p3)*  (1  -  p4) 

+  .2585170512  **(p2**3)  *  (1  -  p2) 

+  .2585170512  *((p2**2)*(l  -p2)  +  (p2*(l  -p2)  +  (l  -  p2)*  p3)  *p3)*  p3 

+  .2215860439  *((p2**2)  *(1  -  p2)  +  (p2  *(1  -  p2)  +  (1  -  p2)*  p3)  *p3)  *(1  -  p3) 

.2215860439* 

((p2  *(1  -  p2)  +  (1  -  p2)  *p3)*  (1  -  p3)  +  (1  -  p2)  *(1  -  p3)*  p4)*  p4  ^ 

.1846550366* 

((p2  *(1  -  p2)  +  (1  -  p2)  *p3)  *(1  -  p3)  +  (1  -  p2)  *(1  -  p3)  *p4)*  (1  -  p4) 

+  .2698519674  *p3**4  +  .2313016863*  (p3**3)  *  (1  -  p3) 

+  .2313016863  *((p3**2)  *(1  -p3)  +  (p3*(l  -p3)  +  (l  -  p3)  *p4)  *p4)*  p4 

+  .1927514053  *((p3**2)  *(1  -  p3)  +  (p3*  (1  -  p3)  +  (1  -  p3)*  p4)  *p4)  *(1  -  p4) 

+  .4111397594; 
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Figure  1.3:  GAMS  Input  File. 


VARIABLE  pi,  p2,  p3,  p4,  DIFl,  DIF2,  DIF3,  DIF4,  OBJ,  Ml,  M2,  M3,  M4; 

EQUATION  El,  E2,  E3,  E4,  E5,  E6,  E7,  E8,  E9,  ElO,  Ell,  E12,  E13,  E14,  E15,  E16, 

E17; 

El  ..Ml  =E=Pl+8; 

E2..  M2=E=9*pl**2  +  8*pl*(l  -pl)  +  8*(l  -pl)*p2  +  7*(l  -pl)*(l  -  p2); 

E3..  M3=E=9*pl**3  +  8*(pl**2)*(l  -  pi)  +  8*(pl*(l  -pl)  +  (l  -pl)*p2)*p2  + 

7*(pl*(l  -  pi)  +  (1  -  pl)*p2)*(l  -  p2)  +  7*(1  -  pl)*(l  -  p2)*p3  +  6*(1  -  pl)*(l  -  p2)*(l  - 

p3); 

E4  ..  M4  =E=  9*pl**4  +  8*(pl**3)*(l  -  pi)  +  8*((pl**2)*(l  -  pi)  +  (pl*(l  -  pi)  +  (1  - 

pl)*p2)*p2)*p2  +  7*((pl**2)*(l  -pl)  +  (pl*(l  -pl)  +  (l  -pl)*p2)*p2)*(l  -  p2)  + 

7*((pl*(l  -pl)  +  (l  -pl)*p2)*(l  -p2)  +  (l  -pl)*(l  -p2)*p3)*p3  + 

6*((pl*(l  -pl)  +  (l  -pl)*p2)*(l  -p2)  +  (l  -pl)*(l  -p2)*p3)*(l  -p3)  +  6*(l  -pl)*(l  - 

p2)*(l  -  p3)*p4  +  5*(1  -  pl)*(l  -  p2)*(l  -  p3)*(l  -  p4); 

E5  ..  DIFl  =E=  9.00 -Ml; 

E6  ..  DIF2  =E=  8.896  -  M2; 

E7..DIF3=E=8.74-M3; 

E8  ..  DIF4  =E=  8.60  -  M4; 

E9  ..  OBJ  =E=  ABS(DIFl)  +  ABS(DIF2)  +  ABS(DIF3)  +  ABS(DIF4); 

ElO..  P1=L=1.0; 

Ell  ..PI  =G=0.0; 

P1.L0  =  0.0001; 

E12..  P2=L=  1.0; 

E13  ..  P2=G=0.0; 

P2.L0  =  0.0001; 

E14..  P3=L=1.0; 

E15..  P3=G=0.0; 

P3.L0  =  0.0001; 

E16..  P4=L=1.0; 

E17..  P4=G=0.0; 

P4.L0  =  0.0001; 

MODEL  TAREK  /  ALL  /; 

option  dnlp  =  minosS; 

SOLVE  TAREK  USING  DNLP  MINIMIZING  OBJ; 

display  pl.l,  p2.1,  p3.1,  p4.1; 
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Figure  1.4:  Input  File  For  Maple  to  get  Last  State  Vector 


maple 

with(linalg): 

P:=matrix(5,5,[0.958,0.042,0,0,0,0,0.0001,0.9999,0,0,0,0,0.0001,0.9999AO,0,0,0.0001,0. 

9999,0,0,0,0,1]); 

S(0)  :=matrix(l,5,[  1,0,0,0,0]); 

P2 

P3 

P4 


=multiply(P,P); 

=multiply(P2,P); 

=multiply(P3,P); 
S(4)  :=multiply(S(0),P4); 
S(4)  :=[.8422907597  ,  .03693100731  ,  .03855028105  .  .04024055068  ,  .04198740126] 
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1.2.4  Application  of  the  Markov  Ctiains  Model 

Once  the  transition  probability  matrix  is  obtained,  the  prediction  of  the  fliture 
condition  by  Markov  chains  becomes  a  matter  of  matrix  multiplication.  Let  us  take  the 
paint  type  (2)  state  road  steel  bridges  performance  curve.  As  mentioned  earlier,  the  initial 
state  vector  of  the  first  range  of  four  years  for  new  paint  is  always  [1,0,0,0,0].  Therefore, 
the  major  problem  is  to  obtain  the  transition  probability  matrix  for  paint  type  (2)  State 
road."^' 

The  values  of  Y(t)  obtained  from  the  regression  function  are  used  to  solve  the 
NLP  equation.  This  solution  provides  the  transition  probability  matrix  for  different  paint 
age  ranges  of  four  years.  The  following  matrix  is  the  transition  probability  matrix  for 
state  road  paint  type  (2)  bridges  in  the  first  range  of  four  years. 


[.93 

5 

.065 

0 

0 

0      ] 

[    0 

.766 

.234 

0 

0       ] 

=    [    0 

0 

.301 

.699 

0       ] 

[    0 

0 

0 

.0001 

.9999] 

[    0 

0 

0 

0 

.1       ] 

p(l)  =  0.935  for  the  first  range  that  indicates  that  the  probability  of  state  road  bridges  of 
paint  type  (2)  of  paint  age  4  years  or  less  transition  from  state  1  (  condition  rating  9  )  to 
state  1  (  remaining  in  state  1  )  in  one  year  period  is  0.935.  and  the  probability  of  transition 
from  state  1  to  state  2  (  condition  rating  8)  is  q(l)  =  0.065.  Similarly,  p(2)  ==  0.766  for  the 
first  range  indicates  that  the  probability  of  state  road  bridges  of  paint  type  (2)  of  paint  age 
4  years  or  less  transition  from  state  2  (  condition  rating  8  )  to  state  2  (remaining  in  state  2) 
in  one  year  period  is  0.766,  and  the  probability  of  transition  from  state  2  to  state  3 
(condition  rating  7)  is  q(2)  =  0.234.  ^'^^ 

To  calculate  the  condition  rating  value  by  using  a  Markov  chains,  the  value  of 
E(t,P)  can  be  calculated  by  using  the  equations  described  before. 
S(0)  =  [1,0,0,0,0] 
E(0,P)  =  S(0)*R  =  9.0 
where:  R  is  the  vector  of  condition  rating  9,8,7,6  and  5. 
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E(1,P)  :=multiply(S(0),P,R);  E(l,  P)  =  [8.935] 

E(2,P)  :=multiply(S(0),P2,R);  E(2,  P)  =  [8.859015] 

E(3,P)  :=multiply(S(0),P3,R);  E(3,  P)  =  [8.765686375] 

E(4,P)  :=multiply(S(0).P4.R);  E(4,  P)  =  [8.647524705]  =^^'  -  '^s) 

Similarly,  the  condition  rating  for  each  year  of  paint  age  can  be  calculated  for  the 
same  type  of  paint  and  road.  Consequently,  the  deterioration  curve  for  that  type  of  paint 
can  be  drawn  using  the  values  of  condition  rating  calculated  by  using  the  Markov  chains. 
Figure  1.5  indicates  this  curve.  Table  1.3  indicates  the  paint  rating  prediction  for  different 
ages.  Table  1.4  indicates  the  predicted  ages  that  corresponding  to  different  states. 

These  calculations  could  be  made  for  all  the  paint  categories:  Interstate  Paint 
types  (1)  &  (2),  State  paint  types  (1)  &  (2)  and  Paint  type  (3).  Figures  from  1.5  to  1.9 
indicate  these  curves  respectively. 

1.3.  DIFFERENCE  BETWEEN  REGRESSION  AND  MARKOV  PROCESS 

Even  though  bridge  paint  performance  curves  have  been  developed  by  using 
regression,  it  is  still  necessary  to  use  the  Markov  chains  model  to  predict  individual 
bridge  paint  conditions.  As  a  matter  of  fact,  both  regression  and  Markov  chains  models 
play  important  roles  in  analyzing  bridge  paint  systems.  The  regression  model  can 
be  used  to  estimate  the  extent  of  condition  improvement  as  a  measure  of  effectiveness  in 
response  to  alternate  rehabilitation  and  repair  strategies.  However,  when  condition 
prediction  is  concerned,  the  Markov  chains  model  provides  more  reasonable  estimates  of 
bridge  paint  conditions.  This  procedure  is  described  in  the  following  section.  ^'^^  ^'^^ 

The  prediction  models  currently  in  use  vary  in  complexit}'  from  simple  straight- 
line  extrapolation,  regression  models,  to  probability-based.  Markov,  models.  Straight-line 
extrapolation  is  used  to  predict  the  condition  of  bridge  paint.  When  sufficient  data  is 
available,  it  is  found  that  the  shape  of  the  deterioration  curves  is  generally  curvilinear, 
rather  than  the  straight  line  that  results  from  straight-line  extrapolation. 
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The  probability-based  Markov  model  was  first  developed  for  the  Arizona  PMS  to 
describe  pavement  condition  changes.  Intuitively,  the  behavior  of  pavement  is  not 
deterministic  but  probabilistic.  Similarly,  the  behavior  of  the  steel  bridge  paint  is  also  not 
deterministic  but  probabilistic.  Consequently,  the  selection  of  an  appropriate  repair 
strategy  for  pavement  or  paint  is  also  an  uncertain  procedure.  Because  of  the  probabilistic 
nature  of  steel  bridge  paint,  it  is  decided  to  develop  probability-based  prediction  models 
to  predict  its  behavior.  ^""'*  *■*' 

The  usage  of  a  Markov  probability  decision  process  has  the  following  advantages: 

1-  Future  decisions  on  preservation  actions  are  not  fixed,  but  depend  on  how 
bridge  paint  actually  perform. 

2-  Actions  to  be  taken  now  can  be  identified.  Also,  actions  likely  to  be  taken  in 
the  next  few  years  can  be  identified  with  a  high  degree  of  probability. 

3-  It  is  possible  to  compare  the  expected  proportions  in  given  condition  states 
with  the  actual  proportions  observed  in  the  field.  In  this  way  possible  defects 
in  construction,  materials,  quality  and  so  on  can  be  identified  for  bridge  paint. 

4-  A  dynamic  decision  model  has  the  potential  for  significant  cost  savings  by 
selecting  less  conservative  rehabilitation  actions  that  will  still  satisfy  the 
prescribed  performance  standards.  ^^^  ^'"^^  ''^^ 

Regression  extrapolation  techniques  are  deterministic  and  do  not  attempt  to 
explain  the  variability  among  the  data  points;  they  merely  fit  the  best  line  to  the  data. 
Regression  techniques  are  powerful  tools,  but  in  many  cases  the  models  are  chosen  for 
the  best  fit  without  regard  to  the  suitability  or  intrinsic  relevance  to  the  variables  selected. 
Polynomials  of  different  degrees  and  mathematical  functions  can  be  manipulated  to  fit 
the  data;  but  when  these  functions  are  projected  beyond  the  bounds  of  the  data  results 
could  be  totally  misleading.  ^^^ 

It  is  known  that  the  rate  of  paint  deterioration  is  uncertain.  Therefore,  the 
predictive  model  should  reflect  this,  rather  than  using  the  erroneous  assumption  of 
deterministic  behavior.  The  Markov  process  imposes  a  rational  structure  on  the 
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deterioration  model.  This  form  of  predictive  methodology  has  the  further  advantages  of 
ensuring  that  projections  beyond  the  limits  of  the  data  will  continue  to  have  the  classic 
pattern  of  worsening  condition  with  age,  something  that  the  regression  models  cannot 
guarantee.  *^^ 

Another  advantage  of  probability  based  models  is  the  ease  with  which  they  can  be 
integrated  into  optimization  processes.  The  Markov  process  is  a  natural  tool  that  is  used 
in  alliance  with  dynamic  programming  to  produce  the  optimal  solutions.  It  is  believed 
that  the  application  of  Markov  process  in  conjunction  with  dynamic  programming  will 
produce  optimal  maintenance  and  rehabilitation  (M&R)  strategies  for  selected  steel 
bridges  quickly  and  efficiently.  ^^^ 

Since  a  performance  curve  for  steel  bridges  paint  represents  the  average  or  mean 
condition  rating  at  any  given  bridge  paint  age,  it  is  obvious  that  both  Markov  chains 
method  and  regression  method  can  be  used  to  predict  the  average  condition  ratings  for 
bridge  paint.  However,  the  Markov  process  has  great  advantages  over  regression  in 
predicting  individual  bridge  paint  condition.  ^'^* 


Table  1.3:  Paint  Rating  Prediction  Using  Regression  and  Markov. 
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Age 

Paint  Rating                                                                 | 

Interstatel 

Statel 

Intersi 

tate2 

State2 

Paints 

Reg 

Markov 

Reg 

Markov 

Reg 

Markov 

Reg 

Markov 

Reg 

Markov 

0 

9.00 

9.00 

9.00 

9.00 

9.00 

9.00 

9.00 

9.00 

9.00 

9.00 

1 

8.98 

8.94 

9.05 

9.00 

9.07 

8.94 

8.95 

8.94 

8.76 

8.79 

2 

8.89 

8.88 

9.03 

8.99 

8.90 

8.88 

8.86 

8.86 

8.63 

8.63 

3 

8.80 

8.80 

8.99 

8.98 

8.74 

8.76 

8.76 

8.77 

8.51 

8.50 

4 

8.70 

8.71 

8.95 

8.95 

8.60 

8.60 

8.65 

8.65 

8.39 

8.39 

5 

8.61 

8.61 

8.90 

8.90 

8.47 

8.47 

8.54 

8.47 

8.27 

8.26 

6 

8.50 

8.51 

8.83 

8.83 

8.35 

8.35 

8.41 

8.35 

8.14 

8.14 

7 

8.40 

8.40 

8.76 

8.76 

8.24 

8.23 

8.28 

8.24 

8.02 

8.02 

8 

8.29 

8.29 

8.67 

8.67 

8.13 

8.12 

8.14 

8.15 

7.90 

7.90 

9 

8.18 

8.18 

8.58 

8.58 

8.03 

8.03 

8.00 

7.99 

7.77 

7.77 

10 

8.06 

8.07 

8.47 

8.47 

7.93 

7.94 

7.84 

7.84 

7.65 

7.65 

11 

7.94 

7.98 

8.36 

8.36 

7.83 

7.83 

7.68 

7.68 

7.53 

7.53 

12 

7.82 

7.88 

8.23 

8.23 

7.73 

7.74 

7.52 

7.52 

7.40 

7.40 

13 

7.69 

7.69 

8.10 

8.10 

7.62 

7.63 

7.34 

7.34 

7.28 

7.28 

14 

7.56 

7.55 

7.95 

7.95 

7.51 

7.51 

7.17 

7.18 

7.16 

7.16 

15 

7.43 

7.42 

7.79 

7.80 

7.39 

7.39 

6.98 

6.98 

7.04 

7.04 

16 

7.29 

7.30 

7.62 

7.63 

7.26 

7.26 

6.79 

6.81 

6.91 

6.92 

17 

7.15 

7.15 

7.45 

7.44 

7.11 

7.10 

6.60 

6.59 

6.79 

6.79 

18 

7.01 

7.00 

7.26 

7.25 

6.95 

6.93 

6.41 

6.39 

6.67 

6.66 

19 

6.86 

6.86 

7.06 

7.06 

6.77 

6.77 

6.20 

6.20 

6.54 

6.54 

20 

6.71 

6.70 

6.85 

6.88 

6.58 

6.61 

6.00 

6.04 

6.42 

6.42 

21 

6.55 

6.54 

6.63 

6.65 

6.36 

6.33 

5.79 

5.63 

6.30 

6.29 

22 

6.39 

6.38 

6.40 

6.38 

6.12 

6.07 

5.58 

5.58 

6.17 

6.17 

23 

6.23 

6.22 

6.16 

6.16 

5.85 

5.85 

5.37 

5.38 

6.05 

6.05 

24 

6.06 

6.10 

5.91 

5.97 

5.56 

5.67 

k5.15 

5.f5i 

5.93 

5.94 

25 

5.90 

5.90 

5.65 

5.54 

'5.2T' 

^^24- 

^:4.94  t 

%^5MM 

5.81 

5.80 

26 

5.72 

5.70 

5.38 

5.26 

ms4S8^ 

m§JiOi& 

5.68 

5.67 

27 

5.55 

5.55 

5.56 

5.56 

28 

5.37 

5.37 

4.Wi 

1^^^ 

5.44 

5.44 

29 

5:i8' 

5:18 

5.31 

5.30 

30 

^imoM 

fc5*fiCfe 

5.19 

5.18 

31 

S5M 

5.05 

32 

:±4M^. 

5.00 

33 


Table  l.4:The  predicted  Age  Corresponding  to  Different  States. 


Paint 
Rating 

Age                                           1 

Interstatel 

Statel 

Interstate2 

State2 

Paints 

9 

0.00 

0.00 

0.00 

0.00 

0.00 

8 

10.50 

13.66 

9.32 

9.00 

7.15 

7 

18.04 

19.30 

17.70 

14.90 

15.28 

6 

24.38 

23.66 

22.45 

20.00 

23.41 

5 

29.97 

27.35 

25.67 

24.70 

31.54 
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CHAPTER    II 

STEEL  BRIDGE  EXISTING  PAINT  LIFE  CYCLE  COST  ANALYSIS 

11.1  EXISTING  SYSTEMS  LIFE  CYCLE  COST  ANALYSIS 

//.  1. 1  Markov  Decision  Process  (MDP)  (Based  on  Stochastic  IVIettiod) 

Managers  need  to  make  decisions  regarding  selecting  the  optimal  rehabilitation 
action  for  each  paint  condition.  The  action  that  is  chosen  from  several  alternatives  affects 
the  transition  probabilities  in  the  Markov  chains  model,  as  well  as  their  immediate  and 
subsequent  costs.  The  decision  process  for  selecting  the  optimal  actions  for  the  respective 
states  of  paint  when  considering  immediate  and  subsequent  costs  is  referred  to  as  the 
Markov  Decision  Process  (MDP).  ^^^^ 

Life  cycle  cost  analysis  for  the  steel  bridges  paint  maintenance  /  Rehabilitation 
problem  can  be  studied  by  using  the  Markov  Decision  Process  (MDP).  Life  cycle  cost 
analysis  is  performed  on  the  basis  of  long-run  behavior  of  the  paint  systems.  As  plarming 
period  n  approaches  infinity  (long-run),  there  is  a  limiting  probability  (called  the 
steady-state  probability)  that  the  system  will  be  in  state  j  after  a  large  number  of 
transitions,  where  this  probability  is  independent  of  the  initial  state  i.  This  long-run 
behavior  holds  under  relatively  general  conditions  and  it  holds  for  steel  bridge  paint 
rehabilitation  problems  as  well. 

Because  the  long-run  behavior  of  a  Markov  chains  exists  for  steel  bridge  paint, 
the  rehabilitation  problem  can  be  solved  using  the  assumption  that  the  Markov  Decision 
Process  (MDP)  will  be  operating  indefinitely.  The  steps  of  that  Markov  Decision  Process 
(MDP)  can  be  summarized  as  follow:  '--^^'"^ 

1-   There  are  five  possible  states  observed  after  each  transition  ( i  =  1,2,3,4,5  )  associated 
with  steel  bridge  paint  condition  rating,  from  1  to  5,  with  5  being  the  worst.  The  index 
i  is  used  for  the  initial  state,  and  the  index  j  is  used  for  the  future  state.  Therefore, 
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state  1  corresponds  to  condition  rating  1  and  state  2  corresponds  to  condition 
rating  2  and  so  on. 
2-   After  each  observation,  a  decision  (action)  k  is  chosen  from  a  set  of  K  possible 
decisions  (k=  1,2,3,4).  Some  of  the  K  decisions  may  not  be  relevant  for  some  of  the 
states.  The  following  Table  II.  1  indicates  the  decisions  and  their  relevant  states. 

Table  11.1:  Decisions  and  relevant  states. 


Decision  k 

Action  Description 

Relevant  to  States 

1 

Do  nothing 

1,2,3,4 

2 

Spot  Repair 

2,3,4 

3 

Over-coating 

3,4 

4 

Complete  Repainting 

5 

3-   If  decision  di  =  k  is  made  in  state  i,  an  immediate  cost  is  incurred  that  has  an  expected 
value  Cik.  The  costs  in  the  following  Table  II. 2  are  the  estimated  unit  cost  ($/ton)  for 
steel  bridges. 


Table  11.2:  Estimated  Unit  Cost  of  Paint  Rehabilitation. 


Decision 

State 

Description 

Rehab.  Cost 

Cr  ( $/ton) 
Paint  Types 

(1)            (2) 

1 

1,2,3,4 

Do  nothing 

0.0           0.0 

2 

2 

Spot  Repair 

25            20 

3 

~     ~      ~ 

50            40 

4 

~     ~      .^ 

90             75 

3 

3 

Over-coating 

110           100 

4 

~     ~      ~ 

180           150 

4 

5 

Complete 
Repainting 

220           180 

Note:  $/ton(English  system)  =  1.1023  $/ton(Metric  system) 
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4-  Decision  di  =  tc  in  state  i  determines  what  the  transition  probabilities  will  be  for  the 
next  transition  from  state  i.  These  transition  probabilities  can  be  denoted  by  Pijfkj  for 
j  =  1,2,3,4,5.  The  parameter  k  in  Pij(k)  is  used  to  indicate  that  the  appropriate 
transition  probability  depends  upon  the  decision  k. 

5-  There  are  several  policies  that  can  be  used  to  make  rehabilitation  for  steel  bridge  paint 
according  to  its  state.  A  policy  is  a  set  of  decisions  (actions)  for  each  state.  Table  III. 3 
indicates  the  proposed  policies  or  scenarios  and  their  decisions  at  each  state.  '•■^■^''■^^' 

TableII.3  lists  the  relevant  policies  for  steel  bridge  paint  problems.  For  example. 
Policy  Yl:  (di,  d2,  d3,  d4,  d5)  =  (1,1,1,1,4)  describes  a  policy  where  decision  1  (do 
nothing)  is  made  in  states  1,2,3,4  and  decision  4  (complete  repainting)  is  made  in 
state  5.  The  following  transition  probability  matrix  is  assumed  to  be  the  matrix  that 
corresponds  to  Policy  Yl.  This  assumption  is  made  for  illustration  purpose  only. 


PI 


[0.80 

0.20 

0 

0 

0       ] 

[    0 

0.30 

0.50 

0.20 

0       ] 

[    0 

0 

0.1 

0.5 

0.4       ] 

[    0 

0.6 

0.4 

0 

0       ] 

[  1 

0 

0 

0 

0       ] 

Table  11.3:  Different  Policies  for  Painting  Rehabilitation. 


Decision 

Policy  (Yl) 

Policy  (Y2) 

Policy  (Y3) 

Policy  (Y4) 

Policy  (Y5) 

dl 

1 

1 

1 

1 

1 

d2 

1 

1 

1 

1 

-> 

d3 

1 

2 

3 

1 

1 

d4 

1 

1 

1 

3 

3 

ds 

4 

4 

4 

4 

4 
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The  long  run  behavior  of  PI  reached  after  seventeen  transitions.  It  is  noticed  that 
the  five  rows  have  identical  entries.  This  implies  that  the  probability  of  being  in  state  j 
after  17  periods  appears  to  be  independent  of  the  initial  paint  condition. 


0.362 

0.226 

0.191 

0.143 

0.078 

0.362 

0.226 

0.191 

0.143 

0.078 

=    [ 

0.362 

0.226 

0.191 

0.143 

0.078 

0.362 

0.226 

0.191 

0.143 

0.078 

0.362 

0.226 

0.191 

0.143 

0.078 

(PI)'" 


Because  the  long  run  behavior  of  Markov  chains  exists  for  all  five  relevant 
policies,  paint  rehabilitation  problem  can  be  solved  using  the  assumption  that  the  Markov 
decision  process  will  be  operating  indefinitely. 

6.  The  objective  is  to  find  an  optimal  Policy  that  minimizes  the  expected  total  discoimted 
cost  (immediate  cost  and  discounted  subsequent  costs  that  result  fi-om  the  future 

processes).  *^'^''*'^^^  .    ,.' 

11.1.2  Minimizing  the  Total  Expected  Discount  Cost 

Given  a  distribution  P{Xo  =  i}  over  the  initial  states  of  the  system  and  a  Policy  R, 
a  system  evolves  over  time  according  to  the  joint  effect  of  the  probabilistic  laws  of 
motion  and  sequence  of  decisions  made  (actions  taken).  In  particular,  when  the  system  is 
in  state  i  and  decision  </,■  (R)  =  A  is  made,  then  the  probability  that  the  system  is  in  state  j 
at  the  next  observed  time  period  is  given  by  pi/k).  Furthermore,  a  known  expected  cost 
C/A-is  incurred.  Denoted  by  F"/(7?)  the  expected  total  discounted  cost  of  a  system  starting 
in  state  i  (at  the  first  observed  time  period)  and  evolving  for  n  time  periods.  Then  V"i(R) 
has  two  components:  (1)  C/a,  the  cost  incurred  at  the  first  observed  time  period  as  a  result 
of  the  current  state  i  and  the  decision  rf,  (R)  =  k  when  operating  under  the  Policy  R,  and 

M 

(2)  oc  Z  Pij(k)  V"''j(R),  the  expected  total  discounted  cost  of  the  process  evolving  over 

j=0 
the  remaining  n-1  time  periods.  A  discount  factor  oc<  1  \s  specified,  so  that  the  present 

value  of  1  unit  of  cost  m  periods  in  the  future  is  cc^ .  cc  can  be  interpreted  as  equal  to 

l/(i+l),  where  i  is  the  current  interest  rate.  Thus  the  recursive  equation 
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M 

V"i(R)  =  Cik+  ocZ Pij(k)  V"'j(R),  for  i=  0,1,2,  ...,M  and  V' i(R)-  Qu  for  all  i  is  obtained. 

j=0 
This  policy  can  be  evaluated  using  the  techniques  associated  with  dynamic  programming. 

It  can  be  shown  that  as  n  approach  infinity,  this  expression  converges  to: 

M 

Vi(R)  =  Cm  +oc2:  pij(k)  Vj(R  >,  for  i=  0, 1 ,2,  . . .  ,M 

j=0 
where  Vi(R)  can  now  be  interpreted  as  the  expected  long  run  total  discounted  cost  for  a 

system  starting  in  state  i  and  continuing  indefinitely.  There  are  M+1  equations  and  M+1 

unknowns,  and  hence  Vi(R)  may  be  obtained  by  standard  methods  of  solving 

simultaneous  equations.  *'°''^^^' 

11.1.3  Policy  Improvement  Technique  Algorithm 

The  steps  for  that  algorithm  are  as  follow: 

(1)  Value  determination:  For  an  arbitrary  chosen  policy  R\,  use pij(kj)  and  Qki  to  solve 
the  set  of  M+1  equations 

M 

Vi(Rj)  =  Qui  +  ^Zpi/ki)  Vj(Rj),  for  i=  0,1,2,  ...,M 

for  all  (M+1)  unknown  values  of  Vi(Ri). 

(2)  Policy  Improvement:  Using  the  current  values  of  Vi(Rj),  find  the  alternative  policy  R2 
such  that,  for  each  state  i,  di  (R2)  =  A2  is  the  decision  that  makes 

M 

Cik2+  ocZpij(k2)  Vj(Ri),  for  i=  0,1,2,  ...,M 
j=0 

a  minimum;  that  is,  for  each  state  i.  find  the  appropriate  value  of  Aj  that 

M 
Minimizes    {Cik2  +  oc  Z  Pij(k2)  Vj(Ri )},  for  i=  0, 1 ,2,  . . .  ,M 
k2=1.2,...k  j=0 

and  then  set  rf,-  (R2)  =  minimizing  the  value  of  A?-  This  procedure  defines  a  new  policy  R2. 

If  R2  does  not  equal  to  Ri,  then  return  to  step  1  by  using  R2  instead  of  Ri,  and  solve  for 

Vi(R2),  i=  0,1,2,  ...,M.  Using  these  values,  go  to  step  (2)  and  find  R3.  Continue  in  this 
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fashion  until  you  find  two  successive  R's  to  be  equal.  When  you  find  them,  the  optimal 
policy  is  achieved,  and  the  algorithm  terminates.  In  fact,  it  can  be  shown  that: 

1-  Vi(Rj^,)^Vi(Rj)  for  i=  0,1,2,  ...,Mandj  =  1,2,.... 

2-  The  algorithm  terminates  with  the  optimal  solution  in  a  finite  number  of  iterations 

3-  The  algorithm  is  valid  without  the  assumption  that  Markov  chain  associated  with 
every  transition  matrix  is  irreducible.  ^    '"'  "' 

11.1.4  Application  of  Markov  Decision  Process  (IVIDP)  for  Steel  Bridge 
Paint 

The  previous  procedure  for  MDP  is  applied  to  steel  bridge  paint  according  to  the 
previous  proposed  policies  and  their  costs.  One  transition  probabihty  matrix  is  estimated 
for  each  4-year  period  or  range  along  paint  age.  This  transition  probability  matrix  is 
changed  or  enhanced  due  to  different  decisions  that  can  be  taken  at  each  paint  state  and 
according  to  different  policies.  Consequently,  the  procedure  that  is  explained  in  the 
previous  section  should  be  repeated  in  each  4-year  range  or  period,  for  the  same  paint 
type,  to  decide  the  optimum  rehabilitation  policy  at  each  range  of  paint  life. 

Based  on  the  transition  probability  matrix  and  the  costs  associated  with  the 
selected  policy,  five  equations  with  five  unknovras  could  be  solved  to  get  the  values  of 
the  estimated  costs  at  each  state  of  the  same  range  of  life  cycle.  Therefore,  for  the 
selected  policy,  the  procedure  for  checking  this  policy  is  applied.  The  results  of  the 
application  indicate  whether  this  policy  is  suitable  for  this  range  or  not.  If  it  is  not, 
another  policy  is  applied  and  checked  and  so  on.  If  it  is  accepted,  then  this  is  the  optimum 
policy  for  that  range.  This  procedure  is  repeated  for  each  range  of  the  paint  type  life 
cycle.  Consequently,  maintenance  plan  based  on  the  optimum  policy  can  be  planned  for 
each  paint  type  along  its  life  cycle. 

Based  on  the  estimated  costs  that  are  indicated  in  the  previous  section,  the 
optimum  policy  is  calculated  using  MDP  for  each  paint  type.  For  Interstate  Paint  Types 
(1)&(2),  State  Paint  Types  (1)&(2)  and  Paint  type  (3),  the  policy  Yl  is  the  optimum  one 
for  all  the  ranges  along  the  life  of  each  type.  This  policy  proposes  doing  nothing  for  states 
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1,2,3  and  4  and  does  complete  repainting  for  state  5.  All  the  Markov  decision  process 
(MDP)  calculations  are  shown  in  Appendix  C. 

11.1.5  Economic  Analysis  for  Different  Paint  Types 

In  the  past  few  years,  the  cost  of  maintenance  and  rehabilitation  of  existing 
coating  systems  on  steel  bridges  has  risen  dramatically.  This  increase  can  be  attributed  to 
stricter  environmental  constraints  and  higher  safety  standards  for  workers.  In  addition, 
reduced  government  funding  and  sub-optimal  maintenance  scheduling  have  contributed 
to  the  increasing  deterioration  of  these  structures.  Eventually,  corrosive  action  can  reduce 
the  cross  section  of  steel  members  and  decrease  load  capacity.  Due  to  these  problems,  a 
research  project  was  jointly  conducted  between  INDOT  and  Purdue  University  to 
investigate  current  rehabilitation  methods  to  select  the  most  economical  method.  *    ^ 

One  of  the  strategies  used  by  bridge  managers  is  to  allow  the  coating  system  to 
deteriorate  without  implementing  any  maintenance  activities.  This  option  is  certainly 
feasible  under  specific  situations.  For  example,  if  the  structure  is  near  its  expected  design 
life,  it  is  more  economic  to  replace  the  structure.  In  most  cases,  this  option  will  reduce  the 
surface  life  of  the  structure.  The  reason  for  the  decreased  life  span  is  attributed  to  the 
durability  of  the  coating  systems.  Most  coating  systems  have  life  spans  that  lie  in  the 
range  between  several  years  to  approximately  30  years,  depending  on  environmental 
conditions.  However,  as  rough  estimate  on  the  basic  service  life  of  a  bridge,  Hiroshi 
(1988)  foimd  that  the  average  life  span  of  steel  bridges  is  approximately  35  years.  There 
may  be  some  differences  between  Japanese,  Canadian,  and  U.S.  bridges,  but  the  same 
guidelines  are  used  by  all  of  them.  Therefore,  most  structures  must  require  at  least  one 
rehabilitation  activity  during  its  ser\dce  life. '"  * 

The  service  life  of  a  steel  bridge  is  limited  by  several  factors  such  as  fatigue, 
loading  capacity,  and  corrosion.  The  first  two  factors,  fatigue  and  loading  capacity', 
depend  on  loading  conditions  and  require  either  replacement  of  certain  components  or  of 
the  entire  structure  to  maintain  service.  The  third  factor,  corrosion,  can  be  controlled  by 
preventive  maintenance.  Application  of  rehabilitation  activities  can  extend  the  service  life 
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of  a  structure.  However,  the  cost  of  rehabilitation  may  sometimes  be  less  economic  than 
replacing  the  structure.  The  most  economic  choice  will  depend  on  the  costs  of  these  two 
options  and  their  relative  service  lives.  Therefore,  some  decision  support  techniques  such 
as  life-cycle  cost  analysis  can  be  used  to  help  bridge  managers  choose  the  best  strategy 
for  an  existing  or  proposed  steel  bridge.  '''*' 

Several  researchers  have  begun  to  incorporate  life  cycle  cost  analysis  ( Al-Subhi 
et  al.  1990;  Markow  1990;  McNeil  and  Finn  1987)  or  discounted  cash  flow  methods 
(Chen  and  John-ston  1990  )  to  minimize  the  cost  of  bridge  maintenance.  Weyers  et  al. 

(1988)  has  applied  this  approach  to  compare  the  equivalent  uniform  annual  cost  of 
various  rehabilitation  activities  during  the  service  life  of  a  bridge.  Rajagopal  and 
George  (1990)  have  applied  this  approach  to  road  pavements  in  order  to  reduce  their 
maintenance  costs.  In  this  research,  a  life-cycle  cost  analysis  using  equivalent  annual  cost 
is  developed  to  compare  the  cost  of  three  maintenance  strategies:  spot  repair,  over- 
coating and  complete  repainting. ''  ' 


Life  cycle  cost  analysis  is  a  relatively  simple  approach  for  minimizing  coating 
maintenance  costs.  The  general  objective  of  this  approach  is  to  determine  all  the  costs 
associated  with  the  corrosion  protection  of  the  structure  throughout  its  remaining  service 
life.  The  total  cost  for  a  combination  of  a  particular  maintenance  strategy  is  compared  to 
the  total  cost  of  another  strategy.  The  strategy  that  yields  the  lowest  cost  is  considered  to 
be  the  optimal  maintenance  strategy  for  the  specific  structure.  ^^°^ 

The  basic  components  required  in  the  life  cycle  cost  analysis  using  equivalent 
annual  cost  are:  ~ 

1-  Database:  deterioration  functions  and  maintenance  costs. 

2-  Formula  for  equivalent  annual  costs. 

3-  Correlation  between  condition  upgrades  and  rehabilitation  activities. 
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4-    Constraints:  condition  limits  for  each  activity. 

The  database  or  deterioration  curves  are  constructed  for  the  INDOT  steel  bridges 
using  the  INDOT  database.  These  functions  are  constructed  for  Interstate  and  State  paint 
types  (1)  and  (2)  and  paint  type  (3). 

The  basic  formula  used  in  the  life  cycle  cost  analysis  can  be  obtained  from  any 
basic  economics  textbook  (  Riggs  1986).  The  equation  for  equivalent  annual  cost  is: 


A  =  F(A/F,i,N)  *L=F 


(1  +  if-l 

u'here:  F  =  future  cost;  i  =  interest  rate;  A  =  annuity;  L  =  inflation  Factor  =  (1+ij'    and 
N  =  maintenance  period. 


Also,  i(l+if*L 


A=P(A/P,i,N)  *L=P 


(1+if-l 


where:  P  =  present  cost;  A  =  annual  payment. 


(20) 


These  formulas  are  used  to  calculate  the  equivalent  annual  cost  for  five  different 
categories  of  paint  deterioration  models  according  to  five  proposed  rehabilitation  and 
maintenance  alternatives.  The  proposed  alternatives  of  rehabilitation  and  maintenance  are 
as  follows:  {Note:  the  $/ton  numbers  use  the  English  system) 

1-  Do  nothing  and  do  complete  repainting  after  reaching  state  5  (approximately  30  years) 
Bridge  life  span  is  60  years. 

Oyrs  30yrs  60yrs 


T 


$220/ton  $220/ton  complete  reconstruction 
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2-  Spot  repairs  are  made  at  state  2,  occurs  every  10  years  until  the  end 
of  the  bridge  life  of  60  years. 

Oyrs        lOyrs     20yrs       30yrs     40yrs       50yrs        60yrs 


i       i       1      i      I       i       i 

$195/ton     $25/ton  $25/ton  $25/ton    $25/ton    $25/ton     reconstruction 
+  $25/ton 


3-  Spot  repairs  are  made  at  state  3,  which  are  repeated  each  18  years  until  the  end  of  the 
bridge  life. 

Oyrs  18yrs  36yrs  60yrs 


i      i 


$170/ton 
+$50/ton 


$50/ton 


$50/ton         complete  reconstruction 


4-  Over-coating  is  done  at  state  3,  repeated  each  18  years  until  the  end  of  the  bridge  life. 
Oyrs  ISyrs  36yrs  60yrs 


$110/ton 
+$110/ton 


i       i 

$  1 1 0/ton  $  1 1 0/ton         complete  reconstruction 


5-  Over-coating  is  done  at  state  4  after  the  first  24  years  and  spot  repairs  are  done  after  18 
years,  that  is  at  42nd  year,  until  the  end  of  the  bridge  life. 

Oyrs  24yrs  42yrs  60yrs 


$220/ton 


i  i 

$180/ton  $50/ton  complete  reconstruction 
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These  options  or  alternatives  are  applied  for  the  five  categories  of  paint  taking 
into  consideration  the  difference  in  cost  between  different  types  of  paint.  The  mentioned 
cost  figures  in  the  five  policies  are  based  on  paint  type  (1)  and  costs  for  other  paint  types 
are  indicated  in  Table  II. 2.  The  estimated  number  of  years  corresponding  to  bridge  states 
is  indicated  in  Table  1.3.  After  calculating  the  equivalent  annual  cost  for  each  alternative 
considering  one  paint  category,  the  minimum  cost  alternative  will  be  the  optimum  one. 
By  analyzing  the  five  categories,  it  is  concluded  that  alternative  number  two  (spot  repairs 
at  state  2  and  repetition  of  spot  repair  along  the  life  span  of  the  bridge)  is  the  optimum 
one.  This  happens  for  all  the  categories.  The  calculations  for  these  alternatives  are 
indicated  in  Figure  II.  1  to  Figure  II.5.  The  remaining  calculations  are  indicated  in 
Appendix  D. 

The  bar  chart  in  Figure  II.  I  indicates  that  the  optimum  alternative  for  Interstate  1 
is  alternative  number  2,  v/here  spot  repairs  every  10  years  should  be  done  for  the  paint 
until  the  end  of  the  bridge  life.  This  optimum  cost  is  $  15.27  /  ton  ($16.83/Metric  ton). 
This  alternative  is  less  than  half  the  cost  of  the  do  nothing  alternative  until  complete 
repainting  again  (alternative  number  1). 

The  bar  chart  in  Figure  II. 2  indicates  that  the  optimum  alternative  for  Interstate  2 
is  alternative  number  2  where  spot  repairs  every  9  years  should  be  done  until  the  end  of 
the  bridge  life.  This  optimum  cost  is  $  12.32  /  ton  ($13.58/Metric  ton).  This  alternative  is 
less  than  half  the  cost  of  the  do  nothing  alternative  until  complete  repainting  again 
(alternative  number  1). 

The  bar  chart  in  Figure  II. 3  indicates  that  the  optimum  alternative  for  State  1  is 
alternative  number  2,  where  spot  repairs  ever>'  14  years  should  be  done  until  the  end  of 
the  bridge  life.  This  optimum  cost  is  $  16.36  /  ton  ($18.03/Metric  ton).  This  alternative  is 
less  than  half  the  cost  of  the  do  nothing  alternative  until  complete  repainting  again 
(alternative  number  1). 


Figure  (II. 1):  EUAC  for  Different  alternatives  for  Interstate  1 


lEUAC 


$33.57 


$15.27 


$19.70 


$23.68 


$28.43 


Alternative  number 


Note  $/ton(English  system)  =  1.1023  $/ton(Metric  system) 


Figure  (11.2):  EUAC  for  Different  Alternatives  for  Interstate  2 


$30.00 


—    S25.00 


^     $20.00 


!£    $15.00 


Q.    $10.00 
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$5.00 


$0.00 


S26.07 

$12.32 

$16.06 

$20.04 
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$22.77 
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EEUAC 


1 


I  EUAC 


$26.07 


$12.32 


$16.06 


$20.04 


$22.77 


Alternative  number 


Note:  $/ton(English  system)  =  1.1023  $/ton(Metric  system) 
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Figure  (11.3):  EUAC  for  Different  Alternatives  for  State  1 


I  EUAC , 


EUAC        32.62790602  16.35941612  20.56072603  24.17487372  28.43286767 

Alternative  Number 


Note:  $/ton{English  system)  =  1.1023  $/ton(Metric  system) 


Figure  (11.4):  EUAC  for  Different  Alternatives  for  State  2 


S30.00 


S25.00 


w      S20.00 


B.     SI  5.00 


S10.00 


o 


S5.00 


so.oo 


I  EUAC 


I  EUAC  S25.81 


$12.32 


$15.22 


$19.92 


$23.39 


Alternative  Number 


Note:  $/ton(English  system)  =  1.1023  $/ton(Metric  system) 
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The  bar  chart  in  Figure  II.4  indicates  that  the  optimum  alternative  for  State  2  is 
alternative  number  2,  where  spot  repairs  every  9  years  should  be  done  until  the  end  of 
the  bridge  life.  This  optimum  cost  is  $  12.32  /  ton  ($13.58/Metric  ton).  This  alternative  is 
less  than  half  the  cost  of  the  do  nothing  alternative  until  complete  repainting  again 
(alternative  number  1). 

The  bar  chart  in  Figure  II. 5  indicates  that  the  optimum  alternative  for  Paint  3  is 
alternative  number  2,  where  spot  repairs  every  7  years  should  be  done  until  the  end  of  the 
bridge  life.  This  optimum  cost  is  $  15.08  /  ton  ($16.62/Metric  ton).  This  alternative  is  less 
than  half  the  cost  of  the  do  nothing  alternative  until  complete  repainting  again  (alternative 
number  1). 

The  bar  chart  in  Figure  II. 6  indicates  that  the  optimum  alternative  for  paint 
Interstate  1  &  2  is  alternative  number  2  for  Interstate  2,  where  spot  repairs  every  9  years 
should  be  done  until  the  end  of  the  bridge  life.  This  optimum  cost  is  $  12.32  /  ton 
($13.58/Metric  ton).  This  alternative  is  less  than  half  the  cost  of  the  do  nothing  alternative 
until  complete  repainting  again  (alternative  number  1). 


The  bar  chart  in  Figure  II. 7  indicates  that  the  optimum  alternative  for  paint  State  1 
&  2  is  alternative  number  2  for  State  2,  where  spot  repairs  every  7  years  should  be  done 
until  the  end  of  the  bridge  life.  This  optimum  cost  is  $  12.32  /  ton  ($13.58/Metric  ton). 
This  alternative  is  less  than  half  the  cost  of  the  do  nothing  alternative  until  complete 
repainting  again  (alternative  number  1). 
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11.2  COMPARISON  OF  EXISTING  WITH  NEW  SYSTEMS 

According  to  life  cycle  cost  analysis,  ZincA^inyl  is  clearly  the  best  existing  paint 
system.  Therefore,  ZincA^inyl  should  be  compared  with  some  new  systems  to  decide 
whether  it  is  truly  the  best  system  for  the  future.  This  comparison  is  only  fi:om  the  cost 
basis.  The  available  information  from  other  states  (Illinois,  Michigan,  Ohio,  Kentucky 
and  Connecticut)  about  new  paint  systems  uses  different  units  of  measurement.  Other 
states  use  $  per  square  foot  as  a  cost  reference,  instead  of  using  $  per  ton.  INDOT  uses  $ 
per  ton  as  a  cost  reference;  therefore  a  conversion  factor  is  very  important  to  compare 
these  different  cost  references.  Our  study  relies  on  the  other  states'  experience  to  evaluate 
some  new  systems  that  did  well  in  these  states.  Consequently,  a  conversion  factor  should 
be  used  to  compare  ZincA^inyl  ($/ton)  with  the  other  new  systems  from  other  states 
($/fr).  Table  II. 4  shows  the  calculation  of  the  conversion  factor  from  a  sample  of  four 
bridges.  This  information  is  based  on  data  from  INDOT  and  some  analysis  from  Mr. 
Maged  Georgy's  report.  Mr.  Ted  Hoppwood*-^^^  of  KentuckyDOT  said  that  the  conversion 
factor  is  approximately  $/ton  =  125  $/  ft^  ($/ton  =  1344$/m~),  based  on  his  experience. 
Therefore,  on  average,  a  conversion  factor  of  125  can  be  used  in  our  study  until  the 
calculation  of  the  accurate  number  for  INDOT  is  possible.  The  new  systems  that  can  be 
studied  here  and  that  are  suitable  to  be  used  by  INDOT  are  a  3 -Coat  System  and 
Metalization.  It  is  proposed  to  compare  these  new  systems  with  the  existing  ZincA/'inyl 
system. 

Based  on  the  information  that  is  included  in  Table  II.4,  the  ZincA^inyl  system  cost 
in  S/ft"  can  be  calculated  as  the  average  of  the  two  available  bridges.  This  calculation 
concludes  a  cost  of  $2.5/6^^  ($26.88/m^)  for  this  system.  Based  on  data  from  INDOT  for 
3 -coat  system,  Inorganic/Organic  Zinc  Epoxy  Urethane  (In/OZEU),  the  unit  cost  per 
square  foot  is  calculated  for  that  system.  This  unit  cost  is  calculated  taking  into 
consideration  two  different  old  paint  systems.  The  first  application  is  done  after  the 
removal  of  Lead  based  paint  as  an  old  system.  The  second  application  is  done  after  the 
removal  of  Zinc  based  paint  as  an  old  system.  Table  II. 5  shows  the  calculation  of  unit 
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Table  11.5:  3C0AT  Sytem  over  Lead  Old  paint 


Original 
paint 

New 
paint 

Cost  /WIG 

Cost/ton 

Cost/ft' 

Lead 

3C0AT 

207 

187.7871 

1.632931 

Lead 

3COAT 

212 

192.323 

1.672374 

Lead 

3C0AT 

394 

357.4305 

3.108091 

Lead 

3C0AT 

143 

129.7273 

1.128064 

Lead 

3C0AT 

324 

293.9276 

2.555892 

Lead 

3COAT 

fUmffis^ 

l^l^P 

W'iiiiyiHiHli^ 

Lead 

3COAT 

218 

197.7661 

i:7fl9705 

Lead 

3COAT 

556 

504.3943 

4.386037 

Lead 

3C0AT 

427 

387.3676 

3.368414 

Lead 

3COAT 

421 

381.9245 

3.321082 

Lead 

3COAT 

505 

458.1279 

3.983721 

Lead 

3C0AT 

562 

509.8374 

4.433369 

Lead 

3C0AT 

562 

509.8374 

4.433369 

Lead 

3C0AT 

585 

530.7026 

4.614806 

Lead 

3C0AT 

585 

530.7026 

4.614806 

Lead 

3C0AT 

563 

510.7446 

4.441257 

Lead 

3C0AT 

563 

510.7446 

4.441257 

Lead 

3C0AT 

556 

504.3943 

4.386037 

Lead 

3COAT 

469 

425.4693 

3.699733 

Lead 

3COAT 

480 

435.4483 

3.786507 

Lead 

3C0AT 

480 

435.4483 

3.786507 

Average 

453.0476 

410.9976 

3.573892 

StDev 

152.0182 

137.9085 

1.199205 

Median 

480 

435.4483 

3.786507 

Confidence 

66.62361 

60.43987 

0.525564 

mmm 

HHHHRHIHH 

Average 

502 

455.4064 

3.960055 

StDev 

78.25258 

70.98949 

0.6173 

Median 

530.5 

481.2611 

4.184879 

Confidence 

39.60048 

35.92493 

0.312391 

Note:$/n''  =  10.76  $/m'^ 
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Table  11.6:  3COAT  Sytem  over  ZINC  Old  paint 


Original 
paint 

New 
paint 

Cost  /WIG 

Cost/ton 

Cost/ft^ 

ZINC 

3COAT 

310 

281.227 

2.445453 

ZINC 

3COAT 

338 

306.6282 

2.666332 

ZINC 

3C0AT 

338 

306.6282 

2.666332 

ZINC 

3C0AT 

342 

310.2569 

2,697886 

ZINC 

3C0AT 

268 

243.1253 

2.114133 

ZINC 

3C0AT 

465 

421.8406 

3.668179 

ZINC 

3C0AT 

465 

421.8406 

3.668179 

ZINC 

3COAT 

419 

380.1101 

3.305305 

ZINC 

3COAT 

337 

305.721 

2.658444 

ZINC 

3COAT 

193 

175.0865 

1.522491 

ZINC 

3COAT 

248 

224.9816 

1 .956362 

ZINC 

3C0AT 

248 

224.9816 

1.956362 

ZINC 

3COAT 

248 

224.9816 

1.956362 

ZINC 

3C0AT 

248 

224.9816 

1 .956362 

ZINC 

3C0AT 

317 

287.5773 

2.500672 

ZINC 

3C0AT 

299 

271.248 

2.358678 

ZINC 

3C0AT 

317 

287.5773 

2.500672 

ZINC 

3C0AT 

312 

283.0414 

2.46123 

ZINC 

3C0AT 

287 

260.3618 

2.264016 

ZINC 

3COAT 

459 

416.3975 

3.620847 

ZINC 

3COAT 

479 

434.541 1 

3.778619 

ZINC 

3C0AT 

479 

434.5411 

3.778619 

ZINC 

3COAT 

422 

382.8316 

3.328971 

ZINC 

3COAT 

531 

481'r7"f47 

4r1888Z4 

ZINC 

3COAT 

s^      531 

481.7147 

4.188824 

ZINC 

3C0AT 

413 

374.657 

3.257974 

Average 

358.1923 

324.9463 

2.82562 

StDev 

96.87229 

87.881 

0.764183 

Median 

337.5 

306.1746 

2.662388 

Confidence 

37.97319 

34.44867 

0.299554 

■■■ill 

I^^^^H 

Average 

304.75 

276.4643 

2.763739 

StDev 

46.6569 

42.3264 

0.643022 

Median 

311 

282.1342 

2.658444 

Confidence 

23.61118 

21.41969 

0.325408 

A/ote:$/ft^  =  10.76  $/m^ 
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cost  of  that  system  application  over  an  old  Lead  based  paint.  It  is  $  3.96/ft~  ($42.58/m^) 
on  average.  Table  II. 6  shows  the  calculation  of  unit  cost  for  that  system  application  over 
an  old  Zinc  based  paint.  It  is  $2.80/ft^  ($30.1 1/m^)  on  average. 

Data  was  collected  from  different  states  for  different  paint  types.  Table  II. 7  shows 
the  cost  data  that  is  collected  and  calculated  for  different  paint  systems  from  different 
departments  of  transportation.  ZincA'^inyl  paint  system  cost  for  INDOT  is  calculated  from 
Table  II.4.  The  3-coat  system  cost  for  INDOT  is  calculated  as  in  Tables  II. 5  and  II. 6. 
Data  was  collected  from  other  state  DOTs,  such  as  ODOT.  MOOT  and  ILDOT,  for  that 
paint  system.  Metalization  data  was  obtained  from  Maged  Georgy's  ^^^  report  for  INDOT. 
Other  data  was  collected  from  ILDOT  and  CTDOT  for  the  cost  per  square  foot  of 
metalization.  Based  on  report  no.  FHWA-RD-96-058,  metalization  cost  per  square  foot  is 
taken  to  be  included  in  this  study. 

Based  on  the  collection  of  cost  data  for  different  paint  systems,  an  economic 
analysis  is  done  to  compare  these  different  systems.  Table  II. 8  shows  the  present  value 
(PV)  calculation  for  the  ZincA^inyl  paint  system.  It  is  $5.50  /  ft^.  The  EUAC  calculated 
for  the  same  paint  type  is  $0.39  /  ft" .  Table  II. 9  shows  the  calculation  of  PV  for  a  3-coat 
system  where  the  old  paint  is  Lead  or  Zinc.  Two  options  are  chosen  for  calculation  based 
on  20  and  25  years  of  service  life.  The  PV  and  EUAC  of  that  paint  type,  based  on  these 
two  different  options,  are  indicated  in  Table  II. 9.  Table  11.10  shows  the  PV  and  EUAC 
calculation  for  metalization.  This  table  shows  that  PV  =  $18.91/  ft^  and  EUAC  =  $1,347/ 
ft"  for  metalization.  Figure  II. 8  shows  the  values  of  PV  for  different  paint  systems.  It 
indicates  that  the  minimum  PV  value  is  for  3-coat  paint  system,  which  is  $4.55/  ft^. 
Figure  II. 9  shows  that  the  EUAC  for  that  paint  system  is  the  minimum,  which  is  $0,324/ 
ft'.  Therefore,  a  conclusion  can  be  made  that  the  best  paint  system  is  the  3-coat  system 
according  to  the  cost  point  of  view. 
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Table  11.7:  Cost  /  ft  and  Service  Life  for  Different  Paint  Systems 
in  Different  States. 


\pors 

Paint  Sys.  \ 

Cost  /  ft^  and  Service  Life  for  Different  Paint  Systems 

INDOT 

ODOT 

MOOT 

ILDOT 

KDOT 

CTDOT 

FHWA 

Service  Life 
(Years) 

Zinc-Vinyl 

$2.50 

X 

X 

X 

X 

X 

X 

15-25 

3Coat/Lead 

$3.96 

$4.0-6.0 

$9.29 

$5 

X 

X 

X 

25-30 

3Coat/Zinc 

$2.80 

$4.0-6.0 

$9.29 

$5 

X 

X 

X 

25-30 

Metalization 

$16.81 

X 

X 

$6.0-9.0* 

X 

$12.0-15.0 

$^4.75 

40-60 

INDOT  Numbers:  calculated  from  data.  Georgy  and  Chang  '^' 
ODOT  Numbers:  collected  from  Mr.  Herald  Schultz'^^'  and  Mr.  R.  Bauer.'^^^ 
MDOT Numbers:  collected  from  Mr.  Sonny  Gduan'^^^ ,  Mr  Brian  Back'^°' , 
Mr  Craig  A.  Russell'^^'  and  Mr  Glenn  Bukosi<y.'^^^ 

ILDOT  Numbers:  collected  from  Mr  Gary  Kowalski.'^^' 
CTDOT  Numbers:  collected  from  Mr  Eric  Lohrey.  '^"^ 
FHWA  Numbers:  collected  from  report  No.  FHWA-RD-96-058.  '"''*'' 
Service  Life  Numbers:  are  collected  from: 

Zinc-Vinyl  from  Mr.  Ted  Hopwood^^^^ 

3Coat-System  from  ILDOT  and  MDOT. 

Metalization  from  ILDOT  and  CTDOT. 


*  Tfiis  number  is  for  new  bridges  only. 
A/ote:$/ft^=  10.76  $/m^ 
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Table  11.8:  ZincA/inyl  Economic  Analysis 


Age 
(yrs) 

Cost 
initial/ft^ 

Interest  Rate 

ref  no.  (32) 

15 

$2.50 

7.00% 

inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost: 

0 

1 

$2.50 

Present  Value  for  Fv  after  15  yrs: 

15 

1.675349 

$1.52 

Present  Value  for  Fv  after  30  yrs: 

30 

2.806794 

$0.92 

Present  Value  for  Fv  after  45  yrs: 

45 

4.702359 

$0.56 

Total  Present  Value  TPV/ft^ 
EUAC 

$5.50 
$0.39 

Note:$fn^  =  10.76  $/m^ 
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Table  11.9:  3-Coat  Sytem  (OZEU)  Economic  Analysis 


Age 
(yrs) 

Cost 

initial/ft^ 

Interest  Rate 

Old  Paint 

Option 

25 
25 

20 
20 

$4.00 
$2.80 

$4.00 
$2.80 

700% 
7.00% 

7.00%, 
7.00% 

Lead 
Zinc 

Lead 
Zinc 

1 
1 

2 
2 

Option  No.  1 

0.035 

Inflation  Factor 

item 

n 

Infl. Factor 

Lead  Cost 

Zinc  Cost 

Present  Value  for  initial  Cost 

0 

1 

$4.00 

2.8 

Present  Value  for  Fv  after  25  yrs 

25 

2.363245 

$1.74 

1.219193 

Present  Value  for  Fv  after  50  yrs 

50 

5.584927 

$0.76 

0.530868 

Total  Present  Value  TPV/ft^ 
EUAC 

$6.50 
$0.46 

$4.55 
$0.32 

Option  No.  2 

0.035 

Inflation  Factor 

Item 

n 

Infl. Factor 

CostLead 

Cost  Zinc 

Present  Value  for  initial  Cost 

0 

1 

$4.00 

2.8 

Present  Value  for  Fv  after  20  yrs 

20 

1.989789 

$2.06 

1.439758 

Present  Value  for  Fv  after  40  yrs 

40 

3.95926 

$1.06 

0.740322 

Total  Present  Value  TPV/ft^ 
EUAC 

$7.11 
$0.51 

$4.98 
$0.35 

Note:$lft  =  10.76$/m^ 
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Table  11.10:  Metalization  Economic  Analysis 


Met.Age 

(yrs) 

Met.Cost 
initial 

Coat  Age 

(yrs) 

Coat  Cost 
initial 

Coat  Cost 
Maintenance 

Coat  Type 

Interest  Rate 

60 

$14.75 

10 

$0.42 

$1.82 

Epoxy/Ur 

7.00% 

Inflation  Factor                                        0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost 
Present  Value  for  Fv  after  10  yrs 
Present  Value  for  Fv  after  20  yrs 
Present  Value  for  Fv  after  30  yrs 
Present  Value  for  Fv  after  40  yrs 
Present  Value  for  Fv  after  50  yrs 

0 

10 

20 

30 

40 

50 

1 
1.410598761 
1.989788863 
2.806793705 
3.959259721 
5.584926856 

$15.17 
$1.31 
$0.94 
$0.67 
$0.48 
$0.35 

Total  Present  Value  TPV/ft^ 
EUAC 

$18.91 
$1.35 

Note-.m  =  10.76  $/m^ 
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Figure  11.8:  Present  Value  ($/ft2)  of  Different  Paint  Systems 
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Figure  11.9:  EUAC  ($/ft2)  of  Different  Paint  Systems 
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CHAPTER   III 

STEEL  BRIDGE  NEW  PAINT  SYSTEM  COST  ANALYSIS 
(3-Coat  Paint  System:  In/OZEU) 

III.1   3-COAT  SYSTEM  MICHIGAN  DATA  ANALYSIS 

Based  on  Michigan's  experience  with  3-coat  paint  system,  data  was  collected 
from  MDOT  to  analyze  that  system.  A  regression  function  was  fit  to  this  data  and  the 
service  life  of  that  paint  system  was  calculated.  This  analysis  was  done  to  support  the 
economic  analysis  and  service  life  estimate  for  the  3-coat  system.  SAS  was  used  to  fit  the 
model  to  the  data.  Table  III.l  shows  the  different  models  that  can  fit  the  data. 

///.  1. 1  Regression  Analysis 

For  this  paint  type  (3-coat  system),  the  performance  function  is  fitted  to  linear, 
quadratic  and  cubic  formulas.  The  r-squared  value  for  each  formula  is  calculated  and 
indicated  in  Table  III.  1 ,  but  the  difference  among  them  is  very  small.  The  lack-of-fit  test 
was  done  for  each  case  and  it  is  not  significant  for  any  formula.  This  means  that  all  the 
models  are  adequate  to  fit  the  data.  The  assumptions  for  the  regression  models  are 
checked  for  the  error  constant  variance,  normality  and  independence.  The  results  indicate 
that  all  the  models  are  good  according  to  these  assumptions,  as  indicated  in  Appendix  A. 
In  addition  to  this  information,  the  significance  of  all  models'  parameters  is  checked  to 
select  the  best  model  that  can  fit  the  data.  Based  on  these  tests  and  arguments,  the  best 
model  is  selected  for  3-coat  system  paint  is  a  linear  model:  *'^' 


Paint  Rating  =  8.08  -  0.  116  *  Age. 


III. 1.2  Markov  Chains  Analysis 

Markov  chains  process  can  be  applied  for  the  data  available  to  draw  the 
probabilistic  deterioration  curve  for  a  3-coat  system.  The  same  steps  for  applying  Marko\" 
chains 
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process  that  was  used  in  the  existing  system  analysis  was  used  for  the  3-coat  system.  The 
steps  are: 
1  -    The  paint  age  is  divided  into  4-year  ranges  with  four  unknown  probabilities  in  each 
range. 

2-  The  Maple  mathematics  program  is  used  to  prepare  the  file  of  matrix  multiplication 
as  illustrated  in  the  previous  section.  The  output  of  this  file  are  the  equations  that 
contain  p(l),  p(2),  p(3)  and  p(4)  as  unknowns.  Therefore,  there  are  four  equations 
with  four  unknowns  that  can  be  solved  to  get  the  values  of  p(l),  p(2),  p(3)  and  p(4). 

3-  The  output  of  the  previous  step  is  substituted  into  the  non-linear  programming 
equation  as  E(t,P),  which  is  subtracted  from  the  part  Y(t).  This  part,  Y(t),  is 
calculated  using  the  regression  formulas  for  the  various  models'  types. 

4-  After  constructing  these  equations,  they  are  put  into  a  GAMS  program  input  file  to 
solve  the  NLP  problem  based  on  the  constraints  that  are  illustrated  in  the  previous 
section.  This  input  file  allows  GAMS  to  minimize  the  absolute  value  of  Y(t)-E(t,P) 
for  each  four  years  of  paint  age.  The  outputs  of  this  step  are  the  values  p(l),  p(2), 
p(3)  and  p(4)  that  minimize  the  NLP  objective  function. 

5-  By  knowing  the  values  of  these  probabilities  in  each  4-year  range,  the  last  year  state 
vector  for  each  range  can  be  calculated  by  using  the  Maple  mathematics  program. 
This  state  vector  is  calculated  to  be  used  as  the  initial  state  vector  for  the  next  range 
of  four  years.  .,  ■•' 

6-  Steps  1  to  5  are  repeated  for  each  4-  year  range  of  the  3-coat  system  paint 

age    (12)  (13)  (23)  (24)  (25) 

After  applying  this  procedure,  transition  probability  matrices  are  constructed  for 
each  4-year  range  of  the  3-coat  system  paint  age.  The  values  of  p(l),  p(2),  p(3)  and  p(4) 
for  each  transition  probability  matrix  are  shovm  in  Appendix  B. 

III. 1.3  Regression  vs.  Markov  Process  Results  for  3-Coat  System 

The  regression  model  provides  the  values  of  paint  rating  at  various  ages  as  shown 
in  Table  in.2.  This  table  indicates  the  paint  rating  values  for  different  ages  using 
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Table  III.2:  Michigan  Paint  Rating  Prediction 

Using  Regression  and  IVIarkov. 


fRaint  RaLirng|H|| 

Mich.  SGOA'H 

iReg    . 

JBi^^M 

0 

8.0819 

8.0819 

1 

7.9657 

8.0001 

2 

7.8496 

7.856 

3 

7.7334 

7.73 

4 

7.6173 

7.618 

5 

7.5011 

7.516 

6 

7.3849 

7.425 

7 

7.2688 

7.341 

8 

7.1526 

7.264 

9 

7.0365 

7.034 

10 

6.9203 

6.918 

11 

6.8041 

6.803 

12 

6.688 

6.692 

13 

6.5718 

6.57 

14 

6.4557 

6.453 

15 

6.3395 

6.34 

16 

6.2233 

6.234 

17 

6.1072 

6.105 

18 

5.991 

5.985 

19 

5.8749 

5.8751 

20 

5.7587 

5.775 

21 

5.6425 

5.632 

22 

5.5264 

5.511 

23 

5.4102 

5.41 

24 

5.2941 

5.33 

25 

5.1779 

5.15 

26 

5.0617 

5.062 

27 

4.9456 

5.0085 

28 

4.8294 

29 

4.7133 

30 

4.5971 

Table  111.1:  3-COAT-System 
from  Michigan 


Age 

Paint  Rating                | 

1st  order 

2nd  order 

3rd  order 

5 

7.5 

7.68675 

7.67153 

10 

6.92 

6.77325 

6.75053 

15 

6.34 

5.29475 

3.88953 

20 

5.76 

3.25125 

-2.30647 

25 

5.18 

0.64275 

-13.2325 

30 

4.6 

-Z53075 

-30.2835 

■■■BB 

WKBM 

R-square 

0.438 

0.461 

0.465 

Parameters  t-tes 

O.K. 

N.O.K. 

N.O.K. 

Model  F-test 

O.K. 

O.K. 

O.K. 

Normality  Assur 
Lack  of  Fit  test 

J     O.K. 
O.K. 

O.K. 
O.K. 

O.K. 
O.K. 

Selection       ^jgggjggg^ 

No 

No 
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regression  and  the  Markov  chains  process.  Consequently,  the  deterioration  curve  for  that 
type  of  paint  can  be  drawn  using  the  values  of  condition  rating  calculated  by  using  the 
Markov  chains  and  regression.  Figure  III.l  shows  these  two  curves.  This  figure  indicates 
that  the  deviation  between  the  two  curves  is  very  small  and  they  look  the  same.  As 
known,  regression  is  based  on  Least  Squared  Method,  which  gives  the  best  fit  for  the 
data.  Therefore,  the  probabilistic  model,  Markov  chains  model,  is  good  for  fitting  this 
data.  Consequently,  the  deterministic  and  probabilistic  fitted  curves  for  this  data  are 
constructed  and  ready  to  be  used  in  any  applications.  Table  III.2  shows  that  at  a  paint 
condition  rating  of  5,  the  age  of  the  3-coat  system  will  be  between  26  and  27  years.  This 
number  supports  the  25  years  paint  age  estimate. 

111.1. 4  Markov  Decision  Process  (MDP)  Application  for3-Coat  System 

The  steps  in  applying  the  Markov  Decision  Process  (MDP)  to  the  3-coat  system 
are  the  same  as  that  of  the  existing  system.  These  steps  are: 

1-  There  are  five  possible  states  (  i  =  1,2,3,4,5  )  observed  after  each  transition 
associated  with  steel  bridge  paint  condition  rating,  from  1  to  5,  with  5  being  the 
worst.  The  index  i  is  used  for  the  initial  state,  and  the  index  j  is  used  for  the  future 
state.  Therefore,  state  1  corresponds  to  condition  rating  1  and  state  2  corresponds  to 
condition  rating  2  and  so  on. 

2-  After  each  observation,  a  decision  (action)  k  is  chosen  from  a  set  of  K  possible 
decisions  (k=  1,2,3,4).  Some  of  the  K  decisions  may  not  be  relevant  for  some  of  the 
states.  Table  III. 3  indicates  the  decisions  and  their  relevant  states. 

Table  III. 3  :  Decisions  and  relevant  states. 


Decision  k 

Action  Description 

Relevant  to  States 

I 

Do  nothing 

1,2,3,4 

2 

Spot  Repair 

3,4 

3 

Complete  Repainting 

5 
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3-    If  decision  di  =  k  is  made  in  state  i,  an  immediate  cost  is  incurred  that  has  an  expected 
value  Cik.  The  costs  in  Table  III.4  are  the  estimated  unit  cost  ($/ton)  for  steel  bridges. 

Table  1 1 1.4  :  Estimated  Unit  Cost  of  3-Coat  Paint  Rehabilitation. 


Rehab.  Cost 

Disruption 

Total  Cost* 

Decision 

State 

Description 

Cr  (  $/ft^) 

Cost  (Cd) 

Cik(  $/ft^) 

Lead    Zinc 

Lead  Zinc 

Lead  Zinc 

1 

1,2,3,4 

Do  nothing 

0.0         0.0 

0.0        0.0 

0.0       0.0 

2 

3 

Spot  Repair 

1.5          1.0 

17.69  17.69 

19.2     18.69 

4 

~     ~      ~ 

2.5          2.0 

26.53  26.53 

29.03  28.53 

3 

5 

Complete 
Repainting 

4.0         2.8 

35.37  35.37 

39.37  38.17 

Note:  $/i 

t"=  10.76 

$/m- 

*  Total  Cost  includes  the  rehabilitation  cost  and  the  disruption  cost  to  travelers  due  to 
the  rehabilitation  work.  Disruption  Cost  Calculations  are  included  in  Appendix  D. 

4-  Decision  di  =  k  in  state  i  determines  what  the  transition  probabilities  will  be  for  the 
next  transition  from  state  i.  These  transition  probabilities  can  be  denoted  by  Pij(k)  for 
j  =  1,2,3,4,5.  The  parameter  k  in  Pij(k)  is  used  to  indicate  that  the  appropriate 
transition  probability  depends  upon  the  decision  k. 

5-  There  are  several  policies  that  can  be  used  to  rehabilitate  the  steel  bridge  paint  based 
on  its  status.  A  policy  is  a  set  of  decisions  (actions)  for  each  state.  Table  in.3 
indicates  the  proposed  policies  or  scenarios  and  their  decisions  at  each  state.  *^^'  ^'°^ 

Table  III.  14  lists  the  relevant  policies  for  steel  bridge  paint  problems.  For 
example,  Policy  Yl:  (di,  d2,  ds,  d4,  ds)  =  (1.1,1,1,4)  describes  a  policy  where  decision  1 
(do  nothing)  is  made  in  states  1,2,3,4  and  decision  4  (complete  repainting)  is  made  in 
state  5. 


6.  The  objective  is  to  fmd  an  optimal  policy  that  minimizes  the  expected  total  discounted 


LIFE  CYCLE  COST  ANAL  YSIS  AND  MAINTENANCE  PLAN 


72 


cost  (immediate  cost  and  discounted  subsequent  costs  that  result  from  the  future 
processes).  ^''^'^^^^ 

Table  1 1 1. 5  :  Different  Policies  for  Painting  Rehabilitation. 


Decision 

Policy  (Yl) 

Policy  (Y2) 

Policy  (Y3) 

di 

1 

/ 

d2 

1 

1 

d3 

2 

1 

d4 

1 

2 

ds 

3 

3 

After  applying  these  steps  to  the  3 -coat  system,  the  calculations  indicate  that 
scenario  or  policy  Yl  is  the  most  economic  one.  This  policy  says  do  nothing  until  the 
bridge  reaches  state  5  and  then  do  a  complete  repaint  for  the  bridge.  All  the  calculations 
are  included  in  Appendix  C. 

III. 1.5  Economic  Analysis  for  3-Coat  System 

The  economic  formulas  are  used  to  calculate  the  present  value  (PV)  and  the 
equivalent  uniform  annual  cost  (EUAC)  for  the  3-coat  system  paint  according  to  three 
proposed  rehabilitation  and  maintenance  alternatives.  The  proposed  alternatives  of 
rehabilitation  and  maintenance  are  as  follows:  (Note:  $/ft"  =  10.76  $/m") 


1-  Do  nothing  and  do  complete  repainting  after  reaching  state  5  (approximately  25  years). 
Bridge  life  span  is  60  years. 

Oyrs  25yrs  50yrs      60yrs 


$4.0/ft^ 


i 


$4.0/ft^ 


1 


$4.0/ft~   reconstruction 
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2-  Spot  repairs  are  made  at  state  3  (rating  7),  occurs  every  15  years  until  the  end 
of  the  bridge  life  of  60  years. 


Oyrs 


15yrs 


3  Oyrs 


45yrs  60yrs 


i  i  i  i  i 

$4.0/ft-  S1.5/ft^  $1.5/ft^  $1.5/ft^   complete  reconstruction 


3-  Spot  repairs  are  made  at  state  4  (rating  6),  which  occurs  every  20  years  until  the  end  of 
the  bridge  life. 

Oyrs  20yrs  40yrs  60yrs 


T 

$4.0/ft^ 


i  i 


$2.5/ft^ 


$2.5/ft       complete  reconstruction 


These  options  or  alternatives  are  applied  for  the  3 -coat  system.  The  mentioned 
cost  numbers  in  the  three  policies  are  based  on  the  3-coat  system  applied  over  Lead  as  a 
base  for  old  paint  and  costs  for  that  paint  over  Zinc  old  paint  are  indicated  in  Table  III.4. 
The  estimated  number  of  years  corresponding  to  bridge  states  are  indicated  in  Table  III.  1 . 
After  calculating  the  PV  and  EUAC  for  each  alternative  considering  one  paint  category, 
the  minimum  cost  alternative  will  be  the  optimvim  one.  By  analyzing  the  five  categories, 
it  is  concluded  that  alternative  number  two  (spot  repairs  at  state  3  and  repetition  of  spot 
repairs  over  the  life  span  of  the  bridge)  is  the  optimum  one.  This  happens  for  all  the 
categories.  The  calculation  for  these  alternatives  PV  and  EUAC,  for  3-coat  over  Lead  as 
an  old  paint,  are  indicated  in  Figures  111.2  and  III. 3.  The  calculations  of  PV  and  EUAC 
for  the  3-coat  over  Zinc  old  paint  are  indicated  in  Figures  III. 4  and  III.5. 


The  bar  chart  in  Figure  III. 2  indicates  that  the  optimum  alternative  for  the  3-coat 
system  is  alternative  number  2,  where  spot  repairs  every  15  years  should  be  done.  This 
optimum  cost  is  PV  =  $  5.8/ft^  ($62.4 1/m^).  This  alternative  is  approximately  two  third 
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the  cost  of  the  do  nothing  alternative  until  complete  repainting  again  (alternative  number 
1).  Figure  III. 3  shows  that  the  best  policy  is  alternative  no.  2  with  EUAC  =  $0.413/ft^ 
($4.44/m-). 

The  bar  chart  in  Figure  III. 4  indicates  that  the  optimum  alternative  for  the  3 -coat 
system  is  alternative  number  2,  where  spot  repairs  every  1 5  years  should  be.  This 
optimum  cost  is  PV  =  $  4.0/ft  ($43.04/m  ).  This  alternative  is  approximately  two  third 
the  cost  of  the  do  nothing  alternative  until  complete  repainting  again  (alternative  number 
1).  Figvire  III. 5  shows  that  the  best  policy  is  alternative  no.  2  with  EUAC  =  $0.285/ft^ 
($3.07W). 

By  adding  the  disruption  cost  to  the  paint  cost,  the  best  scenario  does  not  change 
as  indicated  in  Figures  from  III. 7  to  III.  10  for  3-coat  over  old  paint  Lead  based  or  Zinc 
based.  The  best  scenario  is  making  spot  repairs  every  1 5  years  and/or  the  bridge  reaches 
condition  rating  of  7.  All  disruption  cost  calculations  are  included  in  Appendix  D. 

///.  1.6  SENSITIVITY  ANAL  YSIS 

Sensitivity  analysis  was  carried  out  for  the  3-coat  system  to  check  the  best 
scenario  to  paint  age  and  cost  changes.  Table  III.  6  shows  the  final  results  of  trial  and 
error  method  to  get  the  minimum  paint  age  where  spot  repairs  remain  the  best  scenario. 
The  values  of  PV  and  EUAC  for  spot  repairs  at  state  3  scenario  look  close  to  that  of  spot 
repairs  at  state  4  scenario.  Consequently,  conclusion  can  be  made  that  spot  repairs  at  state 
3  scenario  will  remain  the  best  scenario  if  the  bridge  reaches  state  3  or  condition  rating  7 
after  12  years  old.  Table  III. 6  shows  this  conclusion  for  the  3-coat  system  applied  over 
Zinc  as  an  old  paint.  Table  III. 7  shows  this  conclusion  for  the  3-coat  system  applied  o\er 
Lead  as  an  old  paint.  The  initial  cost  and  spot  repairs  cost  were  changed  by  trial  and  error 
to  check  the  maximum  cost  value  that  maintain  spot  repairs  at  state  3  the  best  scenario. 
The  results  of  this  sensitivity  analysis  are  shown  in  Table  III.8  for  3-coat/Zinc  and  in 
Table  III. 9  for  3-coat/Lead.  Table  III. 8  shows  that  the  3-coat  spot  repairs  at  state  3  \\ill  be 
the  best  scenario  if  the  initial  cost  is  less  than  $8.85/ft'  and  the  cost  e\erv  15  vears  is  less 
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than  $3.16/ft^  ($34.00W)  for  3-coat/Zinc.  Table  III.9  shows  that  the  3-coat  spot  repairs 
at  state  3  will  be  the  best  scenario  if  the  initial  cost  is  less  than  $9.8/ft^($105.45/m^)  and 
the  cost  every  15  years  is  less  than  $3.68/ft^  ($39.60/m")  for  3-coat/Lead. 

III.1.7  Conclusion 

There  is  a  difference  between  the  results  of  the  economic  analysis  and  the  results 
of  Markov  decision  process  (MDP).  The  economic  analysis  method  uses  deterministic 
values  of  cost  and  time  where  MDP  uses  probabilistic  times.  Therefore,  research  was 
done  to  investigate  this  conflict  between  the  results  of  these  two  methods.  The  literature 
of  probability  theory  applications  and  its  relation  with  the  deterministic  theory 
applications  shows  that  it  is  not  necessarily  getting  the  same  results  from  both  methods 
for  the  same  application.  For  example,  in  Industrial  Engineering,  the  TSP  (Traveling 
Salesperson  Problem)  can  be  solved  using  deterministic  and  probabilistic  PTSP 
(Probabilistic  Traveling  Salesperson  Problem)  values.  The  result  of  this  application  is 
shovra  in  the  following  paragraph. 

Bertsimas  (1990)  assigned  the  problem  of  how  heuristic  approaches  to  the 
deterministic  problem  TSP  perform  when  applied  to  the  corresponding  probabilistic 
problem  PTSP.  This  literature  suggested  that  if  the  coverage  probability  p  was  large,  the 
constant  guarantee  heuristics  for  the  deterministic  problem  still  behave  well  for  the 
corresponding  probabilistic  problem.  But  if  p  tends  to  zero,  the  bound  was  not 
informative  and  indeed  one  can  try  with  p  tends  to  zero  or  infinity,  the  optimal 
deterministic  solution  was  an  arbitrarily  bad  approximation  to  the  optimal  a  priori 
solution. 

Jaillet  (1988)  presented  an  example  showing  that  an  optimal  TSP  tour  was  not 
necessarily  a  good  solution  to  the  corresponding  PTSP  problem.  In  this  example.  Figure 
III. 6  contains  24  white  nodes  (and  no  black  nodes)  that  are  positioned  at  the  vertices  of 
two  concentric  12-gons  as  shown  in  Figure  III. 6a.  We  assume  P  =  Pj  (Bemouilli  process 
with  parameter  p).  In  Figure  III. 6b,  two  tours  have  been  designed  through  this  set  of 
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nodes:  tour  a  is  an  optimal  TSP  tour  and  tour  6  is  an  alternative  tour  (see  Jaillet  1 985  for 
the  numerical  derivations).  One  could  then  show  that,  for  a  probability  of  presence  of  0.5, 
the  expected  length  of  tour  a  is  31%  greater  than  the  expected  length  of  tour  h.  This 
numerical  example  raised  the  following  question:  in  general,  how  well  would  a  TSP  tour 
do  as  a  solution  to  a  PTSP  problem?  Actually,  PTSP  introduced  many  features  that  were 
different  from  those  of  its  famous  special  case,  the  TSP.  The  TSP  is  a  special  case  of  the 
PTSP  in  which  all  nodes  are  black.  It  is  then  natural  to  investigate  the  possible  links 
between  the  two  problems.  At  the  end  the  author  came  up  with  a  comment  that  in  fact, 
one  could  show  (using  a  generalization  of  the  star-shaped  example  in  Figure  III.  6)  that 
the  TSP  could  indeed  be  arbifrarily  bad  under  the  condition  of  Figure  III. 6  and  his 
Theorem  2. 

Based  on  the  previous  discussion,  the  deterministic  method  results  are  not 
necessarily  the  same  as  the  probabilistic  method  results  for  the  same  problem. 
Consequently,  there  is  a  conflict  between  the  economic  analysis  as  a  deterministic 
method  and  Markov  decision  process  (MDP)  as  a  probabilistic  method.  But  the  question 
that  rises  now  is:  which  method  should  we  go  with?  In  fact,  the  economic  analysis  (the 
deterministic  method)  is  simple,  straightforward  and  easy  to  be  applied  by  INDOT 
personnel.  Therefore,  the  economic  analysis  results  are  considered  in  our  study. 
Consequently,  the  conclusion  here  is  that  making  spot  repairs  ever>'  1 5  years  and/or  at 
condition  rating  of  7  is  the  best  scenario  that  INDOT  can  apply  for  the  3 -coat  system  in 
the  future. 
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Figure  III.4:  PV  ($/ft2)  Comparison  of  3  COAT  System/ZINC 
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Figure  III.5:  EUAC  ($/ft2)  for  3-COAT  System/  ZINC 
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Figure  III. 2:  PV  Comparison  of  3-COAT  System  /Lead  for  Different  Scenarios 
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Figure  III.3:  EUAC  ($  /ft2)  for  3-COAT  System  /  Lead  for  Different  Scenarios 
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Figure  III.6:  Graph  and  tours  for  the  numerical  example, 
(a)  A  24-node  graph,  (b)  Two  tours  of  the 
24-node  graph. 
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Figure  III.7:  PV  ($/ft2)  Comparison  of  3-COAT  System  /Lead  Including  Disruption  Cost 
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Figure  III.8:  EUAC  ($/ft2)  for  3-COAT  System  /  Lead  Including  Disruption  Cost 
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Figure  III.9:  PV  ($/ft2)  of  3-COAT  System/ZINC  Including  Disruption  Cost 
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Figure  111.10:  EUAC  ($/ft2)  3-COAT  System/  ZINC  Including  Disruption  Cost 
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Table  III.6:     Paint  Age  Sensitivity  Analysis 
3-Coat/Zinc 
Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

Age  (n) 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost 

0 

1 

$2.80 

$20.49 

Present  Value  for  Fv  after  15  yrs 

12 

1.511069 

$1.00 

$12.54 

Present  Value  for  Fv  after  30  yrs 

24 

2.283328 

$1.00 

$8.41 

Present  Value  for  Fv  after  30  yrs 

36 

3.450266 
5.213589 

$1.00 
$1.00 

$5.64 
$3.79 

Present  Value  for  Fv  after  45  yrs 

48 

Total  Present  Value  TPV/ft^ 
EUAC 

$50.86 
$3.62 

State  4  Spot 

$51.54 
$3.67 

Note:  Sin"  =  10.76  $/m^ 

The  values  of  PV  and  EUAC  of  state  3  including  the  disruption  cost  are  close  to  that  of  state  4, 

therefore,  this  scenario  is  the  best  one  if  the  bridge  reaches  state  3  or  condition  rating  7  after  12  years 


Table  ill. 8:     Initial  cost  Sensitivity  Analysis 
3-Coat/Zinc 
Spot  Repairs  at  State  3 


Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

Age  (n) 

Infl. Factor 

Initi.Cost 

Total(disrup.) 

Present  Value  for  initial  Cost 

0 

1 

$8.85 

$26.53 

Present  Value  for  Fv  after  15  yrs 

15 

1.675349 

$3.16 

$12.66 

Present  Value  for  Fv  after  30  yrs 

30 

2.806794 

$3.16 

$7.69 

Present  Value  for  Fv  after  45  yrs 

45 

4.702359 

$3.16 

$4.67 

Total  Present  Value  TPV/ft^ 
EUAC 

$51.54 
$3.67 

State  4  Spot 

$51.54 
$3.67 

Note:$/ff  =  10.76  $/m^ 


The  values  of  PV  and  EUAC  of  state  3  including  the  disruption  cost  are  the  same  as  that  of  state  4, 
therefore,  this  scenario  is  the  best  one  if  the  bridge  is  painted  by  3-Coat  system  where 
the  initial  cost  is  less  than  $8.85/ft2  and  the  rehabilitation  cost  every  15  years  is  S3.16/ft2. 


Table  ill. 7:     Paint  Age  Sensitivity  Analysis 
3-Coat/Lead 
Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

Age  (n) 

Infl. Factor 

Init.Cost 

Total(disrup.) 

Present  Value  for  initial  Cost 

0 

1 

$4.00 

$21.69 

Present  Value  for  Fv  after  12  yrs 

12 

1.511069 

$1.50 

$12.87 

Present  Value  for  Fv  after  24  yrs 

24 

2.283328 

$1.50 

$8.64 

Present  Value  for  Fv  after  36  yrs 

36 

3.450266 

$1.50 

$5.79 

Present  Value  for  Fv  after  48  yrs 

48 

5.213589 

$1.50 

$3.89 

Total  Present  Value  TPV/ft^ 
EUAC 

$52.88 
$3.77 

State  4  Spot 

$53.13 
$3.78 

Note:  $1  fir  =  10.76  $/m^ 

The  values  of  PV  and  EUAC  of  state  3  including  the  disruption  cost  are  close  to  that  of  state  4, 

therefore,  this  scenario  is  the  best  one  if  the  bridge  reaches  state  3  or  condition  rating  7  after  12  years  old. 


Table  111.9:     Initial  cost  Sensitivity  Analysis 
3-Coat/Lead 
Spot  Repairs  at  State  3 
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Inflation  Factor 


Item 


Present  Value  for  initial  Cost 


Present  Value  for  Fv  after  15  yrs 


Present  Value  for  Fv  after  30  yrs 


Present  Value  for  Fv  after  45  yrs 


0.035 


Age  (n)    Infl. Factor  Init.Cost    Total(disrup.) 


0 


15 


30 


45 


1.675349 


2.806794 


4.702359 


$9.80 


$3.68 


$3.68 


$3.68 


$27.49 


$12.97 


$7.88 


$4.78 


Total  Present  Value  TPVIfV 
EUAC 


State  4  Spot 


$53.12 
$3.78 


$53.13 
$3.78 


Note.  $lff=  10.76  $/m^ 

The  values  of  PV  and  EUAC  of  state  3  including  the  disruption  cost  are  close  to  that  of  state  4, 

therefore,  this  scenario  is  the  best  one  if  the  bridge  is  painted  by  3-Coat  system  where 

the  initial  cost  is  less  than  $9.8/ft^  and  the  rehabilitation  cost  every  15  years  is  $3.68/ft2. 
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III.2  INDOT  MAINTENANCE  PLAN 

A  maintenance  plan  is  a  very  important  objective.  A  life  cycle  cost  analysis  was 
done  for  different  paint  systems  to  determine  the  best  scenario  of  rehabilitation  for  the 
maintenance  plan.  A  maintenance  plan  is  a  set  of  rehabilitation  procedure  steps  that 
satisfy  the  minimum  rehabilitation  cost. 

III.2.1  Proposed  Maintenance  Plan  Procedure 

As  well  as  the  implementation  described  in  the  last  section  and  based  on  the  life 
cycle  cost  analysis  for  determination  of  the  best  rehabilitation  scenario,  a  proposed 
maintenance  plan  and  its  procedural  steps  is  summarized  in  Figure  III.l  1. 

The  steel  bridge  paint  should  be  inspected,  as  usual,  every  two  years  by  INDOT 
inspectors.  During  the  inspection,  the  paint  condition  will  be  rated.  If  the  rate  of  paint 
condition  reaches  7  or  below,  spot  painting  must  be  done  on  the  bridge,  regardless  of 
whether  or  not  the  paint  is  new.  In  case  of  new  paint,  the  question  to  the  inspector  is,  Is 
the  paint  life  over  15  years?  If  the  answer  is  Yes,  the  spot  painting  wdll  be  automatically 
performed,  even  if  the  rating  of  the  paint  condition  is  away  above  7.  If  the  bridge  has 
been  painted  within  1 5  years,  the  bridge  still  goes  through  routine  biannually  inspection. 
In  other  words,  nothing  is  done  to  the  paint  job  until  the  rating  reaches  7  or  the  paint  life 
is  over  15  years. 

When  an  old  paint  condition  rating  is  above  7,  the  inspector  should  ask,  "Was  the 
bridge  painted  in  the  last  15  years?"  If  the  bridge  has  not  been  painted  in  the  last  15 
years,  the  spot  painting  should  be  automatically  done.  Based  on  the  recommendation  of 
life  cycle  cost  analysis,  spot  painting  every  15  years  is  the  most  economic  policy,  despite 
the  condition  rating  not  reaching  7. 
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In  case  a  bridge  was  painted  in  the  last  1 5  years,  the  inspector  should  examine  the 
factor  of  bridge  life.  When  the  bridge  life  is  reaching  60  years,  there  is  no  point  in 
repainting  the  bridge,  since  the  bridge  will  be  reconstructed  or  demolished. 

After  being  spot-painted,  any  bridge  can  be  categorized  as  an  "old  bridge". 
Thereafter,  the  maintenance  procedures  should  be  repeated  until  the  bridge  is 
reconstructed. 
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Figure  111.11:  Maintenance  Plan  Procedure  Flowchart 
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Appendix  A: 
Regression  Output 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F 

Value 

Pr  >  F 

Model 

1 

2.5150380 

2.5150380 

7.01 

0.0100 

Error 

71 

25.4849620 

0.3589431 

Corrected  Total 

72 

28.0000000 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.089823 

9.985310 

0.5991 

6.0000 

Source 

DF 

Type  I  SS 

Mean  Square 

F 

Value 

Pr  >  F 

AGE 

1 

2.5150380 

2.5150380 

7.01 

0.0100 

Source 

DF 

Type  II  SS 

Mean  Square 

F 

Value 

Pr  >  F 

AGE 

1 

2.5150380 

2.5150380 

7.01 

0.0100 

T  for  HO 

Pr  >  |T| 

Std 

Error  of 

Parameter 

Estimate 

Parameter 

=  0 

Estimate 

INTERCEPT        7 

.177012484 

15. 

94 

0.0001 

0. 

,45014850 

AGE             -0 

.051327307 

-2. 

65 

0.0100 

0. 

,01939050 
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12:40  Sunday,  May  24,  199i 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Source 

Model 

Error 

Corrected  Total 


DF 
3 

69 

72 
R-Square 
0.127369 


Sum  of 
Squares 

3.5663205 

24.4336795 

28.0000000 

C.V. 

9.917875 


Mean 
Square   F  Value 


1.1887735 
0.3541113 

Root  MSE 
0.5951 


3.36 


Pr  >  F 
0.0237 


RATE  Mean 
6.0000 
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Source 

AGE 

AGE2 

AGE3 

Source 

AGE 

AGE2 

AGE3 


DF 

1 
1 
1 

DF 

1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


2.5150380 
0.0092365 
1.0420460 


2.5150380 
0.0092365 
1.0420460 


7.10 
0.03 
2  .94 


Type  II  SS   Mean  Square   F  Value 


1.1818535 
1.0511788 
1.0420460 


1.1818535 
1.0511788 
1.0420460 


3.34 
2.97 
2.94 


Pr  >  F 

0.0096 
0.8722 
0.0908 

Pr  >  F 

0.0720 
0.0894 
0.0908 


Parameter 

INTERCEPT 
AGE 
AGE2 
AGE  3 


Estimate 

11.21397561 

-0.79338953 

0.03972385 

-0.00065516 


T  for  HO: 
Parameter=0 

4.55 
-1.83 

1.72 
-1.72 


Pr  > 


0.0001 
0.0720 
0.0894 
0.0908 


Std  Error  of 
Estimate 

2.46417003 
0.43428471 
0.02305595 
0.00038192 


Paint  1 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

2 

2.8372088 

1.4186044 

3.95     0.0238 

Error 

70 

25.1627912 

0.3594684 

Corrected  Total 

72 

28.0000000 

R-Sguare 

C.V. 

Root  MSE 

RATE  Mean 

0.101329 

9.992615 

0.5996 

6.0000 

94 


Source 

AGE 
Zl 


DF 

1 
1 


Type  I  SS   Mean  Square   F  Value 


2.5150380 
0.3221708 


2.5150380 
0.3221708 


7.00 
0.90 


Pr  >  F 

0.0101 
0.3470 


Source 


DF 


Type  II  SS   Mean  Square   F  Value 


Pr  >  F 


AGE 

1 

2.7213358 

2, 

.7213358 

7 

.57     0.0075 

Zl 

1 

0.3221708 
T  for  HO: 

0, 

.3221708 
Pr  >  1 T 1 

0 

.90     0.3470 
Std  Error  of 

Parameter 

Estimate 

Parameter=0 

Estimate 

INTERCEPT 

7, 

.404864946 

14.50 

0.0001 

0.51074197 

AGE 

-0, 

.053920039 

-2.75 

0.0075 

0.01959700 

Zl 

-0, 

.195126915 

-0.95 

0.3470 

0.20611257 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Source 

Model 

Error 

Corrected  Total 


DF 

7 

65 

72 

R-Square 

0.214443 


Sum  of 
Squares 

6.0044073 

21.9955927 

28.0000000 

C.V. 

9.695270 


Mean 
Square   F  Value 


0.8577725 
0.3383937 

Root  MSE 
0.5817 


2.53 


Pr  >  F 
0.0227 


RATE  Mean 
6.0000 
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Source 

AGE 

AGE2 

AGE3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 

Source 

AGE 

AGE2 

AGE  3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 

DF 

1 
1 
1 
1 
1 
1 
1 

DF 

1 
1 
1 
1 
1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


2.5150380 
0.0092365 
1.0420460 
0.3058123 
0.0001544 
2.0521323 
0.0799878 


2.5150380 
0.0092365 
1.0420460 
0.3058123 
0.0001544 
2.0521323 
0.0799878 


0.2455279 
0.1849023 
0.1368142 
0.2317617 
0.1719455 
0.1205558 
0.0799878 


0.2455279 
0.1849023 
0.1368142 
0.2317617 
0.1719455 
0.1205558 
0.0799878 


7.43 
0.03 
3  .08 
0.90 
0.00 
6.06 
0.24 


Type  II  SS   Mean  Square   F  Value 


0.73 
0.55 
0.40 
0.68 
0.51 
0.36 
0.24 


Pr  >  F 

0.0082 
0.8693 
0.0840 
0.3453 
0.9830 
0.0165 
0.6285 

Pr  >  F 

0.3975 
0.4624 
0.5271 
0.4109 
0.4785 
0.5527 
0.6285 


Parameter 

INTERCEPT 

AGE 

AGE2 

AGE3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 

Estimate 

38.81026093 
-4.59619975 

0.20218756 

-0.00283768 

-28.11616802 

3 .86121075 
-0.16396185 

0.00217928 


T  for  HO: 
Parameter=0 

1.15 
-0.85 

0.74 
-0.64 
-0.83 

0.71 
-0.60 

0.49 


Pr  > 


0.2561 
0.3975 
0.4624 
0.5271 
0.4109 
0.4785 
0.5527 
0.6285 


Std  Error  of 
Estimate 

33.87278278 
5.39584736 
0.27352339 
0.00446282 

33.97398020 
5.41675567 
0.27470082 
0.00448241 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

4 

3.8721328 

0.9680332 

2.73     0.0362 

Error 

68 

24.1278572 

0.3548216 

Corrected  Total 

72 

28.0000000 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.138290 

9.927817 

0.5957 

6.0000 

96 


Source 

AGE 
AGE2 
AGE3 
Zl 

Source 

AGE 
AGE2 
AGE3 
Zl 


DF 

1 
1 
1 
1 

DF 

1 
1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


2.5150380 
0.0092365 
1.0420460 
0.3058123 


2.5150380 
0.0092365 
1.0420460 
0.3058123 


7.09 
0.03 
2.94 
0.86 


Type  II  SS   Mean  Square   F  Value 


1235090 
0170542 
0342484 


0.3058123 


1.1235090 
1.0170542 
1.0342484 
0.3058123 


3.17 
2.87 
2.91 
0.86 


Pr  >  F 

0.0097 
0.8723 
0.0911 
0.3565 

Pr  >  F 


0796 
0950 
0923 


0.3565 


Parameter 


Estimate 


T  for  HO : 
Parameter=0 


Pr  > 


Std  Error  of 
Estimate 


INTERCEPT 

AGE 

AGE2 

AGE  3 

Zl 


11.25654570 
-0.77441275 
0.03909080 
-0.00065272 
-0.19462403 


4.56 
-1.78 

1.69 
-1.71 
-0.93 


0.0001 
0.0796 
0.0950 
0.0923 
0.3565 


2.46706630 
0.43520034 
0.02308913 
0.00038231 
0.20963993 


T1 
Polynomial  Regression  Paint  (1): 
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Y  =  8.79682  +  6.30E-03X  -  7.81E-03X**2  +  1.16E-04X**3 
R-Sq  =  0.778 

Analysis  of  Variance 

SOURCE  DF  SS         MS          F          P 

Regression  3  29.6770    9.89234    82.0007         0 

Error  70  8.4446    0.12064 

Total  ■•  73  38.1216 


SOURCE  DF  Seq  SS  F  P 

Linear  1  29.1387  233.554  0 

Quadratic  1  0.4013  3.31996  7.27E-02 

Cubic  1  0.1370  1.13573  0.290218 


Macro  is  running  . . .  please  wait 

Regression         ^ 


The  regression  equation  is 
y  =  9.49  -  0.143  x 


Predictor 

Coef 

StDev 

T 

P 

Constant 

9.4888 

0.2154 

44.04 

0.000 

X 

-0.142790 

0.009343 

-15.28 

0.000 

S  =  0.3532       R-Sq  =  76.4%      R-Sq(adj)  =  76.1% 
Analysis  of  Variance  .  :  , 


Source 

DF 

SS 

MS 

F 

P 

Regression 

1 

29.139 

29.139 

233.55 

0.000 

Error 

72 

8.983 

0.125 

Total 

73 

38.122 

Macro  is  running  . . .  please  wait 

Polynomial  Regression 


Y  =  9.06329  -  8.21E-02X  -  1.78E-03X**2 
R-Sq  =  0.775 


Analysis  of  Variance 


SOURCE 

DF 

SS 

MS 

F 

P 

Regression 

2 

29, 

.5400 

14.7700 

122, 

.200 

0 

Error 

71 

8, 

.5816 

0.1209 

Total 

73 

38, 

.1216 

SOURCE 

DF 

Seq 

SS 

F 

P 

Linear 

1 

29.1387 

233  . 

554 

0 

Quadratic 

1 

0.4013 

3.31996   7.27E- 

-02 

Regression  Analysis-lnter-Paint(1 ): 
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The  regression  equation  is 
Paint  Rate  =    9.A$    -    0.143  Age 


Predictor 

Coef 

StDev 

Constant 

9.4888 

0.2154 

Age 

-0.142790 

0.009343 

S  =  0.3532 

R-Sa  - 

=  76.4% 

m 

J. 

? 

44.04 

0.000 

15.28 

0.000 

R-Sq{adj)  =  76.1% 


Analysis  of  Variance 


Source 

DF 

SS 

MS 

77 

•p 

Regress; 

ion     1 

29.139 

29. 

139 

233. 

55 

0.000 

Error 

72 

8.983 

0. 

125 

Total 

73 

38.1 

.22 

Unusual 

Observat 

ions 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Res; 

Ldual 

St  Resid 

14 

26.0 

5.0000 

5. 

.7763 

0.0517 

-0, 

.7763 

-2.22R 

53 

23.0 

7.0000 

6. 

.2047 

0.0412 

0 

.7953 

2.27R 

54 

10.0 

8.0000 

8. 

,0609 

0.1250 

-0 

.0609 

-0.18  X 

58 

23.0 

7.0000 

6, 

.2047 

0.0412 

0 

.7953 

2.27R 

65 

7.0 

8.0000 

8. 

.4893 

0.1517 

-0 

.4893 

-1.53  X 

74 

1.0 

9.0000 

9. 

.3460 

0.2063 

-0 

.3460 

-1.21  X 

R  denotes  an  observation  winh  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  1.50 
9  rows  with  no  reolicates 


0.1415   DF{pure  e: 


Regression  Analysis 


Z.^    F^=  -G-^  ;p^-<^'-- 


The  regression  equation  is 

Paint  Rate  =  9.06  -  0.0821  Age  -  0.00178  Agesq 


Predictor  Coef 
Constant  9.0  63  3 
Age  -0.08210 

Agesq      -0.0  017  33  2 

S  =  0.3477       R-Sq  = 

Analysis  of  Variance 


StDev  T        P 

0.3154  28.73  0.000 

0.03455  -2.38  0,-020^ 

0.0009787  -1.82  01^73,/ 

"  ■ 5 %      R-Sq(adj)  =  76.9% 


Source 

DF 

SS 

MS 

Regression 

2 

29.540 

.  14.770 

Error 

71 

3.582 

0.121 

Total 

73 

38.122 

Source 

DF 

Seq  SS 

Age 

1 

29.139 

Agesq 

1 

0.401 

122, 


P 
0.000 


Unusual  Observations 
Obs        Age    Paint  Ra 
14       26.0      5.0000 


Fit   StDev  Fit   Residual     St  Resid 
5.7231      0.0587     -0.7231       -2.11R 


42 

32.0 

5.0000 

4.6099 

0.1948 

0.3901 

1.35  X 

53 

23.0 

7.0000 

6.2315 

0.0432 

0.7685 

2.23R 

54 

10.0 

8.0000 

8.0639 

0.1230 

-0.0639 

-0.20  X 

58 

23.0 

7.0000 

6.2315 

0.0432 

0.7685 

2.23R 

65 

7.0 

8.0000 

8.4012 

0.1570 

-0.4012 

-1.29  X 

74 

1.0 

9.0000 

8.9794 

0.2858 

0.0206 

0.10  X 
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R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  1.34 
9  rows  with  no  replicates 


0.2172   DF(pure  error) 


58 


Regression  Analysis 


The  regression  equation  is 
Paint  Rate  =  8.80  +  0.0063  Age 


0.00781  Agesq  +0.000116  Agecubic 


Predictor 

Constant 

Age 

Agesq 

Agecubic 

S  =  0.3473 


Coef 

8.7968 

0.00630 

-0.007809 

0.0001158 


StDev 

0.4023 

0.08985 

0.005738 

0.0001086 


21.87 
0.07 

-1.36 
1.07 


0.000 
0.944 
0.178 
0.290 


R-Sq  =  77.8% 


R-Sq(adj)  =  76.9% 


Analysis  of  Variance 


Source 

DF 

ss 

MS 

F 

P 

Regression 

3 

29.6770 

9.8923 

82. 

00 

0.000 

Error 

70 

8.4446 

0.1206 

Total 

73 

38.1216 

Source 

DF 

Seq  SS 

Age 

1 

29.1387 

Agesq 
Agecubi 

c 

1 
1 

0.4013 
0.1370 

Unusual 

Observations 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

3 

13, 

.0 

8.0000 

7 

.8133 

0.1518 

0, 

.1867 

0.60  X 

14 

26, 

.0 

5.0000 

5 

.7163 

0.0590 

-0, 

.7163 

-2.09R 

42 

32, 

.0 

5.0000 

4, 

.7952 

0.2609 

0, 

.2048 

0.89  X 

53 

23, 

.0 

7.0000 

6 

.2191 

0.0447 

0, 

.7809 

2.27R 

54 

10, 

.0 

8.0000 

8 

.1945 

0.1736 

-0, 

.1946 

-0.65  X 

58 

23, 

.0 

7.0000 

6 

.2191 

0.0447 

0 

.7809 

2.27R 

65 

7 

.0 

8.0000 

8 

.4980 

0.1812 

-0, 

.4980 

-1.68  X 

73 

13 

.0 

8.0000 

7 

.8133 

0.1518 

0, 

.1867 

0.60  X 

74 

1 

.0 

9.0000 

8, 

.7954 

0.3337 

0, 

.2046 

2.12RX 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  1.35   P  =  0.2158   DF(pure  error)  =  5S 
9  rows  with  no  replicates 
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Regression  Plot 

Y  =  9.06329  -  8.21  E-02X  -  1 .78E-03X"2 
R-Sq  =  0.775 


10 


20 


30 


Age 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  Paint  Ra) 


-1  0  1 

Normal  Score 


102 


Regression  Plot 

Y  =  9.48882 -0.142790X 
R-Sq  =  0.764 


9    — ' 


CO 

DC 


CO 
CL 


6    — 


5    — 


Age 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  Paint  Ra) 


-1  0  1 

Normal  Score 
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Regression  Plot 

Y  =  8.79682  +  6.30E-03X  -  7.81E-03X**2  +  1.16E-04X"3 
R-Sq  =  0.778 


03 

c 
CL 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  Paint  Ra) 


■1  0 

Normal  Score 


Paint  2 
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12:40  Sunday,  May  24,  1958 

I2-.        ,„, 

1 

General 

Linear  Models  Procedure 

1  C\  1  r-^.  •sTV'  ~i — -^ 

ependent 

Variabl' 

e:  RATE 

ource 

DF 

Sum  of 
Squares 

Mean 
Square 

F  Value     Pr  >  F 

odel 

1 

173.24946 

173.24946 

1401.23     0.0001 

rror 

313 

38.69975 

0.12364 

Drrected 

Total 

314 

211.94921 

R- Square 

C.V. 

Root  MSE 

RATE  Mean 

1 

0.817410 

4.388369 

0.3516 

8.0127 

Durce 

DF 

Type  I  SS   Mean  Square 

F  Value     Pr  >  F 

3E 

1 

173.24946 

173.24946 

1401.23     0.0001 

Durce 

DF 

Type  II  SS   Mean  Square 

F  Value     Pr  >  F 

3E 

1 

173.24946 
T  for  HO: 

173.24946 
Pr  >  |t| 

1401.23     0.0001 
Std  Error  of 

irameter 

Estimate 

Parameter=0 

Estimate 

^TERCEPT 

9 

.154824493 

251.65 

0.0001 

0.03637923 

5E 

-0 

.124487791 

-37.43 

0.0001 

0.00332562 

Paint  2 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

3 

176.08035 

58.69345 

508.90     0.0001 

Error 

311 

35.86886 

0.11533 

Corrected  Total 

314 

211.94921 

R-Sguare 

C.V. 

Root  MSE 

RATE  Mean 

0.830767 

4.238379 

0.3396 

8.0127 
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Source 

AGE 
AGE  2 
AGE  3 

Source 

AGE 
AGE2 
AGE  3 


DF 

1 
1 
1 

DF 

1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


173.24946 
0.51146 
2.31942 


173.24946 
0.51146 
2.31942 


1502.15 

4.43 

20.11 


Type  II  SS   Mean  Square   F  Value 


7.4068635 
1.8075775 
2.3194247 


7.4068635 
1.8075775 
2.3194247 


64.22 
15.67 
20.11 


Pr  >  F 

0.0001 
0.0360 
0.0001 

Pr  >  F 

0.0001 
0.0001 
0.0001 


Parameter 

INTERCEPT 
AGE 
AGE2 
AGE3 


Estimate 

9.258946262 
■0.200983652 

0.010336822 
■0.000348433 


T  for  HO: 
Parameter=0 

140.49 

-8.01 

3.96 

-4.48 


Pr  > 


0.0001 
0.0001 
0.0001 
0.0001 


Std  Error  of 
Estimate 

0.06590502 
0.02507969 
0.00261106 
0.00007770 


Paint  2  194 

12:40  Sunday,  May  24,  1998 

,  General  Linear  Models  Procedure 

\  108 

Dependent  Variable:  RATE 

Sum  of  Mean 

Source  DF        Squares        Square   F  Value     Pr  >  F 

Model                     4  176.74850  44.18715    389.14     0.0001 

Error                   310  35.20061  0.11355 

Corrected  Total         314  211.94921 

R-Square  C.V.  Root  MSE           RATE  Mean 

0.833920  4.205479  0.3370              8.0127 

Source  DF      Type  I  SS   Mean  Square   F  Value     Pr  >  F 

AGE  1  173.24946  173.24946  1525.75  0.0001 

AGE2  1  0.51146  0.51146  4.50  0.0346 

AGE3  1  2.31942  2.31942  20.43  0.0001 

Zl  1  0.66825  0.66825  5.89  0.0158 

Source  DF     Type  II  SS   Mean  Square   F  Value     Pr  >  F 

AGE  1  7.2336428  7.2336428  63.70  0.0001 

AGE2  1  1.6279323  1.6279323  14.34  0.0002 

AGE3  1  2.0429357  2.0429357  17.99  0.0001 

Zl  1  0.6682507  0.6682507  5.89  0.0158 

T  for  HO:      Pr  >  | T |     Std  Error  of 
Parameter  Estimate      Parameter=0  Estimate 

INTERCEPT  9.236278094  139.82  0.0001  0.06605767 

AGE  -0.198754909  -7.98  0.0001  0.02490196 

AGE2  0.009840288  3.79  0.0002  0.00259887 

AGE3  -0.000328803  -4.24  0.0001  0.00007752 

Zl  0.110556949  2.43  0.0158  0.04557743 


Paint  2 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Source 

DF 

Squares 

Model 

7 

179.22781 

Error 

307 

32.72140 

Corrected  Total 

314 

211.94921 

R-Square 

C.V. 

0.845617 

4.074440 

Mean 
Square   F  Value 


25.60397 
0.10658 

Root  MSE 
0.3265 


240.22 


Pr  >  F 
0.0001 


RATE  Mean 
8.0127 


Source 

AGE 

AGE2 

AGE  3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 

Source 

AGE 

AGE2 

AGE3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 

DF 

1 

1 

1 

1- 

1 

1 

1 

DF 

1 
1 
1 
1 
1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


173 .24946 
0.51146 
2.31942 
0.66825 
1.50311 
0.08818 
0.88792 


173 .24946 
0.51146 
2.31942 
0.66825 
1.50311 
0.08818 
0.88792 


4.0387631 
0.3851595 
0.5603813 
0.1662791 
0.6915095 
0.9631843 
0.8879200 


4.0387631 
0.3851595 
0.5603813 
0.1662791 
0.6915095 
0.9631843 
0.8879200 


1625.47 
4.80 

21.76 
6.27 

14.10 
0.83 
8.33 


Type  II  SS   Mean  Square   F  Value 


37.89 
3.61 
5.26 
1.56 
6.49 
9.04 
8.33 


Pr  >  F 

0.0001 
0.0292 
0.0001 
0.0128 
0.0002 
0.3638 
0.0042 

Pr  >  F 

0.0001 
0.0582 
0.0225 
0.2126 
0.0113 
0.0029 
0.0042 


Parameter 

INTERCEPT 

AGE 

AGE2 

AGE3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 


Estimate 

9.220953652 
-0.166557032 

0.005338528 
-0.000189415 

0.208567849 
-0.182641478 

0.026217927 
-0.000879608 


T  for  HO: 
Parameter=0 

130.29 
-6.16 

1.90 
-2.29 

1.25 
-2.55 

3.01 
-2.89 


Pr  > 


0.0001 
0.0001 
0.0582 
0.0225 
0.2126 
0.0113 
0.0029 
0.0042 


Std  Error  of 
Estimate 

0.07076992 
0.02705737 
0.00280833 
0.00008261 
0.16698416 
0.07170461 
0.00872148 
0.00030475 


Regression  Analysis(l2): 
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The  regression  equation  is 
Paint  Rate  =  9.>^  -  0.124  Age 
0 

Predictor  Coef  StDev 
Constant  9.15482  0.03638 
Age         -0.124488     0.003326 


T 
251.65 
-37.43 


P 
0.000 
0.000 


S  =  0.3516       R-Sq  = 
Analysis  of  Variance 


il.7^ 


R-Sq(adj)  =  81.7% 


Source 

DF 

SS 

MS 

F         P 

Regression     1 

173. 

25 

173 

.25 

1401. 

23     0.000 

Error 

313 

38. 

70 

0 

.12 

Total 

314 

211. 

95 

Unusual 

Observations 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

6 

15.0 

8.0000 

7 

.2875 

0.0277 

0.7125 

2.03R 

16 

15.0 

8.0000 

7 

.2875 

0.0277 

0.7125 

2.03R 

29 

16.0 

8.0000 

7 

.1630 

0.0301 

0.8370 

2.39R 

30 

16.0 

8.0000 

7 

.1630 

0.0301 

0.8370 

2.39R 

33 

16.0 

8.0000 

7 

.1630 

0.0301 

0.8370 

2.39R 

34 

16.0 

8.0000 

7 

.1630 

0.0301 

0.8370 

2.39R 

38 

15.0 

8.0000 

7 

.2875 

0.0277 

0.7125 

2.03R 

40 

16.0 

8.0000 

7 

.1630 

0.0301 

0.8370 

2.39R 

75 

15.0 

8.0000 

7 

.2875 

0.0277 

0.7125 

2.03R 

149 

18.0 

6.0000 

6 

.9140 

0.0354 

-0.9140 

-2.61R 

162 

19.0 

6.0000 

6 

.7896 

0.0382 

-0.7896 

-2.26R 

185 

19.0 

6.0000 

6, 

.7896 

0.0382 

-0.7896 

-2.26R 

192 

17.0 

6.0000 

7, 

.0385 

0.0327 

-1.0385 

-2.97R 

202 

18.0 

6.0000 

6, 

.9140 

0.0354 

-0.9140 

-2.61R 

225 

18.0 

6.0000 

6, 

.9140 

0.0354 

-0.9140 

-2.61R 

233 

16.0 

6.0000 

7 

.1630 

0.0301 

-1.1630 

-3.32R 

239 

18.0 

6.0000 

6, 

.9140 

0.0354 

-0.9140 

-2.61R 

244 

18.0 

6.0000 

6, 

.9140 

0.0354 

-0.9140 

-2.61R 

247 

18.0 

6.0000 

6, 

.9140 

0.0354 

-0.9140 

-2. SIR 

250 

24.0 

5.0000 

6, 

.1671 

0.0531 

-1.1671 

-3.36KX 

261 

18.0 

6.0000 

6, 

.9140 

0.0354 

-0.9140 

-2.61R 

273 

24.0 

6.0000 

6, 

.1671 

0.0531 

-0.1671 

-0.48  X 

281 

23.0 

6.0000 

6, 

.2916 

0.0501 

-0.2916 

-0.84  X 

287 

25.0 

6.0000 

6, 

.0426 

0.0562 

-0.0426 

-0.12  X 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  14.91 
2  rows  with  no  replicates 


0.0000   DF(pure  error)  =  289 


A 


egression  Analysis 

^1   F^- 5.sq58,p^  o,ol58 


The  regression  equation  is 

Paint  Rate  =  9^6  -  0.201  Age  +  0.0103  Age2  -0.000348  Age3 


Predictor 


Coef 


StDev 


Constant 
Age 
Age2 
Age3 

S  =  0.3396 


9.25895 

-0.20098 

0.010337 

-0.00034843 

R-Sq 


0.06591 

0.02508 

0.002611 

0.00007770 

=  83  .1% 


140.49 

-8.01 

3  .96 

-4.48 


0.000 
0.000 
0.000 
0.000 


111 


R-Sq(adj)  =  82.9% 


Analysis  of  Variance 


Source 

DF 

SS 

MS 

F         P 

Regression 

3 

176.080 

58. 

693 

508. 

90     0.000 

Error 

311 

35.869 

0. 

115 

Total 

314 

211.949 

Source 

DF 

Seq 

SS 

Age 

1 

173  .249 

Age  2 

1 

0.511 

Age3 

1 

2.319 

Unusual 

Observations 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

29 

16.0 

8.0000 

7, 

.2623 

0.0359 

0.7377 

2.18R 

30 

16.0 

8.0000 

7, 

.2623 

0.0359 

0.7377 

2.18R 

33 

16.0 

8.0000 

7, 

.2623 

0.0359 

0.7377 

2.18R 

34 

16.0 

8.0000 

7, 

.2623 

0.0359 

0.7377 

2.18R 

40 

16.0 

8.0000 

7, 

.2623 

0.0359 

0.7377 

2.18R 

149 

18.0 

6.0000 

6, 

.9583 

0.0391 

-0.9583 

-2.84R 

162 

19.0 

6.0000 

6, 

.7819 

0.0427 

-0.7819 

-2.32R 

185 

19.0 

6.0000 

6, 

.7819 

0.0427 

-0.7819 

-2.32R 

192 

17.0 

6.0000 

7, 

.1177 

0.0372 

-1.1177 

-3.31R 

202 

18.0 

6.0000 

6. 

.9583 

0.0391 

-0.9583 

-2.84R 

225 

18.0 

6.0000 

6, 

.9583 

0.0391 

-0.9583 

-2.84R 

233 

16.0 

6.0000 

7, 

.2623 

0.0359 

-1.2623 

-3.74R  ■ 

239 

18.0 

6.0000 

6, 

.9583 

0.0391 

-0.9583 

-2.84R 

244 

18.0 

6.0000 

6, 

.9583 

0.0391 

-0.9583 

-2.84R 

247 

18.0 

6.0000 

6, 

.9583 

0.0391 

-0.9583 

-2.84R 

250 

24.0 

5.0000 

5. 

.5726 

0.1342 

-0.5726 

-1.84  X 

258 

22.0 

6.0000 

6, 

.1302 

0.0796 

-0.1302 

-0.39  X 

261 

18.0 

6.0000 

6, 

.9583 

0.0391 

-0.9583 

-2.84R 

273 

24.0 

6.0000 

5 

.5726 

0.1342 

0.4274 

1.37  X 

281 

23.0 

6.0000 

5, 

.8651 

0.1037 

0.1349 

0.42  X 

283 

22.0 

6.0000 

6, 

.1302 

0.0796 

-0.1302 

-0.39  X 

287 

25.0 

6.0000 

5, 

.2506 

0.1714 

0.7494 

2.56RX 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  14.12 
2  rows  with  no  replicates 


P  =  0.0000   DF(pure  error)  =  289 
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Regression  Plot 

Y  =  9.25895  -  0.200984X  +  1.03E-02X**2  -  3.48E-04X"3 
R-Sq  =  0.831 


10 


15 


20 


25 


Age 


Paint  1  (State) 

14:15  Thurs 

2 

;day.  May  28,  1998 

s\ 

General 

Linear  Models  Procedure 

st^ 

Dependent 
Source 

Variable:  RATE 

DF 

Sum  of 
Squares 

Mean 
Square 

F  Value     Pr  >  F 

pa\^(V) 

Model 

1 

66.125592 

66.125592 

332.80     0.0001 

Error 

58 

11.524408 

0.198697 

Corrected 

Total         59 

77.650000 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.851585 

6.507359 

0.4458 

6.8500 

Source 

DF 

Type  I  SS   Mean  Square 

F  Value     Pr  >  F 

AGE 

1 

66.125592 

66.125592 

332.80     0.0001 

Source 

DF 

Type  II  SS   Mean  Square 

F  Value     Pr  >  F 

AGE 

1 

66.125592 

66.125592 

332.80     0.0001 

Parameter 

Estimate 

T  for  HO: 
Parameter=0 

Pr  >  |t| 

Std  Error  of 
Estimate 

INTERCEPT 

AGE 

9.500528787 
-0.142373972 

60.79 
-18.24 

0.0001 
0.0001 

0.15627383  - 

0.00780443 

Paint  1  (State)  4 

14:15  Thursday,  May  28,  1998 

General  Linear  Models  Procedure 

dependent  Variable:  RATE 


Source 

DF 

Sum  of 
Squares 

Mean 
Square 

F  Value 

Pr  >  F 

lodel 

2 

66.199891 

33.099946 

164.78 

0.0001 

:rror 

57 

11.450109 

0.200879 

'orrected 

Total 

59 

77.650000 

R-Square 

c.v. 

Root  MSE 

RATE  Mean 

0.852542 

6.542998 

0.4482 

6.8500 

ource 

DF 

Type  I  SS   Mean  Square 

F  Value 

Pr  >  F 

GE 
1 

1 
1 

66.125592 
0.074300 

66.125592 
0.074300 

329.18 
0.37 

0.0001 
0.5455 

ource 

DF 

Type  II  SS  Mean  Square 

F  Value 

Pr  >  F 

GE 

1 

1 

1 

65.335605 
0.074300 

65.335605 
0.074300 

325.25 
0.37 

0.0001 
0.5455 

arameter 

Estimate 

T  for  HO: 
Parameter=0 

Pr  >  |t| 

Std  Error  of 
Estimate 

NTERCEPT 
3E 

1 

9 

-0 

0 

490484836 
143085589 
077639613 

60.07 

-18.03 

0.61 

0.0001 
0.0001 
0.5455 

0. 
0. 
0. 

15799522 
00793393 
12766070 
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Paint  1  (State) 
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14:15  Thursday,  May  28,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Source 

Model 

Error 

Corrected  Total 


DF 

3 

56 

59 

R-Square 

0.898076 


Slim   of 
Squares 

69.735570 
7.914430 

77.650000 

C.V. 

5.488141 


Mean 
Square   F  Value 


23.245190 
0.141329 

Root  MSE 
0.3759 


164.48 


Pr  >  F 
0.0001 


RATE  Mean 
6.8500 
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Source 

AGE 

AGE2 

AGE3 


DF 

1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


66.125592 
3.466438 
0.143541 


66.125592 
3.466438 
0.143541 


467.88 

24.53 

1.02 


Pr  >  F 

0.0001 
0.0001 
0.3179 


Source 

AGE 

AGE2 

AGE3 


DF 

1 
1 

1 


Type  II  SS   Mean  Square   F  Value 


0.1108375 
0.4070476 
0.1435409 


0.1108375 
0.4070476 
0.1435409 


0.78 
2.88 
1.02 


Pr  >  F 

0.3796 
0.0952 
0.3179 


Parameter 


Estimate 


T  for  HO: 
Parameter =0 


Pr  >  |t| 


Std  Error  of 
Estimate 


INTERCEPT 
AGE 
AGE2 
AGE3 


8.904326281 

0.082658870 

-0.012551592 

0.000161970 


40.71 
0.89 

-1.70 
1.01 


0.0001 
0.3796 
0.0952 
0.3179 


0.21871900 
0.09333875 
0.00739592 
0.00016072 


Paint  1  (State) 
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14:15  Thursday,  May  28,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Source 

Model 

Error 

Corrected  Total 


DP 

4 

55 

59 

R-Square 

0.905410 


Suia   of 
Squares 

70.305109 
7.344891 

77.650000 

C.V. 

5.334832 


Mean 
Square   F  Value 


17  .576277 
0.133543 

Root  MSE 
0.3654 


131.61 


Pr  >  F 
0.0001 


RATE  Mean 
6.8500 
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Source 

AGE 
AGE2 
AGE3 
Zl 


DF 

1 
1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


66.125592 
3.466438 
0.143541 
0.569539 


66.125592 
3.466438 
0.143541 
0.569539 


495.16 

25.96 

1.07 

4.26 


Pr  >  F 

0.0001 
0.0001 
0.3044 
0.0436 


Source 

AGE 
AGE2 
AGES 
Zl 


DF 

1 
1 
1 
1 


Type  II  SS   Mean  Square   F  Value 


0.2433591 

0.6117683 
0.2606442 
0.5695387 


0.2433591 
0.6117683 
0.2606442 
0.5695387 


1.82 
4.58 
1.95 
4.26 


Pr  >  F 

0.1826 
0.0368 
0.1680 
0.0436 


Parameter 


Estimate 


T  for  HO: 
Parameter=0 


Pr 


Std  Error  of 
Estimate 


INTERCEPT 

AGE 

AGE2 

AGE3 

Zl 


8.780668720 
0.125666302 
-0.015737932 
0.000222005 
0.223489259 


39.75 
1.35 

-2.14 
1.40 
2.07 


0.0001 
0.1826 
0.0368 
0.1680 
0.0436 


0.22088028 
0.09309071 
0.00735302 
0.00015891 
0.10821971 


Paint  1  (State) 
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14:15  Thursday,  May  28,  199E 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Slim  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

7 

70.464732 

10.066390 

72.85     0.0001 

Error 

52 

7.185268 

0.138178 

Corrected  Total 

59 

77.650000 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.907466 

5.426618 

0.3717 

6.8500 

Source 

AGE 

AGE2 

AGE  3 

Zl 

AGE*! 

Zl 

AGE2 

*Z1 

AGE3 

*Z1 

Source 

AGE 

AGE2 

AGES 

Zl 

AGE*; 

Zl 

AGE2 

*Z1 

AGE3 

*Z1 

DF 

1 
1 

1 
1 
1 
1 

1 

DF 

1 
1 
1 
,  1 
1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


66.125592 
3.466438 
0.143541 
0.569539 
0.043709 
0.001177 
0.114737 


66.125592 
3.466438 
0.143541 
0.569539 
0.043709 
0.001177 
0.114737 


0.2678078 
0.6213582 
0.3119026 
0.0523017 
0.0903061 
0.1095097 
0.1147372 


0.2678078 
0.6213582 
0.3119026 
0.0523017 
0.0903061 
0.1095097 
0.1147372 


478.55 
25.09 
1.04 
4.12 
0.32 
0.01 
0.83 


Type  II  SS   Mean  Square   F  Value 


1.94 
4.50 
2.26 
0.38 
0.65 
0.79 
0.83 


Pr  >  F 

0.0001 
0.0001 
0.3128 
0.0475 
0.5762 
0.9268 
0.3664 

Pr  >  F 

0.1698 
0.0387 
0.1390 
0.5411 
0.4225 
0.3774 
0.3664 


Parameter 

INTERCEPT 

AGE 

AGE2 

AGE3 

Zl 

AGE*Z1 

AGE2*Z1 

AGE3*Z1 


Estimate 

8.765813443 
0.153673721 

-0.018255036 
0.000278222 
0.278062638 

-0.187303091 
0.016775893 

-0.000374774 


T  for  HO: 
Parameter=0 

32.17 
1.39 

-2.12 
1.50 
0.62 

-0.81 
0.89 

-0.91 


Pr  > 


0.0001 
0.1698 
0.0387 
0.1390 
0.5411 
0.4225 
0.3774 
0.3664 


Std  Error  of 
Estimate 

0.27244959 
0.11038443 
0.00860858 
0.00018518 
0.45196492 
0.23168948 
0.01884426 
0.00041128 


Macro  is  running  . . .  please  wait 

Polynomial  Regression(SI): 
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Y   =    8.90433 
R-Sq   =    0.89i 


.27E-02X   -    1.26E-02X**2    +    1.62E-04X**3 


Analysis  of  Variance 


SOURCE 

DF 

SS 

MS 

F 

P 

Regression 

3    69 

7356 

23 

.2452 

164 

476 

0 

Error 

56     7 

9144 

0 

.1413 

Total 

59    77 

6500 

SOURCE 

DF 

Seq  SS 

F 

P 

Linear 

1 

66.1256 

332 

797 

0 

Quadratic 

1 

3.4564 

24.5207 

6.90E- 

-06 

Cubic 

1 

0.1435 

1.01565 

0.317887 

Macro  is  running  . . .  please  wait 


Polynomial  Regression 


Y  =  9.05508  -  7.04E-03X  -  5.17E-03X**2 
R-Sq  =  0.896 

Analysis  of  Variance 


SOURCE 

DF 

SS 

MS          F 

P 

Regression 

2    69 

5920 

34. 

7960    246.138 

0 

Error 

57     8 

0580 

0. 

1414 

Total 

59    77 

6500 

SOURCE 

DF 

Seq  SS 

F 

p   \  y^ 

^ 

Linear 

1 

66.1256 

332 

797, 

-iL.  VX 

Quadratic 

1 

3.4664 

24.5207; 

6.90E-06 

^^ 


Regression  Analysis 

The  regression  equation  is 
Paint  Rate  =  9 .#0  -  0.142  Age 


Predictor 

Constant 

Age 

Coef 

9.5005 

-0.142374 

0 

StDev           T 
0.1563       60.79 
007804      -18.24 

P 

0.000 
0.000 

S  =  0.4458 

R-Sq  = 

=  85 

2%      R-Sq(adj)  = 

84.9% 

Analysis  of  Variance 


Source 

DF 

ss 

Regression     1 

66.126 

Error 

58 

11.524 

Total 

59 

77.650 

Unusual 

Observat 

.ions 

Obs 

Age 

Paint  Ra 

16 

0.0 

9.0000 

33 

1.0 

9.0000 

35 

1.0 

9.0000 

45 

2.0 

9.0000 

46 

2.0 

9.0000 

57  . 

25.0 

5.0000 

MS 

F 

P 

66.126 

332.80 

0.000 

0.199 
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Fit   StDev  Fit 


9.5005 
9.3582 
9.3582 
9.2158 
9.2158 
5.9412 


0.1563 
0.1490 
0.1490 
0.1419 
0.1419 
0.0751 


Residual 
-0.5005 
-0.3582 
-0.3582 
-0.2158 
-0.2158 
-0.9412 


St  Resid 
-1.20  X 
-0.85  X 
-0.85  X 
-0.51  X 
-0.51 


X 


-2.14R 


R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  2.30 
4  rows  with  no  replicates 


0.0181   DF(pure  error)  =  44 


Regression  Analysis 


The  regression  equation  is 
Paint  Rate  =  9.06  -  0.0070  Age 


0.00517  Age2 


Predictor 

Coef 

StDev 

T 

P 

Constant 

9.0551 

0.1596 

56.74 

0.000 

Age 

-0.00704 

0.02811 

-0.25 

0.803 

Age  2 

-0.005173 

0.001045 

-4.95 

'D'.OOO 

S  =  0.3760 

R-Sq  = 

=  89.6% 

R- 

-Sq{ 

adj)  = 

89.3% 

Analysis  of  Variance 


Source 

DF 

SS 

MS 

F 

p 

Regression 

2 

69.592 

34. 

796 

246. 

14 

0.000 

Error 

57 

8.058 

0. 

141 

Total 

59 

77.650 

,,.  . 

Source 

DF 

Seq  SS 

Age 

1 

66.126 

Age2 

1 

3  .466 

Unusual 

Observations 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

16 

0. 

.0 

9. 

0000 

9 

.0551 

0. 

.1596 

-0, 

.0551 

-0.16  X 

54 

20. 

.0 

6. 

,0000 

6 

.8452 

0. 

,0628 

-0. 

.8452 

-2.28R 

55 

20. 

.0 

6. 

0000 

6 

.8452 

0. 

,0628 

-0. 

.8452 

-2.28R 

56 

29, 

.0 

5. 

0000 

4 

.5006 

0. 

.1949 

0, 

.4994 

1.55  X 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  0.71   P  =  0.7404   DF(pure  error) 
4  rows  with  no  replicates 


44 


Regression  Analysis 


The  regression  equation  is 
Paint  Rate  =  8.90  +  0.0827  Age 


0.0126  Age2  +0.000162  Age3 


Predictor 

Coef 

StDev 

T 

P 

Constant 

8.9043 

0.2187 

40 

71 

0 

000 

Age 

0.08266 

0.09334 

0 

89 

0 

380 

Age2 

-0.012552 

0.007396 

-1 

70 

0 

095 

Age3 

0.0001620 

0.0001607 

1 

01 

0 

318 
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S  =  0.3759       R-Sq  =  89.8%      R-Sq(adj)  =  89.3^ 
Analysis  of  Variance 


Source 

DF 

SS 

MS 

F         P 

Regression 

3 

69.736 

23 

245     164. 

48     0.000 

Error 

56 

7.914 

0 

141 

Total 

59 

77.650 

Source 

DF 

Seq  SS 

Age 
Age2 

1 

1 

66.126 
3.466 

Age3 

1 

0.144 

Unusual 

Observations 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

16 

0 

0 

9 

0000 

8 

9043 

0.2187 

0.0957 

0.31  X 

54 

20 

0 

6 

0000 

6 

8326 

0.0640 

-0.8326 

-2.25R 

55 

20 

0 

6 

0000 

6 

8326 

0.0640 

-0.8326 

-2.25R 

56 

29 

0 

5 

0000 

4 

6958 

0.2747 

0.3042 

1.19  X 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 

Pure  error  test  -  F  =  0.69   P  =  0.7486   DF(pure  error)  =  44 
4  rows  with  no  replicates 


Regression  Plot 

Y  =  9.50053- 0.1 42374X 
R-Sq  =  0.852 
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CC 

DC 

c 

'co  7 
Q- 


Age 


Regression  Plot 

Y  =  9.05508  -  7.04E-03X  -  5.17E-03X**2 
R-Sq  =  0.896 
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CO 

DC 

•4—' 

c 

'cO 
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Regression  Plot 

Y  =  8.90433  +  8.27E-02X  -  1 .26E-02X"2  +  1 .62E-04X**3 
R-Sq  =  0.898 


10 


20 


30 


Age 


Paint  2  238 

12:40  Sunday,  May  24,  1998 


General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

1 

254.52948 

254.52948 

1598.24     0.0001 

Error 

275 

43.79543 

0.15926 

Corrected 

Total 

276 

298.32491 

R- 

-Square 

C.V. 

Root  MSE 

RATE  Mean 

0 

853196 

5.087075 

0.3991 

7.8448 

Source  DF  Type  I  SS   Mean  Square   F  Value     Pr  >  F 

AGE  1  254.52948     254.52948   1598.24     0.0001 

Source  DF  Type  II  SS   Mean  Square   F  Value     Pr  >  F 

AGE  1  254.52948     254.52948   1598.24     0.0001 

T  for  HO:  Pr  >  | T |     Std  Error  of 

Parameter  Estimate  Parameter=0                   Estimate 

INTERCEPT  9.276541532  215.23       0.0001       0.04309968 

AGE  -0.156019671  -39.98       0.0001       0.00390264 
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Paint  2 
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12:40  Sunday,  May  24,  199S 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Stun  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

.2  . 

254.68150 

127.34075 

799.46     0.0001 

Error 

274 

43.64341 

0.15928 

Corrected  Total 

275 

298.32491 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.853705 

5.087496 

0.3991 

7.8448 

125 


Source 

AGE 
AGE  2 

Source 


DF 

1 
1 

DF 


Type  I  SS   Mean  Square   F  Value 


254.52948 
0.15202 


254.52948   1597.98 
0.15202      0.95 


Type  II  SS   Mean  Square   F  Value 


Pr  >  F 

0.0001 
0.3295 

Pr  >  F 


AGE 
AGE2 

1 
1 

202.46711 
0.15202 

202 
0 

.46711   : 
.15202 

L271 
0 

.12     0.0001 
.95     0.3295 

Parameter 

Estimate 

T  for  HO: 
Parameter=0 

Pr  >  1 T 

Std  Error  of 
Estimate 

INTERCEPT 

AGE 
AGE2 

9. 
-0. 
-0. 

.268165327 
.154199101 
.000056974 

210.89 

-35.65 

-0.98 

0.0001 
0.0001 
0.3295 

0.04394771 
0.00432503 
0.00005832 

Paint    2 
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12:40    Sunday,    May   24,    1998 
General    Linear   Models    Procedure 
Dependent   Variable:    RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

yiodel 

3 

258.05711 

86.01904 

583.18     0.0001 

Srror 

273 

40.26780 

0.14750 

:orrected  Total 

276 

298.32491 

1 

R-Sguare 

C.V. 

Root  MSE 

PATE  Mean 

1 

0.865020 

4.895732 

0.3841 

7.8448 
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^GE2 
^GE3 


iGE 

^GE2 

^GE3 


DF 

1 
1 
1 

DF 

1 

1 
1 


Type  I  SS   Mean  Square   F  Value 


254.52948 
0.15202 
3.37561 


254.52948 
0.15202 
3.37561 


1725.61 

1.03 

22.89 


Type  II  SS   Mean  Square   F  Value 


2.8874888 
3.4450888 
3.3756105 


2.8874888 
3 .4450888 
3 .3756105 


19.58 
23.36 
22.89 


Pr  >  F 

0.0001 
0.3109 
0.0001 

Pr  >  F 

0.0001 
0.0001 
0.0001 


'arameter 

INTERCEPT 

iGE 

^GE2 

VGE3 


Estimate 

9.027775168 
-0.075278707 
-0.004885682 

0.000053545 


T  for  HO: 
Parameter=0 

137.45 

-4.42 

-4.83 

4.78 


Pr  > 


0.0001 
0.0001 
0.0001 
0.0001 


Std  Error  of 
Estimate 

0.06567823 
0.01701412 
0.00101093 
0.00001119 


Paint  2 
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12:40  Sunday,  May  24,  199? 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

2 

254.77530 

127.38765 

801.48     0.0001 

Error 

274 

43.54961 

0.15894 

Corrected  Total 

276 

298.32491 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.854020 

5.082026 

0.3987 

7.8448 
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Source 

AGE 
Zl 

Source 

AGE 
Zl 


DF 

1 

1 

DF 

1 
1 


Type  I  SS   Mean  Square   F  Value 


254.52948 
0.24582 


254.52948   1601.42 
0.24582      1.55 


Type  II  SS   Mean  Square   F  Value 


248.74084 
0.24582 


248.74084   1565.00 
0.24582      1.55 


Pr  >  F 

0.0001 
0.2147 

Pr  >  F 

0.0001 
0.2147 


Parameter 

INTERCEPT 

AGE 

Zl 


Estimate 

9.248280634 

-0.155418060 

0.062366077 


T  for  HO: 
Parameter=0 

189.96 

-39.56 

1.24 


Pr  > 


0.0001 
0.0001 
0.2147 


Std  Error  of 
Estimate 

0.04868574 
0.00392866 
0.05014843 
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12:40  Sunday,  May  24,  1958 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

yiodel 

4 

258.14737 

64.53684 

436.91     0.0001 

Srror 

272 

40.17754 

0.14771 

:orrected 

Total 

276 

298.32491 

R- 

-Square 

C.V. 

Root  MSE 

RATE  Mean 

0 

865323 

4.899224 

0.3843 

7.8448 
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Source 

\GE 

^GE2 

^GE3 

11 

Source 

\GE 
^GE2 
\GE3 
11 


DF 

1 

1 
1 
1 

DF 

1 
1 

1 
1 


Type    I    SS      Mean    Square      F   Value 


254.52948 
0.15202 
3.37561 
0.09026 


254.52948 
0.15202 
3 .37561 
0.09026 


2.9472396 
3 .2970008 
3.2322704 
0.0902590 


2.9472396 
3 .2970008 
3 .2322704 
0.0902590 


1723 .15 

1.03 

22.85 

0.61 


Type  II  SS  Mean  Square   F  Value 


19.95 

22.32 

21.88 

0.61 


Pr  >  F 

0.0001 
0.3112 
0.0001 
0.4351 

Pr  >  F 

0.0001 
0.0001 
0.0001 
0.4351 


Parameter 

[NTERCEPT 
^GE 
^GE2 
^GE3 

11 


Estimate 

9.014731777 

-0.076260043 

-0.004804587 

0.000052662 

0.038002863 


T  for  HO: 
Parameter=0 

132.94 

-4.47 

-4.72 

4.68 

0.78 


Pr  > 


0.0001 
0.0001 
0.0001 
0.0001 
0.4351 


Std  Error  of 
Estimate 

0.06781008 
0.01707248 
0.00101696 
0.00001126 
0.04861589 


Paint  2 
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12:40  Sunday,  May  24,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Source 

Model 

Error 

Corrected  Total 


DF 
6 

270 

276 

R-Square 

0.865414 


Sum  of 
Squares 

258.17460 
40.15031 

298.32491 

C.V. 

4.915669 


Mean 
Square   F  Value 


43.02910 
0.14870 

Root  MSE 
0.3856 


289.36 


Pr  >  F 
0.0001 


RATE  Mean 
7.8448 
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Source 

AGE 

AGE2 

AGE3 

Zl 

AGE2*Z1 

AGE3*Z1 

Source 

AGE 

AGE2 

AGES 

Zl 

AGE2*Z1 

AGE3*Z1 

DF 

1 

1 
1 
1 
1 

1 

DF 

1 
1 
1 
1 
1 
1 


Type  I  SS   Mean  Square   F  Value 


254.52948 
0.15202 
3 .37561 
0.09026 
0.00953 
0.01770 


254.52948 
0.15202 
3.37561 
0.09026 
0.00953 
0.01770 


2.0735807 
2.5513219 
2.5021819 
0.1013223 
0.0237829 
0.0176980 


2.0735807 
2.5513219 
2.5021819 
0.1013223 
0.0237829 
0.0176980 


1711.64 
1.02 
22.70 
0.61 
0.06 
0.12 


Type  II  SS   Mean  Square   F  Value 


13.94 

17.16 

16.83 

0.68 

0.16 

0.12 


Pr  >  F 

0.0001 
0.3129 
0.0001 
0.4366 
0.8004 
0.7304 

Pr  >  F 

0.0002 
0.0001 
0.0001 
0.4098 
0.6895 
0.7304 


Parameter 

INTERCEPT 

AGE 

AGE2 

AGE3 

Zl 

AGE2*Z1 

AGE3*Z1 


Estimate 

9.001146261 
-0.073517380 
-0.004921304 

0.000053927 

0.062124733 
-0.000570534 

0.000022192 


T  for  HO: 
Parameter=0 

118.72 

-3  .73 

-4.14 

4.10 

0.83 

-0.40 

0.34 


Pr  > 


0.0001 
0.0002 
0.0001 
0.0001 
0.4098 
0.6895 
0.7304 


Std  Error  of 
Estimate 

0.07581545 
0.01968757 
0.00118812 
0.00001315 
0.07526178 
0.00142663 
0.00006433 
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Regression  Plot 

Y  =  9.07656  -  0.101248X  -  2.03E-03X"2  -  2.87E-05X**3 
R-Sq  =  0.865 


10 


20 


Age 


30 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  Paint  Ra) 


3- 


2- 


■g 
"co 

<D 
DC 

■D 

CD         1 

N 

73 


0- 


T3 

C 

CO 


(f)     -1  - 


-2- 


-1  0  1 

Normal  Score 


Regression  Analysis(S2): 


/ 


Hr 


/\ 


i?  V- 


O'"^- 


A^\ 
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The  regression  equation  is 
Paint  Rate  =  9.03  -  0.0753  Age 


0.00489  Age2  +0.000054  Age3 


Predictor 
Constant 
Age 
Age2 
Age  3 


Coef 

9.02778 

-0.07528 

-0.004886 

0.00005354 


StDev 

0.06568 

0.01701 

0.001011 

0.00001119 


T 

137.45 

-4.42 

-4.83 

4.78 


P 

0.000 
0.000 
0.000 
0.000 


S  =  0.3841       R-Sq  =  86.5% 
Analysis  of  Variance 


R-Sq(adj)  =  86.4% 


Source 

DF 

ss 

MS 

F         P 

Regression 

3 

258. 

D57 

86. 

019 

583. 

18     0.000 

Error 

273 

40.. 

268 

0. 

148 

Total 

276 

298.. 

325 

Source 

DF 

Seq 

SS 

Age 

1 

254.. 

529 

Age2 

1 

0.. 

L52 

Age3 

1 

3.: 

376 

Unusual 

Observat 

ions 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

21 

16.0 

6.0000 

6 

.7919 

0.0350 

-0.7919 

-2.07R 

31 

23.0 

5.0000 

5 

.3633 

0.0891 

-0.3633 

-0.97  X 

52 

15.0 

8.0000 

6 

9800 

0.0333 

1.0200 

2.67R 

81 

16.0 

6.0000 

6 

7919 

0.0350 

-0.7919 

-2.07R 

86 

27.0 

5.0000 

4 

4875 

0.1396 

0.5125 

1.43  X 

106 

16.0 

6.0000 

6 

7919 

0.0350 

-0.7919 

-2.07R 

135 

14.0 

8.0000 

7 

1632 

0.0329 

0.8368 

2.19R 

137 

19.0 

7.0000 

6 

2010 

0.0502 

0.7990 

2.10R 

161 

26.0 

5.0000 

4 

9935 

0.3840 

0.0065 

1.52  X 

215 

14.0 

8.0000 

7 

1632 

0.0329 

0.8368 

2.19R 

218 

14.0 

8.0000 

7 

1632 

0.0329 

0.8368 

2.19R 

229 

20.0 

7.0000 

5 

9963 

0.0583 

1.0037 

2.64R 

231 

20.0 

7.0000 

5 

9963 

0.0583 

1.0037 

2.64R 

235 

16.0 

6.0000 

6 

7919 

0.0350 

-0.7919 

-2.07R 

237 

16.0 

6.0000 

6 

7919 

0.0350 

-0.7919 

-2.07R 

268 

10.0 

9.0000 

7 

8400 

0.0344 

1.1600 

3.03R 

270 

12.0 

9.0000 

7 

5134 

0.0338 

1.4866 

3.89R 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 


Pure  error  test  -  F  =  11.11 
3  rows  with  no  replicates 


P  =  0.0000   DF(pure  error)  =  252 


Regression  Analysis 


The  regression  equation  is 
Paint  Rate  =  9.2  8  -  0.156  Age 

Predictor       Coef       StDev  T        P 

Constant      9.27654      0.04310      215.23     0.000 
Age        -0.156020     0.003903      -39.98     0.000 
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S  =  0.3991       R-Sq  =  85.3%      R-Sq(adj)  =  85.3% 
Analysis  of  Variance 


Source 

DF 

SS 

MS 

F         P 

Regression    1 

254. 

,53 

254 

.53 

1598. 

24     0.000 

Error 

275 

43. 

.80 

0 

.16 

Total 

276 

298. 

,32 

Unusual 

Observations 

Obs 

Age 

Paint  Ra 

Fit 

StDev  Fit 

Residual 

St  Resid 

31 

23.0 

5.0000 

5, 

.6881 

0.0590 

-0.6881 

-1.74  X 

52 

15.0 

8.0000 

6, 

.9362 

0.0330 

1.0638 

2.67R 

86 

27.0 

5.0000 

5, 

.0640 

0.0736 

-0.0640 

-0.16  X 

135 

14.0 

8.0000 

7, 

.0923 

0.0305 

0.9077 

2.28R 

161 

26.0 

5.0000 

5, 

.2200 

0.0699 

-0.2200 

-0.56  X 

215 

14.0 

8.0000 

7, 

.0923 

0.0305 

0.9077 

2.28R 

218 

14.0 

8.0000 

7, 

.0923 

0.0305 

0.9077 

2.28R 

223 

22.0 

5.0000 

5, 

.8441 

0.0555 

-0.8441 

-2.14R 

229 

20.0 

7.0000 

5, 

.1561 

0.0486 

0.8439 

2.13R 

231 

20.0 

7.0000 

5, 

.1561 

0.0486 

0.8439 

2.13R 

268 

10.0 

9.0000 

7. 

.7163 

0.0242 

1.2837 

3.22R 

270 

12.0 

9.0000 

7. 

.4043 

0.0264 

1.5957 

4.01R 

R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 

Pure  error  test  -  F  =  11.99   P  =  0.0000   DF(pure  error)  =  252 
3  rows  with  no  replicates 


X 


y 


y 


Regression  Analysis(Paint  3): 


The  regression  equation  is 
Paint  Rate  =  8.88  -  0.123_Age 


p5 
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Predictor 

Coef 

StDev 

T 

P 

Constant 

8.8752 

0.1642 

54.06 

0.000 

Age 

-0.12252 

0.01098 

-11.16 

0.000 

S  =  0.3766       R-Sq  =  91.9% 
Analysis  of  Variance 


SS 

17.671 

1.560 

19.231 


Unusual  Observations 

Obs  Age    Paint  Ra 

8  33.0  5.000 

9  21.0  7.000 


R-Sq(adj)     =    91.2% 


Source 

DF 

Regression 

1 

Error 

11 

Total 

12 

MS 

F 

P 

17.671 

124.60 

0.000 

0.142 

Fit 
4.832 
6.302 


StDev  Fit 
0.258 
0.147 


Residual 
0.168 
0.698 


St  Resid 
0.61  X 
2.01R 


R  denotes  an  observation  with  a  large  standardized  residual 

X  denotes  an  observation  whose  X  value  gives  it  large  influence. 

Sum  of  squares  for  pure  error  is  (nearly)  zero. 
Cannot  do  pure  error  test. 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  Paint  Ra) 
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Normal  Score 
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Regression  Plot 

Y  =  8.87520 -0.122517X 
R-Sq  =  0.919 


Paint  3 


18 


14:15  Thursday,  May  28,  1998 
General  Linear  Models  Procedure 
Dependent  Variable:  RATE 


Sum  of 

Mean 

Source 

DF 

Squares 

Square 

F  Value     Pr  >  F 

Model 

2 

66.199891 

33.099946 

164.78     0.0001 

Error 

57 

11.450109 

0.200879 

Corrected  Total 

59 

77.650000 

R-Square 

C.V. 

Root  MSE 

RATE  Mean 

0.852542 

6.542998 

0.4482 

6.8500 
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Source 

AGE 
Zl 

Source 


DF 

1 
1 

DF 


Type  I  SS   Mean  Square   F  Value 


66.125592 
0.074300 


66.125592 
0.074300 


329.18 
0.37 


Type  II  SS   Mean  Square   F  Value 


Pr  >  F 

0.0001 
0.5455 

Pr  >  F 


AGE 
Zl 

1 
1 

65.335605 
0.074300 

65 
0 

335605 
074300 

325 
0 

.25     0.0001 
.37     0.5455 

Parameter 

Estimate 

T  for  HO: 
Parameter=0 

Pr  >  |T| 

Std  Error  of 
Estimate 

INTERCEPT 

AGE 

Zl 

9 

-0 

0 

490484836 
143085589 
077639613 

60.07 

-18.03 

0.61 

0.0001 
0.0001 
0.5455 

0.15799522 
0.00793393 
0.12766070 

Worksheet  size:  100000  cells 
Macro  is  running  . . .  please  wait 

Regression 


The  regression  equation  is 
V  =  8.08  -  0.116  'x 
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Predictor 
Constant 

X 

Coef 

8.08193 

-0.11616 

StDev 
0.08592 
0.01686 

T 
94.07 
-q_.89 

? 
0.000 
0.000 

S  =  0.5195 

R-Sq  = 

=  43.8% 

R-Sq(adj)  = 

42.8% 

Source 

DF 

SS 

MS 

F 

P 

Regression 

1 

12.809 

12.809 

47.47 

0.000 

Error 

51 

16.461 

0.270 

Total 

62 

29.270 

Macro  is  running  . . .  please  wait 

Polynomial  Regression 


Y  =  8.03525  -  1.32E-02X  -  1.13E-02X-"2 
R-Sq  =  0.461 

.-j?-alysis  of  Variance 

SOUR.CE  DF         SS         MS          F          P 

Regression  2  13.4399    6.74496    25.6464   8.92Z-09 

Error  60  15.7799    0.25300 

Total  62  29.2598 


SOUR.CE       D?     Seq  SS  F  P 

Linear        1    12. 8*0  9  3    47.4688   3  .58E-09 
Quadratic     1     0.6807    2.58810   0.112919   / 


Macro  is  running  . . .  please  wait 

Polynomial  Regression 


Y  =  8.04753  -  0.113598X  +  1.70E-02X**2  -  1.86E-03X'-3 
R-Sq  =  0.465 

Analysis  of  Variance 

SS        MS         F  P 

13.5979  4.53263    17.0640   4.28E-08 

15.6719  0.26563 
29.2598 


SOURCE 

DF 

Regression 

3 

Error 

59 

Total 

52 
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SOtlRCE 

DF 

Seq  SS 

F 

P 

Linear 

1 

12.8093 

47.4638 

3^_58Er09. 

Quadracic 

1 

0.6807 

2.58810 

0.112919 

Cubic 

1 

0.1080 

0.406518 

0.526209 

Regression  Analysis 


The  regression  equation  is  L^''^ 

Paint  Racing  =  8.08  -  0.116  3CoatAge 

Predictor       Coef       StDev  T        p 

Constanc      8.08193      0.08592       94.07     0.000 
3CoatAge     -0.11616      0.01686       -6.89     0.000 

S  =  0.5195       R-Sq  =  43.8%      R-Sq(adj)  =  42.8% 

Analysis  of  Variance  ■  '  -  , 

Source  DF  SS  MS          F         P 

Regression  1  12.809  12.809      47.47     0.000 

Error  61  16.461  0.270 

Total  62  29.270 

Unusual  Observations 

Obs   3CoatAge    Paint  Ra        Fit   StDev  Fit    Residual     St  Resid 

1  11.0      5.0000      6.8041      0.1454     -1.8041       -3.62R 

2  7.0      6.0000      7.2688      0.0904     -1.2688       -2.48R 
4       7.0      6.0000      7.2688      0.0904     -1.2688       -2.48R 

13        0.0      7.0000      8.0819      0.0859     -1.0819       -2.11R 

R  denotes  an  observation  with  a  large  standardized  residual 

Pure  error  test  -  F  =  8.31   P  =  0 .  00P_0^  DF(pure  error)  =  53 
2  rows  with  no  reolicates       '  - 
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CO 

DC 

'cO 

CL 
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Regression  Plot 

Y  =  8.08193- 0.11 61 63X 
R-Sq  =  0.438 


•  • 


10 


SCoatAge 
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Normal  Probability  Plot  of  the  Residuals 

(response  is  Paint  Ra) 
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Residuals  Versus  SCoatAge 

(response  is  Paint  Ra) 
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Residuals  Versus  the  Fitted  Values 

(response  is  Paint  Ra) 
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7.0  7.5  8.0 

Fitted  Value 
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Appendix  B: 
Non-Linear  Programming  (NLP) 

Calculations 
(Transition  Probability  Matrix 

Values) 
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Transition  Probabilities  for  Paint  Type  (1):  Interstate 


-:-*^^^ 

P(2) 

P(3) 

r   P(4) 

0.0  -  4.0 

0.943 

0.768 

0.701 

0.999 

5.0 '  8.0 

0.938 

0.759 

0.671 

0.496 

9.0-12.0 

0.887 

0.923 

0.899 

0.928 

13.0-16.0 

0.829 

0.802 

0.887 

0.961 

17.0-20.0 

0.891 

0.84 

0.816 

0.815 

21.0-24.0 

0.809 

0.748 

0.746 

0.794 

25.0-28.0 

0.612 

0.559 

0.59 

0.694 

29.0-32.0 

0.0001 

0.0001 

0.0001 

0.12 
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Transition  Probabilities  for  Paint  Type  (1):  State 


Age 

P(1) 

P{2} 

-  ^-^-. 

m 

0.0  -  4.0 

0.998 

0.479 

0.999 

0.476 

5.0  -  8.0 

0.922 

0.943 

0.001 

0.927 

9.0-12.0 

0.933 

0.862 

0.554 

0.692 

13.0-16.0 

0.917 

0.804 

0.647 

0.788 

17.0-20.0 

0.842 

0.724 

0.673 

0.818 

21.0-24.0 

0.683 

0.632 

0.633 

0.677 

25.0-28.0 

0.0001 

0.064 

0.062 

0.265 

29.0-32.0 

0.0001 

0.0001 

0.0001 

0.0001 
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Transition  Probabilities  for  Paint  Type  (2):  interstate 


Age 

lillMiii 

P{1) 

P(2) 

P(3) 

iiiiiii m 

P(4) 

0,0  -  4.0 

0.958 

0.0001 

0.0001 

0.0001 

5.0  -  8.0 

0.892 

0.909 

0.505 

0.648 

9.0-12.0 

0.831 

0.942 

0.999 

0.999 

13.0-16.0 

0.933 

0.865 

0.83 

0.571 

17.0-20.0 

0.885 

0.819 

0.749 

0.661 

21.0-24.0 

0.629 

0.684 

0.49 

0.466 

25.0-28.0 

0.0001 

0.0001 

0.0001 

0.0001 
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Transition  Probabilities  for  Paint  Type  (2):  State 


0.0  -  4.0 


0.935 


0.766 


0.301 


0.0001 


5.0  -  8.0 


0.923 


0.815 


0.553 


0.703 


9.0-12.0 


0.888 


0.797 


0.611 


0.879 


13.0-16.0 


0.844 


0.74 


0.674 


0.85 


17.0-20.0 


0.799 


0.617 


0.548 


0.824 


21.0-24.0 


0.0001 


0.0001 


0.999 


0.21 


25.0-28.0 


0.0001 


0.0001 


0.0001 


0.0001 
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Transition  Probabilities  for  Paint  Type  (3) 


Age 

P(1) 

P{2)V 

P(3) 

:  m 

P(4) 

0.0  -  4.0 

0.794 

0.996 

0.0001 

0.0001 

5.0  -  8.0 

0.748 

0.952 

0.633 

0.053 

9.0-12.0 

0.999 

0.851 

0.872 

0.744 

13.0-16.0 

0.923 

0.889 

0.851 

0.772 

17.0-20.0 

0.864 

0.863 

0.837 

0.797 

21.0-24.0 

0.841 

0.807 

0.784 

0.81 

25.0-28.0 

0.83 

0.694 

0.656 

0.782 

29.0-32.0 

0.219 

0.283 

0.336 

0.528 
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Transition  Probabilities  for  3-Coat  System  OZEU 


"iPi^ 

'    P(2) 

P{4)   1 

0.0  -  4.0 

0.0001 

0.856 

0.981 

0.33 

5.0  -  8.0 

0.944 

0.782 

0.829 

9.0-12.0 

0.897 

0.826 

0.866 

13.0-16.0 

0.867 

0.811 

0.863 

17.0-20.0 

0.808 

0.766 

0.834 

21.0-24.0 

0.632 

0.639 

0.758 

25.0-28.0 

0.0001 

0.129 

0.36 
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Appendix  C: 
Markov  Decision  Process  (MDP) 

Calculations 
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Application  of  Markov  Decision  Process  (MDP)  for  Interstate  Paint  type  (1) 

(0-4  years  Range) 


V1 
V2 
V3 
V4 
V5 

0.3284 
0.9686 
1.4326 
1.965 
220.3 

0.943 
0 
0 
0 

1 

0.057 
0.768 

0 

0 

0 

0 
0.232 
0.701 

0 

0 

0 

0 
0.299 
0.999 

0 

0 

0 

0 

0.001 

0 

State  1 
K 

C1k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

0.309681 

0.05521 

0 

0 

0 

0.364891 

State  2 
K 

C2k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
25 

0 
0 

0.743885 
0.9686 

0.332363 
0 

0 
0 

0 
0 

1 .076248 
25.9686 

"* 

State  3 
K 

C3k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 
25 

110 

0 
0 
0 

0 
0 
0 

1 .004253 
1.21771 
1.4326 

0.587535 

0.29475 

0 

0 
0 
0 

1.591788 
26.51246 
1 1 1 .4326 

^^m 

State  4 
K 

C4k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 

90 
180 

0 
0 
0 

0 
0 
0 

0 
0 
0 

1 .963035 
1.965 
1.965 

0.2203 
0 
0 

2.183335 
91.965 
181.965 

^^ 

State  5 
K 

C5k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

4 

220 

0.3284 

0 

0 

0 

0 

220.3284 

^^^K 

The  Result 

State  1 

Decision  K 

Policy  Y1 

1 

2 

3 

4 

5 
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Application  of  IVlarkov  Decision  Process  (IVIDP)  for  Interstate  Paint  type  (1) 

(5-8  years  Range) 


V1 
V2 
V3 
V4 
V5 


38.1933 
105.54 
155.8056 
208.4249 
254.3739 


0.938 
0 
0 
0 

1 


0.062 
0.759 

0 

0 

0 


0 
0.241 
0.671 

0 

0 


0 

0 
0.329 
0.496 

0 


0 
0 
0 
0.504 
0 


State  1 
K 


Clk 


P1*V1 


P2*V2 


P3*V3 


P4*V4 


P5*V5 


Value       Selection 


1 


35.82532 


6.54348 


0 


0 


0 


42.3688 


State  2 
K 


C2k 


P1*V1 


P2*V2 


P3*V3 


P4*V4 


P5*V5 


Value       Selection 


0 
25 


80.10486 
105.54 


37.54915 
0 


0 
0 


0 
0 


117.654 
130.54 


State  3 
K 


C3k 


P1*V1 


P2*V2 


P3*V3 


P4*V4 


P5*V5 


Value       Selection 


1 
2 
3 


0 

25 
110 


0 

0 
0 


104.5456 
132.4348 
155.8056 


68.57179 

31.26374 

0 


0 
0 
0 


173.1173 
188.6985 
265.8056 


State  4 
K 


C4k 


prvi 


P2*V2 


P3*V3 


P4*V4 


P5*V5 


Value       Selection 


231.5832  i»«§ 

298.4249 

388.4249 


1 
2 
3 


0 

90 
180 


0 
0 
0 


103.3788 
208.4249 
208.4249 


128.2044 
0 
0 


State  5 
K 


C5k 


prvi 


P2*V2 


P3*V3 


P4*V4 


P5*V5 


Value       Selection 


4 


220 


38.1933 


0 


0 


0 


0 


258.1933 


The  Result 

State  I 

Decision  kI 

Policy  Y1 
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Application  of  Markov  Decision  Process  (MDP)  for  Interstate  Paint  type  (1) 

(9-12  years  Range) 


V1 

8.8077 

0.887 

0.113 

0 

0 

0 

V2 

17.4681 

0 

0.923 

0.077 

0 

0 

V3 

42.6748 

0 

0 

0.899 

0.101 

0 

V4 

89.6218 

0 

0 

0 

0.928 

0.072 

V5 

227.9269 

1 

0 

0 

0 

0 

State  1 

K 

C1k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

7.81243 

1.973895 

0 

0 

0 

9.786325 

State  2 

K 

C2k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

0 

16.12306 

3.28596 

0 

0 

19.40902  IHIBHI 

2 

25 

0 

17.4681 

0 

0 

0 

42.4681 

State  3 

K 

C3k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

0 

0 

38.36465 

9.051802 

0 

47.41645  mgimgl 

2 

25 

0 

0 

36.27358 

13.44327 

0 

74.71685 

3 

110 

0 

0 

42.6748 

0 

0 

152.6748 

State  4 

K 

C4k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
90 

0 
0 

0 
0 

0 
0 

83.16903 
89.6218 

16.41074 
0 

99.57977 
179.6218 

::..:^'im 

3 

180 

0 

0 

0 

89.6218 

0 

269.6218 

State  5 

K 

C5k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

4 

220 

8.8077 

0 

0 

0 

0 

228.8077 

^"^^^^i^ 

The  Result 

State  1 

1 

2 

3 

4 

5 

Decision  K 

MB 

MHHi 

■HBH 

■HnMH 

Policy  Y1 

FV 

^™ 

^" 

^™ 
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Application  of  l\/larl<ov  Decision  Process  (IVIDP)  for  Interstate  Paint  type  (1) 

(13  - 16  years  Range) 


VI 
V2 
V3 
V4 
V5 

77.7264 
19.3459 
30.2022 
59.8996 
230.5538 

0.829 
0 
0 
0 

1 

0.171 
0.802 

0 

0 

0 

0 
0.198 
0.887 

0 

0 

0 

0 
0.113 
0.961 

0 

0 
0 
0 
0.039 
0 

State  1 
K 

C1k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

9.721186 

3.308149 

0 

0 

0 

13.02933 

State  2 
K 

C2k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
25 

0 
0 

15.51541 
19.3459 

5.980036 
0 

0 

0 

0 
0 

21.49545 
44.3459 

State  3 
K 

C3k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 

25 
110 

0 
0 
0 

0 
0 
0 

26.78935 
25.67187 
30.2022 

6.768655 

8.98494 

0 

0 
0 
0 

33.55801 
59.65681 
140.2022 

-^m 

State  4 
K 

C4k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5'V5 

Value 

Selection 

1 
2 
3 

0 

90 
180 

0 
0 
0 

0 
0 
0 

0 
0 
0 

57.56352 
59.8996 
59.8996 

8.991598 
0 
0 

66.55511 
149.8996 
239.8996 

1 

State  5 
K 

C5k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

4 

220 

1 1 ,7264 

0 

0 

0 

0 

231 .7264 

^mi 

The  Result 

State  1 
Decision  K 
Policy  Y1 

1 

2 

3 

4 

5 
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Application  of  Markov  Decision  Process  (MDP)  for  Interstate  Paint  type  (1) 

(17-20  years  Range) 


VI 
V2 
V3 
V4 
V5 

27.9045 

56.3495 

95.4811 

153.1388 

245.1141 

0.891 
0 
0 
0 

1 

0.109 
0.84 

0 

0 

0 

0 
0.16 
0.816 

0 

0 

0 

0 
0.184 
0.815 

0 

0 
0 
0 
0.185 
0 

State  1 
K 

C1k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

24.86291 

6.142096 

0 

0 

0 

31.00501 

State  2 
K 

C2k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
25 

0 
0 

47.33358 
56.3495 

15.27698 
0 

0 
0 

0 

0 

62.61056 
81.3495 

■^ 

State  3 
K 

C3k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 
25 

110 

0 
0 
0 

0 
0 

0 

77.91258 
81.15894 
95.4811 

28.17754 
22.97082 

0 

0 
0 
0 

106.0901 
129.1298 
205.4811 

^am 

State  4 
K 

C4k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 

90 
180 

0 
0 
0 

0 
0 
0 

0 
0 
0 

124.8081 
153.1388 
153.1388 

45.34611 
0 
0 

170.1542 
243.1388 
333.1388 

MM 

State  5 
K 

C5k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Seiect'O- 

4 

220 

27.9045 

0 

0 

0 

0 

247.9045 

.i^#m 

The  Result 

State  1 

Decision  K 

Policy  Y1 

1 

2 

3 

4 

5 

wm 
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Application  of  Marl<ov  Decision  Process  (MDP)  for  Interstate  Paint  type  (1) 

(21  -24  years  Range) 


V1 

53.0988 

0.809 

0.191 

0 

0 

0 

V2 

83.9882 

0 

0.748 

0.252 

0 

0 

V3 

121.02 

0 

0 

0.746 

0.254 

0 

V4 

173.9596 

0 

0 

0 

0.794 

0.206 

V5 

267.7889 

1 

0 

0 

0 

0 

State  1 

K 

Clk 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

42.95693 

16.04175 

0 

0 

0 

58.99868 

;  ,  ^,1JMJ 

State  2 

K 

C2k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
25 

0 
0 

62.82317 
83.9882 

30.49703 
0 

0 
0 

0 
0 

93.3202 

WBSBKM 

108.9882 

State  3 

K 

C3k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

0 

0 

90.28089 

44.18574 

0 

134.4666 

2 

25 

0 

0 

102.867 

26.09394 

0 

153.9609 

3 

110 

0 

0 

121.02 

0 

0 

231.02 

State  4 

K 

C4k 

P1*V1 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
90 

0 
0 

0 
0 

0 

0 

138.1239 
173.9596 

55.16451 
0 

193.2884 
263.9596 

B|i=44!aaSi: 

3 

180 

0 

0 

0 

173.9596 

0 

353.9596 

State  5 

K 

C5k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

4 

220 

53.0988 

0 

0 

0 

0 

273.0988 

iST^l^i 

The  Result 

State  1 

1 

2 

3 

4 

5 

Decision  K 
Policy  Y1 

■^F 

^1 

^B 

',t^tl[l^ 
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Application  of  Markov  Decision  Process  (MDP)  for  Interstate  Paint  type  (1) 

(25  -  28  years  Range) 


VI 
V2 
V3 
V4 
V5 

116.3918 
149.7228 
187.4459 
238.2442 
324.7526 

0.612 
0 
0 
0 

1 

0.388 
0.559 

0 

0 

0 

0 
0.441 
0.59 

0 

0 

0 

0 
0.41 
0.694 

0 

0 
0 
0 
0.306 
0 

State  1 
K 

C1k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

71.23178 

58.09245 

0 

0 

0 

129.3242 

■■" 

State  2 
K 

C2k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
25 

0 
0 

83.69505 
149.7228 

82.66364 
0 

0 
0 

0 
0 

166.3587 
174.7228 

g^'-'v^^H^^I 

State  3 
K 

C3k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 
25 

110 

0 
0 
0 

0 
0 
0 

110.5931 
159.329 
187,4459 

97.68012 

35.73663 

0 

0 
0 
0 

208.2732 
220.0656 
297.4459 

"^ 

State  4 
K 

C4k 

prvi 

P2'^V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 

90 
180 

0 

0 
0 

0 
0 
0 

0 
0 
0 

165.3415 
238.2442 
238.2442 

99.3743 
0 
0 

264.7158 
328.2442 
418.2442 

^.i^-~- 

State  5 
K 

C5k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

4 

220 

116.3918 

0 

0 

0 

0 

336.3918 

— 

The  Result 

State  1 

Decision  K 

Policy  Y1 

1 

2 

3 

4 

5 

apBKjijM^xjiiit^ 

f^?f'^~^^^K 
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Application  of  IVIarkov  Decision  Process  (IVIDP)  for  Interstate  Paint  type  (1) 

(29  -  32  years  Range) 


VI 
V2 
V3 
V4 
V5 

341.0915 
378.9943 
421.109 
467.9036 
526.9823 

0.0001 
0 
0 
0 

1 

0.9999 
0.0001 

0 

0 

0 

0 
0.9999 
0.0001 

0 

0 

0 

0 

0.9999 

0.12 

0 

0 
0 
0 
0.88 
0 

State  1 
K 

C1k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 

0 

0.034109 

378.9564 

0 

0 

0 

378.9905 

State  2 
K 

C2k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 

0 
25 

0 

0 

0.037899 
378.9943 

421.0669 
0 

0 
0 

0 
0 

421.1048 

403.9943 

^SicsSv.^-;"' 

State  3 
K 

C3k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 
25 

110 

0 
0 
0 

0 
0 
0 

0.042111 
357.9427 
421.109 

467.8568 

70.18554 

0 

0 
0 
0 

467.8989 

453.1282 

531.109 

^■^^^^ 

State  4 
K 

C4k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

1 
2 
3 

0 

90 
180 

0 
0 
0 

0 
0 
0 

0 
0 
0 

56.14843 
467.9036 
467.9036 

463.7444 
0 
0 

519.8929 
557.9036 
647.9036 

^^^SHH 

State  5 
K 

C5k 

prvi 

P2*V2 

P3*V3 

P4*V4 

P5*V5 

Value 

Selection 

4 

220 

341.0915 

0 

0 

0 

0 

561.0915 

'-.  \4''*f:; 

The  Result 

State  1 

Decision  K 

Policy  Y1 

1 

3 

4 

5 
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Appendix  D: 
Economic  Analysis  Calculations 


Interstate  paint  type  (1)  Scenario  #  1 


161 


n  =               1  = 

F  =                                                                /  *  L 

30              0.035 

''''n          A  -  AfA/F  i  A/1  +A(A/P  i  N)  -  F 

A(AJP.i,N) 

(1*i)''N-1 

Inflation  factor  = 

2.806794 

AI11=              $33.57 

includes  the  first  value. 

Interstate  paint  type  (1)  Scenario  #  2 


n  = 

1  = 

F  = 

10 

0.035 

25 

A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

P  = 

195 

Inflation  factor  = 

1.410599 

AI12  = 

$15.27  includes  thf 

3  first  value. 

Interstate  paint  type  (1)  Scenario  #  3 


n  = 

i  = 

F  = 

18 

0.035 

50 

A=A1(A/P,i 

,N1)  +A2(A/P,i 

,N2) 

P  = 

170 

Inflation  factor  = 

1.857489 

AI13  = 

$19.70 

includes  the  first  value. 

Interstate  paint  type  (1)  Scenario  #  4 


n  = 
18 
P  = 

0.035 
110 

F  = 
110 

A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

Inflation  factor  = 

1.857489 

AI14  = 

$23.68  includes  thf 

3  first  value. 

r 

Interstate  paint  type  (1)  Scenario  #  5 
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n=  i=  F= 

24  0.035  180        A=A1(A/P,i,N1)+A2(A/P,i,N2) 

P=  $310.62 

Inflation  factor  =  2.283328 

PV1  =  $78.83 

PV2=  $11.79  ,  .  ' 


AI15  = 


$28.43  includes  the  first  value. 


Interstate  paint  type  (2)  Scenario  #  1 


163 


n  =               1  = 

F=                                                                 i  *L 

26              0.035 

1 RO          A  -  AfA/F  i  N)  +A/^A/P  /  A/I  =  P 

—     + 

A(A/P,i,N) 

(1  +i)''N-1 

inflation  factor  = 

2.445959 

AI21  =              $26.07  includes  the  first  vaiue. 

Interstate  paint  type  (2)  Scenario  #  2 


n  = 
9 
P  = 

0.035 
160 

F  = 

20         A=A1(A/P,i 

,N1)  +A2(A/P,i 

,N2) 

Inflation  factor  = 

1.362897 

AI22  = 

$12.32  includes  the  first  value. 

Interstate  paint  type  (2)  Scenario  #  3 


n  = 
18 
P  = 

i  = 
0.035 
140 

F  = 

40        A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

Inflation  factor  = 

1.857489 

AI13  = 

$16.06  includes  the  first  value. 

Interstate  paint  type  (2)  Scenario  #  4 


n  = 

i  = 

F  = 

18 

0.035 

100 

A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

P  = 

80 

Inflation  factor  = 

1.857489 

A114  = 

$20.04 

includes  the  first  value. 

Interstate  paint  type  (2)  Scenario  #  5 
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n  = 

i  = 

F  = 

23 

0.035 

150        A 

=  A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

P  = 

$257.42 

Inflation  factor  = 

2.206114 

PV1  = 

$67.99 

- 

PV2  = 

$9.43 

A115  = 

$22.77 

includes  the  first  value. 

state  paint  type  (1)  Scenario  #  1 
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n  =               1  = 

F=                                                             /■  *L 

28              0.035 

"'''fl          A  -  AfA/F  i  N)  +A(A/P  i  N)  -  F 

—    + 

A(A/P,i,N) 

(1  *i)''N-1 

Inflation  factor  = 

imoMi 

AS11  =             $32.63  includes  the  first  value. 

State  paint  type  (1)  Scenario  #  2 


n  = 

i  = 

F  = 

14 

0.035 

25         A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

P  = 

195 

Inflation  factor  = 

1.618695 

AS12  = 

$16.36  includes  the  first  value. 

State  paint  type  (1)  Scenario  #  3 


n  = 
20 
P  = 

0.035 
170 

F  = 

50         A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

Inflation  factor  = 

1.989789 

AS13  = 

$20.56  includes  the  first  value. 
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State  paint  type  (1)  Scenario  #  4 


n  = 
20 
P  = 

0.035 
110 

F  = 
110 

A=A1(A/P,i 

,N1)  +A2(A/P,i 

,N2) 

inflation  factor  = 

1.989789 

AS14  = 

$24.17  includes  the  first  value. 

State  paint  type  (1)  Scenario  #  5 


n=  i=  F= 

24  0.035  180        A=A1(A/P,i,N1)+A2(A/P,i,N2) 

P=  $310.62 

Inflation  factor  =  2.283328 

PV1  =  $78.83 

PV2=  $11.79 


AS15  = 


$28.43  includes  the  first  value. 
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State  paint  type  (2)  Scenario  #  1 


n  =               i  =               F  =                                                              i  *L 
25  0.035  180        A  =  A(A/F,i,N) +A(A/P,i,N)  =  F     


(1  +i)''N-1 


Inflation  factor  =  2.363245 

AS21  =  $25.81  includes  the  first  value. 


A(A/P,i,N) 


State  paint  type  (2)  Scenario  #  2 


n  = 

i  = 

F  = 

9 

0.035 

20 

A=A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

P  = 

160 

Inflation  factor  = 

1.362897 

AS22  = 

$12.32 

includes  the  first  value. 

State  paint  type  (2)  Scenario  #  3 


n  = 

i  = 

F  = 

15 

0.035 

40 

A=A1(A/P,i 

,N1)  +A2(A/P,i 

,N2) 

P  = 

140 

Inflation  factor  = 

1.675349 

AS23  = 

$15.22  includes  the  first  value. 
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State  paint  type  (2)  Scenario  #  4 


n  = 
15 
P  = 

i  = 
0.035 
80 

F  = 

100        A 

=  A1(A/P,i 

,N1) 

+A2(A/P,i 

,N2) 

Inflation  factor  = 

1.675349 

AS24  = 

$19.92 

includes  the  first  value. 

State  paint  type  (2)  Scenario  #  5 


n  = 
20 
P  = 

0.035 
$30.00 

F  = 
150 

A=A1(A/P,i 

,N1)  +A2(A/P,i 

,N2) 

Inflation  factor  = 

1,989789 

AS25  = 

$23.39  includes  the  first  value. 

3-Coat/Lead 
Do  Nothing  Untill  State  5  and  then  Complete  Repaint 
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Age 
(yrs) 

Cost 
initial/ft^ 

Interest  Rate 

25 

$4.00 

7.00% 

Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost: 

0 

1 

$4.00 

Present  Value  for  Fv  after  25  yrs: 

25 

2.363245 

$1.74 

Present  Value  for  Fv  after  50  yrs: 

50 

5.5849269 

$0.76 

Total  Present  Value  TPV/ft^ 
EUAC 

$6.50 
$0.46 

3-Coat/Lead 

Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$4.00 

$4.00 

Present  Value  for  Fv  after  15  yrs: 

15 

1.675348831 

$1.50 

$0.91 

Present  Value  for  Fv  after  30  yrs: 

30 

2.806793705 

$1.50 

$0.55 

Present  Value  for  Fv  after  45  yrs: 

45 

4.702358551 

$1.50 

$0.34 

Total  Present  Value  TPV/ft^ 
EUAC 

$5.80 
0.413112 

3-Coat/Lead 

Spot  Repairs  at  State  4 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$4.00 

$4.00 

Present  Value  for  Fv  after  20  yrs: 

20 

1.989789 

$2.50 

$1.29 

Present  Value  for  Fv  after  40  yrs: 

40 

3.95926 

$2.50 

$0.66 

Total  Present  Value  TPV/ft^ 
EUAC 

$5.95 
0.423565 

3-Coat/Zinc 
Do  Nothing  Until!  State  5  and  then  Complete  Repaint 
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Age 

(yrs) 

Cost 
initial/ft^ 

Interst  Rate 

25 

$2.80 

7.00% 

Inflation  Factor 

0.035 

Item 

n 

Infl.Factor 

Cost 

Present  Value  for  initial  Cost: 

0 

1 

$2.80 

Present  Value  for  Fv  after  25  yrs: 

25 

2.363245 

$1.22 

Present  Value  for  Fv  after  50  yrs: 

50 

5.584927 

$0.53 

Total  Present  Value  TPV/ft' 
EUAC 

...,>. .  ....  .  . 

$4.55 
$0.32 

3-Coat/Zinc 

Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$2.80 

$2.80 

Present  Value  for  Fv  after  15  yrs: 

15 

1.675348831 

$1.00 

$0.61 

Present  Value  for  Fv  after  30  yrs: 

30 

2.806793705 

$1.00 

$0.37 

Present  Value  for  Fv  after  45  yrs: 

45 

4.702358551 

$1.00 

$0.22 

Total  Present  Value  TPV/ft^ 
EUAC 

$4.00 
0.284905 

3-Coat/Zinc 

Spot  Repairs  at  State  4 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$2.80 

$2.80 

Present  Value  for  Fv  after  20  yrs: 

20 

1.989789 

$2.00 

$1.03 

Present  Value  for  Fv  after  40  yrs: 

40 

3.95926 

$2.00 

$0.53 

Total  Present  Value  TPV/ft^ 
EUAC 

$4.36 
0.31036 
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Disruption  Cost  Calculations 


User  Cost 
Lane/day 

Road  of 
4-lanes 

Full  Paint 
4  days 

Spot  Repairs  for 
3  days 

Spot  Repairs  for 
2  days 

(1) 
$45,000 

$180,000 

$720,000 

$540,000 

$360,000 

Average  Bridge  Weight  =                                      177  tons     from  data 
Conversion  factor  =                                            115  ft^ /ton 

User  Cost 

Lane/day 

$/ton 

Road  of 

4-lanes 

$/ton 

Full  Paint 

4  days 

$/ton 

Spot  Repairs  for 

3  days 

$/ton 

Spot  Repairs  for 
2  days 

$/ton 

$254.24 

$1,016.95 

$4,067.80 

$3,050.85 

$2,033.90 

User  Cost 
Lane/day 

$/ft' 

Road  of 
4-lanes 

Full  Paint 
4  days 

$/ft' 

Spot  Repairs  for 
3  days 

S/ft' 

Spot  Repairs  for 
2  days 

$/ft^ 

$2.21 

$8.84 

$35.37 

$26.53 

$17.69 

(1)  Mr.  Mike  Long'^^',  INDOT,  gave  this  approximated  number  as: 
the  user  cost  =  $45,000/lane/day. 
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3-Coat/Lead 
Do  Nothing  Untill  State  5  and  then  Complete  Repaint 


Age 

(yrs) 

Cost 
initial/ft^ 

Interest  Rate 

25 

$4.00 

7.00% 

Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost 

0 

1 

$39.37 

Present  Value  for  Fv  after  25  yrs 

25 

2.363245 

$17.14 

Present  Value  for  Fv  after  50  yrs 

50 

5.5849269 

$7.46 

Total  Present  Value  TPV/ft^ 
EUAC 

$63.98 
$4.56 

3-CoatyLead 

Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost 

0 

1 

$4.00 

$21.69 

Present  Value  for  Fv  after  15  yrs 

15 

1.675348831 

$1.50 

$11.65 

Present  Value  for  Fv  after  30  yrs 

30 

2.806793705 

$1.50 

$7.07 

Present  Value  for  Fv  after  45  yrs 

45 

4.702358551 

$1.50 

$4.30 

Total  Present  Value  TPV/ft^ 
EUAC 

$44.71 
$3.18 

3-Coat/Lead 

Spot  Repairs  at  State  4 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl.  Factor 

Cost 

Total 

Present  Value  for  initial  Cost 

0 

1 

$4.00 

$30.53 

Present  Value  for  Fv  after  20  yrs 

20 

1.989789 

$2.50 

$14.93 

Present  Value  for  Fv  after  40  yrs 

40 

3.95926 

$2.50 

$7.68 

Total  Present  Value  TPV/ft^ 
EUAC 

$53.13 
$3.78 

3-Coat/Zinc 
Do  Nothing  Until!  State  5  and  then  Complete  Repaint 
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Age 

(yrs) 

Cost 

initial/ft^ 

Interst  Rate 

25 

$2.80 

7.00% 

Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost 

0 

1 

$38.17 

Present  Value  for  Fv  after  25  yrs 

25 

2.363245 

$16.62 

Present  Value  for  Fv  after  50  yrs 

50 

5.584927 

$7.24 

Total  Present  Value  TPV/ft^ 
EUAC 

$62.03 
$4.42 

3-Coat/Zinc 

Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl.  Factor 

Cost 

Total 

Present  Value  for  initial  Cost 

0 

1 

$2.80 

$20.49 

Present  Value  for  Fv  after  15  yrs 

15 

1.675348831 

$1.00 

$11.35 

Present  Value  for  Fv  after  30  yrs 

30 

2.806793705 

$1.00 

$6.89 

Present  Value  for  Fv  after  45  yrs 

45 

4.702358551 

$1.00 

$4.18 

Total  Present  Value  TPV/ft' 
EUAC 

$42.91 
$3.06 

3-Coat/Zinc 

Spot  Repairs  at  State  4 


181 


Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost 

0 

1 

$2.80 

$29.33 

Present  Value  for  Fv  after  20  yrs 

20 

1.989789 

$2.00 

$14.67 

Present  Value  for  Fv  after  40  yrs 

40 

3.95926 

$2.00 

$7.54 

Total  Present  Value  TPV/ft^ 
EUAC 

$51.54 
$3.67 
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3-Coat/Lead 
Do  Nothing  Untill  State  5  and  then  Complete  Repaint 


Age 
(yrs) 

Cost 
initial/ft^ 

Interst  Rate 

25 

$4.00 

7.00% 

Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost: 

0 

1 

$720,000.00 

Present  Value  for  Fv  after  25  yrs: 

25 

2.363245 

$313,506.68 

Present  Value  for  Fv  after  50  yrs: 

50 

5.584927 

$136,508.94 

Total  Present  Value  TPV/ft^ 
EUAC 

$1,170,015.62 
$83,339.31 

3-Coat/Lead 

Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$4.00 

$360,000.00 

Present  Value  for  Fv  after  15  yrs: 

15 

1.675348831 

$1.50 

$218,600.47 

Present  Value  for  Fv  after  30  yrs: 

30 

2.806793705 

$1.50 

$132,739.34 

Present  Value  for  Fv  after  45  yrs: 

45 

4.702358551 

$1.50 

$80,602.45 

Total  Present  Value  TPV/ft^ 
EUAC 

$791,942.26 
$56,409.43 

3-Coat/Lead 

Spot  Repairs  at  State  4 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$4.00 

$540,000.00 

Present  Value  for  Fv  after  20  yrs: 

20 

1.989789 

$2.50 

$277,667.60 

Present  Value  for  Fv  after  40  yrs: 

40 

3.95926 

$2.50 

$142,776.47 

Total  Present  Value  TPV/ft^ 
EUAC 

$960,444.07 
$68,411.69 

3-Coat/Zinc 
Do  Nothing  Untill  State  5  and  then  Complete  Repaint 


185 


Age 

(yrs) 

Cost 

initial/ft^ 

Interst  Rate 

25 

$2.80 

7.00% 

Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Present  Value  for  initial  Cost: 

0 

1 

$720,000.00 

Present  Value  for  Fv  after  25  yrs: 

25 

2.363245 

$313,506.68 

Present  Value  for  Fv  after  50  yrs: 

50 

5.584927 

$136,508.94 

Total  Present  Value  TPV/ft^ 
EUAC 

$1,170,015.62 
$83,339.31 

3-Coat/Zinc 

Spot  Repairs  at  State  3 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$2.80 

$360,000.00 

Present  Value  for  Fv  after  15  yrs: 

15 

1.675348831 

$1.00 

$218,600.47 

Present  Value  for  Fv  after  30  yrs: 

30 

2.806793705 

$1.00 

$132,739.34 

Present  Value  for  Fv  after  45  yrs: 

45 

4.702358551 

$1.00 

$80,602.45 

Total  Present  Value  TPV/ft^                                                 :.      ; 
EUAC 

$791,942.26 
$56,409.43 

3-Coat/Zinc 

Spot  Repairs  at  State  4 
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Interest  Rate 

7.00% 
Inflation  Factor 

0.035 

Item 

n 

Infl. Factor 

Cost 

Total 

Present  Value  for  initial  Cost: 

0 

1 

$2.80 

$540,000.00 

Present  Value  for  Fv  after  20  yrs: 

20 

1.989789 

$2.00 

$277,667.60 

Present  Value  for  Fv  after  40  yrs: 

40 

3.95926 

$2.00 

$142,776.47 

Total  Present  Value  TPV/ft^ 
EUAC 

$960,444.07 
$68,411.69 
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ff  ". 


Figure  A.1 :  PV  ($/ft2)  Comparison  of  3-COAT  System/ZINC  for  Disruption  Cost  Only 
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Figure  A.2:  EUAC  ($/ft2)  for  3-COAT  System/  ZINC  for  Disruption  Cost  Only 
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Figure  A.3:  PV  ($/ft2)  Comparison  of  3  COAT  System  /Lead  for  Disruption  Cost  Only 
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Figure  A.4:  EUAC  ($/ft2)  for  3-COAT  System  /  Lead  for  Disruption  Cost  Only 
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Appendix  E: 
Condition  Rating  Standards 
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